Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 582,152 C→T pseudogene (653/653 nt) ybcY ← pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913582,1520CT100.0% 242.8 / NA 74pseudogene (653/653 nt)ybcYpseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/40);  total (34/40)

CGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582023‑582280
                                                                                                                                  |                                                                                                                                
cGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTc                                                                                                                        >  1:573252/1‑139 (MQ=255)
   gtgtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTg                                                                                                                     <  2:573252/139‑1 (MQ=255)
     gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGaa                                                                                                                   <  2:486316/139‑1 (MQ=255)
     gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGaa                                                                                                                   >  1:593417/1‑139 (MQ=255)
      ttGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAg                                                                                                                  >  2:831882/1‑139 (MQ=255)
           gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  1:831882/139‑1 (MQ=255)
           gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1006108/139‑1 (MQ=255)
           gAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1484474/139‑1 (MQ=255)
                gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  2:1007125/139‑1 (MQ=255)
                gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:453008/139‑1 (MQ=255)
                 tgtTGATACATCAACTGCACCTGATATTGCGTGGCCTACGAACCCTGTCAGGGAGTAATCAGT‑AGGAGTATGCCAGCGCACGTCTTAAGCAAGAACATGCTGCGGGT‑GATGTTATTTTTTTCCTGCACCGGAAAACAtt                                                                                                       >  2:444387/1‑139 (MQ=255)
                 tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       >  2:1100143/1‑139 (MQ=255)
                 tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGCAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       <  1:1145256/139‑1 (MQ=255)
                  gtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                      <  2:1627600/139‑1 (MQ=255)
                           tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 >  2:876144/1‑135 (MQ=255)
                           tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 <  1:876144/135‑1 (MQ=255)
                                     cTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                   <  1:1709785/139‑1 (MQ=255)
                                                  gCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTAAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa                                                                         <  2:67563/136‑1 (MQ=255)
                                                  gCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTAAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa                                                                         >  1:67563/1‑136 (MQ=255)
                                                   ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             >  2:1197664/1‑99 (MQ=255)
                                                   ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  1:1197664/99‑1 (MQ=255)
                                                   ccTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                     >  1:323190/1‑139 (MQ=255)
                                                        gAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                <  1:444387/139‑1 (MQ=255)
                                                             cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           >  2:1486282/1‑123 (MQ=255)
                                                             cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           <  1:1486282/123‑1 (MQ=255)
                                                             cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctc                                                           >  2:1505105/1‑139 (MQ=255)
                                                             cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            >  2:1582529/1‑138 (MQ=255)
                                                             cTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            <  1:1582529/138‑1 (MQ=255)
                                                              tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  2:1048535/139‑1 (MQ=255)
                                                              tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  1:895058/1‑139 (MQ=255)
                                                              tGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  1:1066641/139‑1 (MQ=255)
                                                                  aGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTgaga                                                      <  2:323190/139‑1 (MQ=255)
                                                                   gggAGTAATCATTGGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTgaga                                                      <  1:1630801/139‑1 (MQ=255)
                                                                   gggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGagag                                                     <  2:895058/139‑1 (MQ=255)
                                                                    ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAGGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGt                                                    >  2:141497/1‑139 (MQ=255)
                                                                    ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTAttt                                                                                                                             >  1:1229219/1‑66 (MQ=255)
                                                                    ggAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTAttt                                                                                                                             <  2:1229219/66‑1 (MQ=255)
                                                                      aGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtat                                                  <  1:1505105/139‑1 (MQ=255)
                                                                        tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 <  2:197750/90‑1 (MQ=255)
                                                                        tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                 >  1:197750/1‑90 (MQ=255)
                                                                        tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 >  2:736518/1‑138 (MQ=255)
                                                                        tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  1:736518/138‑1 (MQ=255)
                                                                        tAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTatat                                                >  1:104750/1‑139 (MQ=255)
                                                                          aTCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                >  1:1704290/1‑121 (MQ=255)
                                                                          aTCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaaa                                                                <  2:1704290/121‑1 (MQ=255)
                                                                               t‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTaa                                         >  1:805127/1‑139 (MQ=255)
                                                                                 gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                        >  2:1166496/1‑139 (MQ=255)
                                                                                 gggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAg                                        >  2:345164/1‑139 (MQ=255)
                                                                                  ggATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGc                                       <  2:805127/139‑1 (MQ=255)
                                                                                   gcttATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGcc                                      >  2:831435/3‑139 (MQ=255)
                                                                                     ttATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCctt                                    <  1:141497/139‑1 (MQ=255)
                                                                                        tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  2:1209156/114‑1 (MQ=255)
                                                                                        tGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  1:1209156/1‑114 (MQ=255)
                                                                                            gcagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGa                             <  1:831435/139‑1 (MQ=255)
                                                                                                  cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             <  2:1574772/53‑1 (MQ=38)
                                                                                                  cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                             >  1:1574772/1‑53 (MQ=38)
                                                                                                  cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        >  2:1712286/1‑58 (MQ=255)
                                                                                                  cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:1712286/58‑1 (MQ=255)
                                                                                                   gTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtttt                      >  2:54968/1‑139 (MQ=255)
                                                                                                     cTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           >  2:1355014/1‑52 (MQ=38)
                                                                                                     cTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           <  1:1355014/52‑1 (MQ=38)
                                                                                                              gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaaa           >  2:1400375/1‑139 (MQ=255)
                                                                                                                aCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGa         <  1:1100143/139‑1 (MQ=255)
                                                                                                                    gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCt                                                                                           <  2:622367/53‑1 (MQ=38)
                                                                                                                    gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCt                                                                                           >  1:622367/1‑53 (MQ=38)
                                                                                                                    gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGcc     >  1:210146/1‑139 (MQ=255)
                                                                                                                      tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           <  1:7536/67‑1 (MQ=255)
                                                                                                                      tgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           >  2:7536/1‑67 (MQ=255)
                                                                                                                       gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt  <  1:1400375/139‑1 (MQ=255)
                                                                                                                       gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt  <  1:640921/139‑1 (MQ=255)
                                                                                                                        cGGTT‑GATGTTATTTTTTTCCTGAAGCGG‑AAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 >  1:106624/1‑90 (MQ=255)
                                                                                                                        cGGTT‑GATGTTATTTTTTTCCTGAAGCGG‑AAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  2:106624/90‑1 (MQ=255)
                                                                                                                           gttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGATATGCTCTGAGAGTATATGTGTTAAGcc                                      >  1:1054619/4‑100 (MQ=255)
                                                                                                                           gttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGcc                                      <  2:1054619/97‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                
CGGGTGTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATT‑GGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582023‑582280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑