Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 339,228 C→G 6.3% M12I (ATG→ATC)  PP_0281 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947339,2280CG 6.3% 297.3 / 8.9 96M12I (ATG→ATC) PP_0281amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (40/50);  new base G (3/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GCCGCGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGCTTT  >  NC_002947/339081‑339377
                                                                                                                                                   |                                                                                                                                                      
gccgcGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCcatcat                                                                                                                                                      <  2:168966/150‑1 (MQ=255)
  cgcgAACCAGGGTCGAGTAGAAGCCGCCAAGTCATTGCAGCGGCCTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATgg                                                                                                                                                    >  2:1486562/1‑150 (MQ=255)
  cgcgAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATgg                                                                                                                                                    <  1:10320/150‑1 (MQ=255)
  cgcgAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATgg                                                                                                                                                    >  1:2104854/1‑150 (MQ=255)
      aaCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGCGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATGATGGCCgg                                                                                                                                                <  1:1730077/150‑1 (MQ=255)
      aaCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATGATGGCCgg                                                                                                                                                <  1:1841295/150‑1 (MQ=255)
      aaCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCgg                                                                                                                                                >  1:1082687/1‑150 (MQ=255)
            ggTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGaa                                                                                                                                         >  1:1296763/1‑150 (MQ=255)
             gTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAAt                                                                                                                                        <  1:1506120/150‑1 (MQ=255)
             gTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAAt                                                                                                                                        <  2:2109448/150‑1 (MQ=255)
                  gTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCGGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTg                                                                                                                                   >  1:640157/1‑150 (MQ=255)
                     gAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAg                                                                                                                                >  1:1196954/1‑150 (MQ=255)
                     gAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAg                                                                                                                                <  1:1793426/150‑1 (MQ=255)
                       aGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCAGTGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGt                                                                                                                              <  2:1288973/150‑1 (MQ=255)
                        gccgccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTc                                                                                                                             <  1:801747/150‑1 (MQ=255)
                           gccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGt                                                                                                                                    >  1:898418/1‑140 (MQ=255)
                           gccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGt                                                                                                                                    <  2:898418/140‑1 (MQ=255)
                            ccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGata                                                                                                                         <  1:671449/150‑1 (MQ=255)
                            ccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGata                                                                                                                         <  1:2191645/150‑1 (MQ=255)
                            ccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGAGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGata                                                                                                                         >  2:1749233/1‑150 (MQ=255)
                                 ccATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCATGTCTATATTCAt                                                                                                                    <  2:1296763/150‑1 (MQ=255)
                                  cATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATg                                                                                                                   <  1:2170469/150‑1 (MQ=255)
                                   aTTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATgg                                                                                                                  <  2:1549779/150‑1 (MQ=255)
                                    ttGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATggg                                                                                                                 >  2:182804/1‑150 (MQ=255)
                                        aGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCAACAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTg                                                                                                                                   <  2:1430046/128‑1 (MQ=255)
                                        aGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCAACAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTg                                                                                                                                   >  1:1430046/1‑128 (MQ=255)
                                             cTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAg                                                                                                              <  2:1837593/144‑1 (MQ=255)
                                             cTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAg                                                                                                              >  1:1837593/1‑144 (MQ=255)
                                             cTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAtt                                                                                                        <  2:1566693/150‑1 (MQ=255)
                                              ttGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTc                                                                                                           <  1:1595294/146‑1 (MQ=255)
                                              ttGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTc                                                                                                           >  2:1595294/1‑146 (MQ=255)
                                              ttGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAtgt                                                                                                       >  1:722956/1‑148 (MQ=255)
                                                gAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTa                                                                                                     <  1:2112605/150‑1 (MQ=255)
                                                      ggAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACt                                                                                               <  1:240223/150‑1 (MQ=255)
                                                         ggTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGgc                                                                                            <  2:2140494/150‑1 (MQ=255)
                                                           tCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGgcgc                                                                                          >  2:415880/1‑150 (MQ=255)
                                                              tGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGcgcc                                                                                       <  2:227346/150‑1 (MQ=255)
                                                                  cacggCGCCGCGCCGGCCAAGGACGAGACCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCaa                                                                                   <  2:1127980/147‑1 (MQ=255)
                                                                  cAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCaa                                                                                   >  1:1595913/1‑150 (MQ=255)
                                                                    gggCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAgg                                                                                 <  2:2294431/150‑1 (MQ=255)
                                                                    gggCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAgg                                                                                 <  2:125819/150‑1 (MQ=255)
                                                                      gcgcCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGgc                                                                               >  1:1087370/1‑150 (MQ=255)
                                                                       cgcCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGgcg                                                                              >  2:1746770/1‑150 (MQ=255)
                                                                           gAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGcc                                                                          >  1:1962518/1‑150 (MQ=255)
                                                                                 gCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCa                                                                                                               <  1:1174230/107‑1 (MQ=255)
                                                                                 gCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCa                                                                                                               <  1:1174300/107‑1 (MQ=255)
                                                                                 gCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCa                                                                                                               >  2:1174300/1‑107 (MQ=255)
                                                                                 gCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCa                                                                                                               >  2:1174230/1‑107 (MQ=255)
                                                                                  cccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAgg                                                                   >  1:1681554/1‑150 (MQ=255)
                                                                                   ccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGgcg                                                                                           <  2:1318229/125‑1 (MQ=255)
                                                                                   ccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGgcg                                                                                           >  1:1318229/1‑125 (MQ=255)
                                                                                    catgACCAGCCCCAGCAGCAAGTCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACGGAATCCGTGCAGGTCGATATTCCTGGGCAGGGCATTTTAGGGCCTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGc                                                                 >  2:629152/4‑150 (MQ=255)
                                                                                    cAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGgcgcgc                                                                           <  1:1305951/140‑1 (MQ=255)
                                                                                    cAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGgcgcgc                                                                           >  2:1305951/1‑140 (MQ=255)
                                                                                    cAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGc                                                                 >  2:273911/1‑150 (MQ=255)
                                                                                     aGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCg                                                                >  2:602143/1‑150 (MQ=255)
                                                                                          cAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCa                                                           >  2:1563405/1‑150 (MQ=255)
                                                                                          cAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCa                                                           >  1:2215305/1‑150 (MQ=255)
                                                                                            gCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCAtt                                                         <  1:415880/150‑1 (MQ=255)
                                                                                            gCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGAAAGTGCGCGCCGGTCAGGGCGGATCAtt                                                         <  1:933075/150‑1 (MQ=255)
                                                                                             ccccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAg                                                                                                    >  1:429104/1‑106 (MQ=255)
                                                                                             ccccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAg                                                                                                    <  2:429104/106‑1 (MQ=255)
                                                                                             ccccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTc                                                        <  1:1613203/150‑1 (MQ=255)
                                                                                              cccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATGAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCg                                                       >  2:1393196/1‑150 (MQ=255)
                                                                                               ccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGGCGG‑ACCGAATCCGTGCAGTTTGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGa                                                      <  2:1173228/150‑1 (MQ=255)
                                                                                                cagcagCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGAt                                                     >  1:1020553/1‑150 (MQ=255)
                                                                                                    agcaAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTg                                                 <  1:273911/150‑1 (MQ=255)
                                                                                                     gcaAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGa                                                >  2:1376177/1‑150 (MQ=255)
                                                                                                        aGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTc                                                        <  2:1431287/139‑1 (MQ=255)
                                                                                                        aGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTc                                                        >  1:1431287/1‑139 (MQ=255)
                                                                                                        aGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACg                                             <  1:182804/150‑1 (MQ=255)
                                                                                                          gcgcAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACggg                                           <  1:2174595/150‑1 (MQ=255)
                                                                                                           cgcAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGa                                          <  1:1746770/150‑1 (MQ=255)
                                                                                                            gcAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGaa                                         >  1:2194201/1‑150 (MQ=255)
                                                                                                               aaGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGGGCAGTTCGTTATTCGGGGGCAGGTCA‑TTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTAc                                     <  2:2104854/150‑1 (MQ=255)
                                                                                                                  cgcCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGTGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACtt                                   <  2:1196954/150‑1 (MQ=255)
                                                                                                                  cgcCAGTTTTACGGTCATCAGGGTGCCGGCCAGCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACtt                                   <  1:1486562/150‑1 (MQ=255)
                                                                                                                   gcCAGTTTTACGGTCATCAGGGTGCCGTCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCCAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAgatcggg                                  >  1:2357058/1‑144 (MQ=255)
                                                                                                                    ccAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGt                                 >  1:460320/1‑150 (MQ=255)
                                                                                                                     cAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTc                                <  1:1393196/150‑1 (MQ=255)
                                                                                                                       gTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGt                              <  1:938980/150‑1 (MQ=255)
                                                                                                                         tttACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTg                            <  2:2215305/150‑1 (MQ=255)
                                                                                                                           tACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAt                          <  2:1595913/150‑1 (MQ=255)
                                                                                                                                tcggcaGGGTGCCGGACATCCTGGCCGG‑ACACCAACAGTGAAGCTCGATAATACAGGGCAGGCCATTTTAGCGACTTGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCt                     <  2:1706286/146‑1 (MQ=255)
                                                                                                                                tcatcaGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCt                     <  2:417473/150‑1 (MQ=255)
                                                                                                                                  atcaGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGa                           >  2:1309248/1‑142 (MQ=255)
                                                                                                                                  atcaGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGa                           <  1:1309248/142‑1 (MQ=255)
                                                                                                                                  atcaGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAtt                         <  2:2111544/144‑1 (MQ=255)
                                                                                                                                  atcaGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAtt                         >  1:2111544/1‑144 (MQ=255)
                                                                                                                                        gTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATCCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTgtacgta             <  1:2359464/150‑1 (MQ=255)
                                                                                                                                         gtccGGCCATGATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGGACTTGTCGTTGATTTTCTTgtacgtac            >  1:1438363/3‑150 (MQ=255)
                                                                                                                                          tccGGCCATGATGGCCGGCACCGAATCCGTGCAGGTCGATATTCATGGGCAGGGCATTTTAGGGACTGGCGCGCCGCAAGGgc                                                                               >  1:1138347/2‑83 (MQ=38)
                                                                                                                                          tccGGCCATGATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTAcc           >  1:1087329/2‑150 (MQ=255)
                                                                                                                                          tccGGCCAGGATGGCCGGCACCGAATCCGTGCAGGTCGATATTCATGGGCAGGGCATTTTAGGGACTGGCGCGCCGCAAGGgc                                                                               <  2:1138347/82‑1 (MQ=255)
                                                                                                                                          gCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTAcc           >  2:501805/1‑150 (MQ=255)
                                                                                                                                            cGGCCATCATGGCCGG‑ACCTAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAtt                         >  1:2075317/1‑134 (MQ=255)
                                                                                                                                            cGGCCATCATGGCCGG‑ACCTAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAtt                         <  2:2075317/134‑1 (MQ=255)
                                                                                                                                                   catGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGCttt  <  1:750336/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCCGCGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGG‑ACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGCTTT  >  NC_002947/339081‑339377

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