Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,254,518 T→G 5.4% P279P (CCT→CCG)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,254,5180TG 5.4% 235.3 / 4.4 74P279P (CCT→CCG) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base G (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

GTGCGTGCCCCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGAC  >  NC_002947/2254372‑2254665
                                                                                                                                                  |                                                                                                                                                   
gtgcgtgcCCCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCtt                                                                                                                                                    <  8:379803/148‑1 (MQ=255)
gtgcgtgcCCCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCt                                                                                                                                                     <  8:200804/147‑1 (MQ=255)
         cccAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTc                                                                                                                                          >  2:199112/1‑149 (MQ=255)
              aCAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCGTGGCAAGGCTCgtaca                                                                                                                                     <  8:347628/149‑5 (MQ=255)
                aGTTCGACGAGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGATCCTTTCCACGGCTCGGCACCg                                                                                                                                   >  5:233028/1‑149 (MQ=255)
                  ttCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGa                                                                                                                                 >  5:407961/1‑149 (MQ=255)
                    cGACGTGATGGTGACCGCCAACATGTTCGGGGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATTCTGCCTTCGGCGTCCCTGGATGCGGACAACATCGGCATGTACGAGCCTTGCCACGGCTGGGCACCGGAc                                                                                                                                <  6:233028/148‑1 (MQ=255)
                      aCGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCGTGGCAAGGCTCgtaca                                                                                                                                     >  6:72549/1‑137 (MQ=255)
                      aCGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCGTGGCAAGGCTCgtaca                                                                                                                                     <  5:72549/141‑5 (MQ=255)
                            tGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGc                                                                                                                       >  4:161364/1‑149 (MQ=255)
                            tGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGc                                                                                                                       <  3:186494/149‑1 (MQ=255)
                             ggTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCCGCATCCTGCCTACCGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGcc                                                                                                                      >  2:107373/1‑149 (MQ=255)
                               tGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCtt                                                                                                                    >  7:340909/1‑149 (MQ=255)
                               tGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCtt                                                                                                                    >  3:334065/1‑149 (MQ=255)
                                 aCCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTg                                                                                                                   <  6:407961/148‑1 (MQ=255)
                                   cGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCa                                                                                                                >  1:118785/1‑149 (MQ=255)
                                      caacaTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGTCCTTGGCATCg                                                                                                             >  8:158248/1‑149 (MQ=255)
                                      caacaTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCg                                                                                                             <  4:389550/149‑1 (MQ=255)
                                      caacaTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCg                                                                                                             <  5:196575/149‑1 (MQ=255)
                                          aTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCaa                                                                                                         >  4:107375/1‑149 (MQ=255)
                                                          tGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTg                                                                                                  >  1:163419/1‑140 (MQ=255)
                                                          tGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTg                                                                                                  <  2:163419/140‑1 (MQ=255)
                                                               gATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTc                                                                                    >  4:122104/1‑149 (MQ=255)
                                                               gATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATAGCCAACCCACTGGCGACCATCCTGTc                                                                                    <  1:216743/149‑1 (MQ=255)
                                                                  gAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGgt                                                                                 >  1:280738/1‑149 (MQ=255)
                                                                    aGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGgtgt                                                                               >  8:365857/1‑149 (MQ=255)
                                                                            tGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGatgatg                                                                       >  2:199559/1‑149 (MQ=255)
                                                                                gACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                   >  5:206651/1‑149 (MQ=255)
                                                                                gACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                   >  3:66536/1‑149 (MQ=255)
                                                                                gACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                   >  4:363362/1‑149 (MQ=255)
                                                                                 aCCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                   <  3:363362/148‑1 (MQ=255)
                                                                                 aCCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                   <  1:249994/148‑1 (MQ=255)
                                                                                 aCCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGAGGTTGCg                                                                  >  3:321772/1‑149 (MQ=255)
                                                                                   cGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGt                                                                 >  4:42160/1‑148 (MQ=255)
                                                                                   cGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGACTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGt                                                                 >  8:27476/1‑148 (MQ=255)
                                                                                      ttCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACa                                                             <  7:367038/149‑1 (MQ=255)
                                                                                          aTCGGCATGCTGCCTTCGGCGTCGCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCtt                                                         <  2:210307/149‑1 (MQ=255)
                                                                                          aTCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCtt                                                         >  1:189506/1‑149 (MQ=255)
                                                                                             ggCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCaa                                                      <  3:161364/149‑1 (MQ=255)
                                                                                             ggCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCaa                                                      <  8:340909/149‑1 (MQ=255)
                                                                                                  gctgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGt                                                 >  1:174909/1‑149 (MQ=255)
                                                                                                   ctgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGtt                                                                     <  7:242350/128‑1 (MQ=255)
                                                                                                   ctgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGtt                                                                     >  8:242350/1‑128 (MQ=255)
                                                                                                   ctgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGt                                                 <  7:27476/148‑1 (MQ=255)
                                                                                                   ctgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTc                                                <  3:122104/149‑1 (MQ=255)
                                                                                                    tgcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCg                                               >  5:334341/1‑149 (MQ=255)
                                                                                                     gcCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCgg                                              <  4:139801/149‑1 (MQ=255)
                                                                                                       cTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCt                                            <  1:107373/149‑1 (MQ=255)
                                                                                                        ttCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTg                                           >  5:74159/1‑149 (MQ=255)
                                                                                                        ttCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTg                                           >  2:311871/1‑149 (MQ=255)
                                                                                                        ttCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTg                                           >  4:170357/1‑149 (MQ=255)
                                                                                                               tCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGa                                    <  2:97328/149‑1 (MQ=255)
                                                                                                                 ccTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAgg                                  <  6:334341/149‑1 (MQ=255)
                                                                                                                     gATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGa                               <  1:311871/148‑1 (MQ=255)
                                                                                                                     gATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGa                               <  1:199559/148‑1 (MQ=255)
                                                                                                                       ggCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCg                            <  8:322705/148‑1 (MQ=255)
                                                                                                                         ccGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCga                           <  8:233561/148‑1 (MQ=255)
                                                                                                                          cGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCgag                          >  3:386467/1‑148 (MQ=255)
                                                                                                                             gtggcaaGGCTCGTACATGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAgg                      >  7:143445/6‑149 (MQ=255)
                                                                                                                              aacaaAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGc                     >  3:79957/1‑149 (MQ=255)
                                                                                                                                caaAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATc                             >  4:53953/1‑139 (MQ=255)
                                                                                                                                caaAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATc                             <  3:53953/139‑1 (MQ=255)
                                                                                                                                 aaaGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGt                  <  3:170357/149‑1 (MQ=255)
                                                                                                                                  aaGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTg                 <  2:174909/149‑1 (MQ=255)
                                                                                                                                     gCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGc              >  1:262551/1‑149 (MQ=255)
                                                                                                                                      cATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGcc             <  7:365857/149‑1 (MQ=255)
                                                                                                                                       aTGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGcc                    >  6:7905/1‑141 (MQ=255)
                                                                                                                                       aTGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGcc                    <  5:7905/141‑1 (MQ=255)
                                                                                                                                       aTGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCt            <  7:158248/149‑1 (MQ=255)
                                                                                                                                           aCGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGcc                    <  2:96953/137‑1 (MQ=255)
                                                                                                                                           aCGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGcc                    >  1:96953/1‑137 (MQ=255)
                                                                                                                                           aCGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGtt        >  2:232554/1‑149 (MQ=255)
                                                                                                                                            cGAGCCGTGGCAAGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTc       >  3:397531/1‑149 (MQ=255)
                                                                                                                                             gAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCt      >  6:205428/1‑149 (MQ=255)
                                                                                                                                                ccTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGa   >  3:289767/1‑149 (MQ=255)
                                                                                                                                                ccTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGa   >  5:335284/1‑149 (MQ=255)
                                                                                                                                                  ttGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGAc  >  5:102578/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGCGTGCCCCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGAC  >  NC_002947/2254372‑2254665

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