Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,132 (C)8→7 pseudogene (3955/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,1250C.96.8% 250.9 / 2.2 62pseudogene (3962/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (25/35);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCC  >  NC_000913/1593988‑1594253
                                                                                                                                         |                                                                                                                                
cGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccc                                                                                                                                <  1:288842/140‑1 (MQ=255)
   catcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccc                                                                                                                             <  1:360569/140‑1 (MQ=255)
      caGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccag                                                                                                                          >  2:1266429/1‑138 (MQ=255)
        ggCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                                                        >  1:1390693/1‑136 (MQ=255)
         gCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCa                                                                                                                      <  2:1309719/140‑1 (MQ=255)
          ccATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACGTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCAt                                                                                                                     <  1:429683/140‑1 (MQ=255)
          ccATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCAt                                                                                                                     <  2:853567/140‑1 (MQ=255)
          ccATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCAt                                                                                                                     >  2:278129/1‑140 (MQ=255)
           cATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCTGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATa                                                                                                                    <  2:1262350/140‑1 (MQ=255)
           cATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATa                                                                                                                    >  1:78789/1‑140 (MQ=255)
             ttGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAAc                                                                                                                  >  2:922102/1‑140 (MQ=255)
              tGTTTCCGGTATAGTTGCCGTTAACGGTAAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACg                                                                                                                 <  1:1166467/140‑1 (MQ=255)
               gTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGt                                                                                                                <  2:739908/140‑1 (MQ=255)
                      gTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATc                                                                                                         <  1:868021/140‑1 (MQ=255)
                          aGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCt                                                                                                     >  2:699388/1‑140 (MQ=255)
                              gCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                                                                 >  2:305524/1‑140 (MQ=255)
                                    aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCat                                                                                                            <  2:369155/124‑1 (MQ=255)
                                    aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCat                                                                                                            >  1:369155/1‑124 (MQ=255)
                                    aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGa                                                                                           <  1:734342/140‑1 (MQ=255)
                                      cGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACa                                                                                         >  2:1305099/1‑140 (MQ=255)
                                          cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                                                                     >  1:574838/1‑140 (MQ=255)
                                           aatACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                                                    <  1:305524/138‑1 (MQ=255)
                                           aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                                                    <  1:1005751/140‑1 (MQ=255)
                                                aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                                                               <  1:699388/140‑1 (MQ=255)
                                                 ttCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                                              >  2:103891/1‑140 (MQ=255)
                                                  tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                                             >  1:198000/1‑140 (MQ=255)
                                                  tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                                             >  1:84470/1‑140 (MQ=255)
                                                  tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                                             >  1:1361019/1‑140 (MQ=255)
                                                     ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                                                          <  2:190402/140‑1 (MQ=255)
                                                      cGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                         <  1:1266429/140‑1 (MQ=255)
                                                       ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                                                        >  2:613848/1‑140 (MQ=255)
                                                       ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                       >  1:739453/1‑140 (MQ=255)
                                                        gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                       <  2:883239/140‑1 (MQ=255)
                                                        gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                       <  1:1305099/140‑1 (MQ=255)
                                                        gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                       <  1:1403469/140‑1 (MQ=255)
                                                           tCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                                                                    >  2:1074455/1‑140 (MQ=255)
                                                            ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                         >  2:892654/1‑134 (MQ=255)
                                                            ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                         <  1:892654/134‑1 (MQ=255)
                                                             cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                                                  <  1:882532/140‑1 (MQ=255)
                                                              ccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                                                                 <  2:198000/140‑1 (MQ=255)
                                                               cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                                                                <  1:103891/140‑1 (MQ=255)
                                                                   tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGat                                                            >  2:40106/1‑140 (MQ=255)
                                                                     gCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                                                         >  1:387527/1‑140 (MQ=255)
                                                                       tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTa                                                        <  2:574838/140‑1 (MQ=255)
                                                                              ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt                                                 >  1:826263/1‑140 (MQ=255)
                                                                              ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGtt                                                 <  2:527707/140‑1 (MQ=255)
                                                                                 aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                                              <  2:1274709/140‑1 (MQ=255)
                                                                                    cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg                                           <  1:1253162/140‑1 (MQ=255)
                                                                                      tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                         >  2:1410435/1‑140 (MQ=255)
                                                                                      tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                         <  2:1361019/140‑1 (MQ=255)
                                                                                      tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                         <  1:201145/140‑1 (MQ=255)
                                                                                      tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                         >  2:1237196/1‑140 (MQ=255)
                                                                                       gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCa                                        <  2:1279093/140‑1 (MQ=255)
                                                                                                 tcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                              <  2:214964/140‑1 (MQ=255)
                                                                                                  ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCag                            <  2:387527/140‑1 (MQ=255)
                                                                                                     ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCagag                          <  1:613848/140‑1 (MQ=255)
                                                                                                          ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGa                     <  2:1390693/140‑1 (MQ=255)
                                                                                                            gccgATAGTGGTTTTAGCGACACGCAGTG‑TCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAAc                   <  1:1005338/140‑1 (MQ=255)
                                                                                                               gATAGTGGTTTTAGCGACACGCAGTG‑CACCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                              >  2:466964/1‑126 (MQ=255)
                                                                                                               gATAGTGGTTTTAGCGACACGCAGTG‑CACCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                              <  1:466964/126‑1 (MQ=255)
                                                                                                                           aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACGTCACACCGCTGTTATGGATAtt                                                         >  2:964970/1‑87 (MQ=255)
                                                                                                                           aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                                                         <  1:964970/87‑1 (MQ=255)
                                                                                                                           aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc    >  2:1150334/1‑140 (MQ=255)
                                                                                                                           aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc    <  2:1084521/140‑1 (MQ=255)
                                                                                                                             cGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc  >  1:584397/1‑140 (MQ=255)
                                                                                                                                 acGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc     <  2:1172796/133‑1 (MQ=255)
                                                                                                                                 acGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGACGGAACCGAAAGCAATTTcc     >  1:1172796/1‑133 (MQ=255)
                                                                                                                                         |                                                                                                                                
CGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCC  >  NC_000913/1593988‑1594253

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