Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,182 A→C 12.3% M166L (ATG→CTG) ‡ PP_2340 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,673,1820AC12.3% 254.3 / 27.1 97M166L (ATG→CTG) ‡PP_2340hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base C (5/7);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

CTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCT  >  NC_002947/2673038‑2673329
                                                                                                                                                 |                                                                                                                                                   
cTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                                 <  8:279225/149‑1 (MQ=255)
 tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                                 >  4:117448/1‑148 (MQ=255)
   gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                                 <  3:117448/146‑1 (MQ=255)
   gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                               <  4:244399/148‑1 (MQ=255)
        cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                         <  7:342957/149‑1 (MQ=255)
        cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                         <  7:146867/149‑1 (MQ=255)
        cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCAtcgcg                                                                                                                                          >  4:327986/1‑148 (MQ=255)
           tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                        <  3:386822/147‑1 (MQ=255)
           tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                       <  6:50457/148‑1 (MQ=255)
           tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                       <  1:187193/148‑1 (MQ=255)
             gCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaa                                                                                                                                    <  3:327986/149‑3 (MQ=255)
                     ttGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGAACGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGc                                                                                                                            <  1:335939/149‑1 (MQ=255)
                        ctgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCTGCGTTATATCGCCCCCAGCTTCCTGCGCTATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCg                                                                                                                          <  7:259509/148‑1 (MQ=255)
                         tgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCggg                                                                                                                        >  7:387423/1‑149 (MQ=255)
                          gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                             <  6:332014/127‑1 (MQ=255)
                          gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                             >  5:332014/1‑127 (MQ=255)
                          gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                       >  1:184536/1‑149 (MQ=255)
                           ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCg                                                                                                                      <  6:375034/149‑1 (MQ=255)
                                   ccGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTc                                                                                                              >  6:42842/1‑149 (MQ=255)
                                         cTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                                <  6:996/125‑1 (MQ=255)
                                         cTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                                >  5:996/1‑125 (MQ=255)
                                            ggAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                       >  8:331667/1‑115 (MQ=255)
                                             gAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                    >  7:241678/1‑149 (MQ=255)
                                             gAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                    >  3:131471/1‑149 (MQ=255)
                                              aaCTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCATGAGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCCGTCGCCGGTGTAtt                                                                                                   >  8:167035/1‑149 (MQ=255)
                                                   ccccAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                              >  2:19023/1‑149 (MQ=255)
                                                   ccccAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                              <  8:387423/149‑1 (MQ=255)
                                                       aaTTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGc                                                                                          >  6:409110/1‑149 (MQ=255)
                                                        aactcgcgATCGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCt                                                                                         >  1:154276/5‑149 (MQ=255)
                                                            cgcgATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTg                                                                                     >  5:171040/1‑149 (MQ=255)
                                                             gcgatgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                       <  5:2771/114‑1 (MQ=255)
                                                             gcgatgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                       >  6:2771/1‑114 (MQ=255)
                                                               gatgcgatCGATCGCGAGTTCTTCCTGCTTTCCGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg                                                                                  <  3:242327/149‑1 (MQ=255)
                                                                 tgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                                   <  3:62726/114‑1 (MQ=255)
                                                                 tgcgatCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                                   >  4:62726/1‑114 (MQ=255)
                                                                    gatcgatcGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTc                                                                             >  2:279175/1‑149 (MQ=255)
                                                                     atcgatcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCg                                                                            >  6:340809/1‑149 (MQ=255)
                                                                      tcgatcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGCCAGCGCGCAGCGCCAGTGAGCCGGGCGAATCAGTCGCCCGTGTAt                                                                                                    >  6:210777/1‑124 (MQ=255)
                                                                      tcgatcgCGAGTGCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                    <  5:210777/124‑1 (MQ=255)
                                                                           tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCTGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                       <  8:241678/149‑1 (MQ=255)
                                                                           tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                       >  6:278358/1‑149 (MQ=255)
                                                                           tcgCGAGTTCATCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                       <  4:29490/149‑1 (MQ=255)
                                                                               gAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCg                                                                    >  6:112593/1‑148 (MQ=255)
                                                                                   ccttcCTGCTTTCAGTGTCGCACTACCAGCGCTATATCGCCCCCAGGTTCCTGCGCGATGCCATGGGCAGCGCGGAGGGCCAGTGAGGCGGGCGAATCATTCGCCGGTGTATTCGCAACCGTTGGTGCAGGTTTCGTGGATTCGAACCt                                                               <  1:219745/148‑1 (MQ=255)
                                                                                    cttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCt                                                               <  8:325667/148‑1 (MQ=255)
                                                                                     ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                             >  3:292433/1‑149 (MQ=255)
                                                                                     ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                             >  1:350120/1‑149 (MQ=255)
                                                                                     ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                             >  6:367394/1‑149 (MQ=255)
                                                                                     ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                                   <  5:422798/79‑4 (MQ=21)
                                                                                     ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                                   >  6:422798/1‑76 (MQ=21)
                                                                                       cttGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCCGCTTCCTGCGCTATGCCTTGGGCAGGGCGCAGGGCCAGTGAGCCTGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                              <  2:184536/144‑1 (MQ=255)
                                                                                       ccTGGTTTCAGAGTCGCACTACCAGCGGTATATCGCCCCGAGCTTCCTGCGCTTAGCCATGGGCAGCGCGAAGGGCCAGTGAGCGGTGCGAATCAGTCGCCGGTGTATTCGCAACCGGTGGTGCAGGTTTCGt                                                                           <  7:232913/133‑1 (MQ=255)
                                                                                       ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                              <  4:131471/146‑1 (MQ=255)
                                                                                       ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                             <  4:292433/147‑1 (MQ=255)
                                                                                       ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                             <  1:19023/147‑1 (MQ=255)
                                                                                       ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTACTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                           >  8:232913/1‑133 (MQ=255)
                                                                                        cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAc                                                          <  4:380811/149‑1 (MQ=255)
                                                                                         tGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                           >  3:349276/1‑147 (MQ=255)
                                                                                           cTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGc                                                       >  4:127750/1‑149 (MQ=255)
                                                                                                 gagGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                                 <  6:171040/149‑1 (MQ=255)
                                                                                                 gagGCGCACTACCAGCGCTATATCGCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                                >  7:387922/1‑149 (MQ=255)
                                                                                                  agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                                 <  2:350120/148‑1 (MQ=255)
                                                                                                  agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                                >  6:2643/1‑149 (MQ=255)
                                                                                                   ggCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCAt                                               >  7:195906/1‑149 (MQ=255)
                                                                                                    gcgcACTACCAGCGATATATCCCCCCCCGCGTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATc                                              <  2:349671/149‑1 (MQ=255)
                                                                                                     cgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                             <  8:195906/149‑1 (MQ=255)
                                                                                                      gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAg                                            >  1:222671/1‑149 (MQ=255)
                                                                                                            ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGc                                        <  1:279175/147‑1 (MQ=255)
                                                                                                            ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                      <  3:127750/149‑1 (MQ=255)
                                                                                                              aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                                    >  4:17232/1‑149 (MQ=255)
                                                                                                               gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                         >  5:240185/1‑127 (MQ=255)
                                                                                                               gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGAc                                   <  5:340809/149‑1 (MQ=255)
                                                                                                                   tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTgg                               >  2:175996/1‑149 (MQ=255)
                                                                                                                     tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                      <  3:32912/140‑1 (MQ=255)
                                                                                                                     tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                      >  4:32912/1‑140 (MQ=255)
                                                                                                                       tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                            >  3:212559/1‑148 (MQ=255)
                                                                                                                         gCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                            >  2:141569/1‑146 (MQ=255)
                                                                                                                          cccccAGCTTCCTGCGCGTTGCCATGGGAAGCGCGCAGGGCCAGTGAGCCCGGCGAATCAGTCCCCGGTGTATTAGCCACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTGCGGCATCAGCGGCTTGACCTGGTcc                            >  5:84316/1‑145 (MQ=255)
                                                                                                                          cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                        <  5:278358/149‑1 (MQ=255)
                                                                                                                          cccccAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                        <  5:367394/149‑1 (MQ=255)
                                                                                                                           ccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGAt                       <  1:175996/149‑1 (MQ=255)
                                                                                                                             ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                        <  2:142972/82‑1 (MQ=255)
                                                                                                                             ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                        >  1:142972/1‑82 (MQ=255)
                                                                                                                               aGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCAc                   <  5:409110/149‑1 (MQ=255)
                                                                                                                                   cccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTgg               >  2:61501/2‑149 (MQ=255)
                                                                                                                                     cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                  >  2:232139/1‑64 (MQ=21)
                                                                                                                                     cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                  <  1:232139/64‑1 (MQ=21)
                                                                                                                                     cTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc              <  4:212559/148‑1 (MQ=255)
                                                                                                                                      tGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGa            >  8:46422/1‑149 (MQ=255)
                                                                                                                                        cgcgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt           <  5:42842/148‑1 (MQ=255)
                                                                                                                                          cgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCAc        <  7:46422/149‑1 (MQ=255)
                                                                                                                                           gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc          >  5:251632/1‑146 (MQ=255)
                                                                                                                                           gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc          <  6:251632/146‑1 (MQ=255)
                                                                                                                                               ccATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGc   >  3:367669/1‑149 (MQ=255)
                                                                                                                                               ccATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGc   >  5:285927/1‑149 (MQ=255)
                                                                                                                                                cATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCt  <  5:270363/149‑1 (MQ=255)
                                                                                                                                                cATGGGCAGCGAGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAAAGCTCCGGCATCATCGGCTTTACCTGGTCCCAGCTACACTTGGCGATCAATTctct  >  6:326193/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCGATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCT  >  NC_002947/2673038‑2673329

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: