Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,254,524 C→A 5.1% H281Q (CAC→CAA leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,254,5240CA 5.1% 256.4 / 3.2 79H281Q (CAC→CAAleuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base A (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

CCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCAGGGG  >  NC_002947/2254381‑2254671
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                             <  3:186494/149‑1 (MQ=255)
                   tGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGc                                                                                                                             >  4:161364/1‑149 (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                                                                                                                   <  5:196575/149‑1 (MQ=255)
                             caacaTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCg                                                                                                                   <  4:389550/149‑1 (MQ=255)
                                 aTGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCaa                                                                                                               >  4:107375/1‑149 (MQ=255)
                                                 tGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTg                                                                                                        <  2:163419/140‑1 (MQ=255)
                                                 tGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTg                                                                                                        >  1:163419/1‑140 (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                                                                         >  4:363362/1‑149 (MQ=255)
                                                                       gACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                         >  3:66536/1‑149 (MQ=255)
                                                                        aCCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                         <  1:249994/148‑1 (MQ=255)
                                                                        aCCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGc                                                                         <  3:363362/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                                                            <  8:340909/149‑1 (MQ=255)
                                                                                    ggCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCaa                                                            <  3:161364/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                                                 >  4:170357/1‑149 (MQ=255)
                                                                                               ttCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTg                                                 >  2:311871/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         >  5:335284/1‑149 (MQ=255)
                                                                                                                                       ccTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGa         >  3:289767/1‑149 (MQ=255)
                                                                                                                                         ttGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGAc        >  5:102578/1‑148 (MQ=255)
                                                                                                                                           gCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCa      <  6:74159/148‑1 (MQ=255)
                                                                                                                                            ccACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCAgg    <  3:42160/149‑1 (MQ=255)
                                                                                                                                            ccACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCAgg    <  3:302370/149‑1 (MQ=255)
                                                                                                                                              aCGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCAgggg  <  2:280738/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCCAAACAGTTCGACGTGATGGTGACCGACAACATGTTCGGTGACATTCTGTCGGATGAAGCTTCCATGTTGACCGGTTCCATCGGCATGCTGCCTTCGGCGTCCCTGGATGCCGACAACAAAGGCATGTACGAGCCTTGCCACGGCTCGGCACCGGACATCGCCGGCCTTGGCATCGCCAACCCACTGGCGACCATCCTGTCGGTGTCGATGATGTTGCGTTACAGCTTCAATCAGTCGGCTGCCGCCGAGGCGATCGAGAAGGCCGTGAGCCTGGTTCTGGACCAGGGG  >  NC_002947/2254381‑2254671

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: