Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,956 T→C 5.5% intergenic (+354/+260) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9560TC 5.5% 563.6 / 14.5 202intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (54/126);  new base C (6/4);  total (67/135)
Fisher's exact test for biased strand distribution p-value = 7.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 205 total aligned reads displayed.

CGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGG  >  NC_000913/3269807‑3270105
                                                                                                                                                     |                                                                                                                                                     
     cATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                                  <  1:2803325/150‑1 (MQ=255)
     cATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                                  <  2:696251/150‑1 (MQ=255)
      aTCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                                 <  2:3624835/150‑1 (MQ=255)
          ggCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                             <  1:3656025/150‑1 (MQ=255)
          ggCATGAATAAAGAGCACTCAACAATCTGAAACCCGCCGGTACTGTTTATTACGTACATCATGTAAAACTGACCGATATGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                             <  2:3951574/150‑1 (MQ=21)
                actAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                       <  2:2104380/148‑1 (MQ=255)
                aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                       <  2:3294568/150‑1 (MQ=255)
                  tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCAGCTTTCTTTTGGGTATATTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGTAAAAAACCCCCCAAGAgg                                                                                                                                     >  1:3235524/1‑150 (MQ=18)
                      gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                 <  2:1881843/150‑1 (MQ=255)
                      gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                 <  2:1733737/150‑1 (MQ=255)
                      gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCACCCCCAACAAAAGAAAGAACCAGCCGAAGCGGGtgg                                                                                                                                 >  2:776495/1‑148 (MQ=18)
                        gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                               <  2:2747007/150‑1 (MQ=255)
                        gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                               >  2:258900/1‑150 (MQ=255)
                        gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                               <  1:174038/150‑1 (MQ=255)
                        gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                               <  1:894548/150‑1 (MQ=255)
                         cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                              <  2:1070664/150‑1 (MQ=255)
                          aCTCAACAATCTGAAACCCGCCGGAACGTTTTATTACGTACATCAGTAAAAACTGACCGATACGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                             <  1:1385921/150‑1 (MQ=18)
                           cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTAct                                                                                                                            <  1:2089270/150‑2 (MQ=255)
                              aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                                         >  2:3805200/1‑150 (MQ=255)
                                              ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                         <  2:2836434/150‑1 (MQ=255)
                                              ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                         <  1:2564564/150‑1 (MQ=255)
                                                 gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCAATCATCATTCTGCCCCTCCAAAAGTAAAAACCCCCCGAAGCCGGGTTTTACGGCCAACAGGTGCAACGGAccca                                                                                                      >  1:3119332/1‑148 (MQ=17)
                                                  aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAGAACAGGTGAAACTGACCGAt                                                                                                     <  1:2133918/150‑1 (MQ=255)
                                                   aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                    <  2:1181428/150‑1 (MQ=255)
                                                    cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTCATCTTTCTGCCCCTCCCAAATTAAAAACCCGCCGAAGCGGGTTTTTCAGTCCAACAGGTGAAACGGACCGGTac                                                                                                   >  1:2874493/1‑149 (MQ=14)
                                                    cGGTTTATTACCTACATCAGGTAAAACAGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCATTAAAATTACAGCCCCTCCAAAAGTAAAAAACCCCCCAACCGGGTTTTTTAGTAAAAAAGCTTAAACAGAcccaaaa                                                                                                   >  2:3790582/1‑145 (MQ=9)
                                                         tattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                                                                                              <  1:647201/150‑1 (MQ=255)
                                                           ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt                                                                                            >  1:2610198/1‑150 (MQ=255)
                                                           ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCCATCATATTTAAGCACCAACAAAAGTAAAAACCCGCCGAAGAGGGTTTTTTAGTAAAACAGGTGAAACAGACCGATAAGACGcgt                                                                                            >  2:1472929/1‑148 (MQ=14)
                                                             aCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                                                                                                                                 >  2:2751728/1‑95 (MQ=37)
                                                              cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAATAACCCGCCGAAGCGGGTTGTTACGTAAGTCAGGTGAAACTGACCGATGAGCCGcgttct                                                                                         >  2:3359152/1‑150 (MQ=14)
                                                              cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcgttct                                                                                         >  1:2066898/1‑150 (MQ=255)
                                                               gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAAAGGTGAAACTGACCGATAAGCCGctttctt                                                                                        >  2:429714/1‑150 (MQ=255)
                                                                tACATCATGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                       <  2:78170/150‑1 (MQ=255)
                                                                  cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                     <  2:2955727/150‑1 (MQ=39)
                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                    <  1:3312873/150‑1 (MQ=38)
                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                    <  2:1236377/150‑1 (MQ=38)
                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                    <  1:1857574/150‑1 (MQ=38)
                                                                   aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                    <  1:739241/150‑1 (MQ=38)
                                                                        gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATAGGGTTTTTACGGAAAACAGGTGAAACAGACCGATAAGACGCTGTCTTTTTGTtata                                                                               >  2:1404932/1‑150 (MQ=11)
                                                                         tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                              <  2:4029599/150‑1 (MQ=35)
                                                                         tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATAAACTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAg                                                                              >  1:3358291/1‑150 (MQ=255)
                                                                               tGACCGATAATCAGCTTTCTTTTGTGTATAGTGTCGCGGACAGTCTGTCATCTTACTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                        <  2:871613/150‑1 (MQ=11)
                                                                               tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                        <  1:1002978/150‑1 (MQ=37)
                                                                                gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                       <  1:3805717/150‑1 (MQ=37)
                                                                                gACCGATAAGCCGCTTTCGTTTGGGTATAGTGTCGTGGACAGTCATTCAACTTTATGCCCCACCAAAAATAAAAAACCGCCGAAGCGGGTGTTTACGTAAAAAAGGTGTAACAGACCGATAAGCCCCTTTCGTTTGGGGATAGTGTCGTg                                                                       >  1:723345/1‑150 (MQ=11)
                                                                                  ccGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTCATTCATCTTTCCGCCCCTCCAACAGTAAAACCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGTAACTGACCGATAAGCAGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                     <  2:3825579/150‑1 (MQ=11)
                                                                                   cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                                                                                             >  2:2421053/1‑125 (MQ=35)
                                                                                          ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                             <  1:36174/150‑1 (MQ=37)
                                                                                          ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCCCCTACAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGTGGACAGtcaat                                                             >  2:2177836/1‑148 (MQ=255)
                                                                                          ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                             <  1:1018891/150‑1 (MQ=37)
                                                                                           cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                            <  1:3369937/150‑1 (MQ=37)
                                                                                                 ctttTGGGTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                      <  2:1763738/150‑1 (MQ=37)
                                                                                                 ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                      <  1:2676134/150‑1 (MQ=37)
                                                                                                  ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTAAGTAAAACAGGTGAAACAGACCGATAAGCCGCGTTCTTTTGGGGATAGTGTAGTGTACAGtcattcaacatga                                                     >  1:3040653/1‑144 (MQ=11)
                                                                                                  atttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTc                                                     <  1:1702371/149‑1 (MQ=37)
                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTg                                                   <  2:2227667/150‑1 (MQ=37)
                                                                                                     tGGGTTTAGTGTCGTCGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                  <  2:2308645/150‑1 (MQ=255)
                                                                                                     tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTTCGTATAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTTTCGTGGACAGTCATTCATCTTTCTGc                                                  <  1:3790582/150‑1 (MQ=255)
                                                                                                      gggTATAGTGTCGTGGACAGTCCTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                 <  2:1317914/150‑1 (MQ=37)
                                                                                                        gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                               <  2:2081897/150‑1 (MQ=37)
                                                                                                        gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                               >  2:558319/1‑150 (MQ=37)
                                                                                                          ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                     <  1:3396970/62‑1 (MQ=35)
                                                                                                          ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                     >  2:3396970/1‑62 (MQ=35)
                                                                                                          ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTTTCTGCCCCTc                                             <  2:2518838/150‑1 (MQ=37)
                                                                                                              tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaa                                         <  2:2610198/150‑1 (MQ=37)
                                                                                                               gtCGTGGACAGTCATTCTTCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                        <  1:2879796/150‑1 (MQ=255)
                                                                                                               gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                        <  1:2703860/150‑1 (MQ=37)
                                                                                                                  gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTGGTGGCAAGTCAATCAAGGTTCAGCCCCTCCAAAAGcc                                     >  2:2380934/1‑149 (MQ=11)
                                                                                                                   tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa                                    <  2:2017033/150‑1 (MQ=37)
                                                                                                                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaaa                                 <  1:2199411/150‑1 (MQ=37)
                                                                                                                      aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaaa                                 <  2:3003797/150‑1 (MQ=37)
                                                                                                                         gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                              <  1:3152673/150‑1 (MQ=37)
                                                                                                                         ctcattcatCTATCTCCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGAATTTACGTAAAACAGCTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAACCAAAAAccc                              <  1:111669/149‑1 (MQ=11)
                                                                                                                          tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                             >  1:3940635/1‑150 (MQ=37)
                                                                                                                          tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                             >  1:2376555/1‑150 (MQ=37)
                                                                                                                               catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACATTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                        <  1:2121028/150‑1 (MQ=255)
                                                                                                                                  cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                >  1:271832/1‑123 (MQ=17)
                                                                                                                                  cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                <  2:271832/123‑1 (MQ=17)
                                                                                                                                   tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                    <  1:4018700/150‑1 (MQ=37)
                                                                                                                                   tttCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGGTTTTAAGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTATTGTGGGTATAGTGTCGTGGACAGTCATACAACTTTCTGCCCCTCCAAAAGCACAAACACGCCGAAGCgg                    >  1:3554718/1‑150 (MQ=11)
                                                                                                                                   ttccggCCCCGCCCACAGAACAAACCCGCCGCAGCGGGTTTTTACGTCAAACCGCTGAAACGGACCGATAAGCCGCTGTCTTTCGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                    <  1:2325230/145‑1 (MQ=11)
                                                                                                                                     tCTGCCCCTCCACCAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                  <  2:3235248/150‑1 (MQ=255)
                                                                                                                                       tGCCCCTCCACACGGACTAACGCGCCGCAGCGGGTTTTTACGTAACACCGGTGAAACTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                <  2:3119332/150‑1 (MQ=11)
                                                                                                                                       tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                >  2:3061932/1‑150 (MQ=37)
                                                                                                                                        gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt               >  1:324219/1‑150 (MQ=37)
                                                                                                                                        gCCCCACCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttt               <  1:1125518/150‑1 (MQ=11)
                                                                                                                                         ccccTCCAAAAGTAAAAACCCGCCGACGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt              <  1:1827982/150‑1 (MQ=255)
                                                                                                                                          cccTCCACAAGTAAACACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa             <  1:592000/150‑1 (MQ=12)
                                                                                                                                          ccatCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGAGGACAGTCATTCATAATAGAGCCCCACCAAAATCAAAAACCCGCCGCAGCCGGttagta             >  2:3532045/4‑146 (MQ=11)
                                                                                                                                           ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCGTTTGGGGCTAGTGTCGTGGACAGTCAATCATCTTTCTGCCCCGCCAAAACCAAAAAACCGCCCAAGCGCGTTTTTAc            >  1:569627/1‑150 (MQ=11)
                                                                                                                                             tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                                                                                                >  2:1396177/1‑48 (MQ=255)
                                                                                                                                               cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGATAGTCATTCATCATTCTTCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGGTTTTACGTat        >  2:3780692/1‑149 (MQ=11)
                                                                                                                                               cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa        >  1:1285525/1‑150 (MQ=37)
                                                                                                                                               cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAAACCGCCGAACCGGGTTTTTACGTaa        >  1:412146/1‑150 (MQ=255)
                                                                                                                                                aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa       <  1:1234713/150‑1 (MQ=37)
                                                                                                                                                aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa       >  1:2415413/1‑150 (MQ=37)
                                                                                                                                                aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa       >  1:221768/1‑150 (MQ=37)
                                                                                                                                                 aaaGCACAAACCCGCCGCAGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc      <  2:1608695/150‑1 (MQ=255)
                                                                                                                                                    gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg   <  1:3567795/150‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                     
CGCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGG  >  NC_000913/3269807‑3270105

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