Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,178,483 G→C 6.1% V13L (GTT→CTT) ‡ rplK → 50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,178,4830GC 6.1% 450.4 / 14.7 132V13L (GTT→CTT) ‡rplK50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base G (60/64);  new base C (3/5);  total (63/69)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 132 total aligned reads displayed.

GATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGTTT  >  NC_000913/4178335‑4178630
                                                                                                                                                    |                                                                                                                                                   
gaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGt                                                                                                                                                    <  1:1190172/150‑1 (MQ=255)
gaTTGGAATACAATTTCGCCCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGt                                                                                                                                                    >  2:1384145/1‑150 (MQ=255)
    ggAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCa                                                                                                                                                <  2:768845/150‑1 (MQ=255)
         aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTgg                                                                                                                                           <  1:851946/150‑1 (MQ=255)
         aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAACCTGCAGCttg                                                                                                                                           <  1:529264/150‑3 (MQ=255)
          cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCCAGCTGCAACCTGCAGCttga                                                                                                                                          <  2:540364/150‑4 (MQ=255)
          cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGGAggtggt                                                                                                                                          >  1:668198/1‑150 (MQ=255)
          cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGt                                                                                                                                          >  2:1397581/1‑150 (MQ=255)
          cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGt                                                                                                                                          >  2:1285759/1‑150 (MQ=255)
          cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAACCTGCAGCttga                                                                                                                                          <  1:1089836/150‑4 (MQ=255)
            aTTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTAt                                                                                                                                        >  2:315216/1‑150 (MQ=255)
             tttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATg                                                                                                                                       <  1:922630/150‑1 (MQ=255)
               tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGc                                                                                                                                     <  1:265945/150‑1 (MQ=255)
               tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGc                                                                                                                                     <  2:1226264/150‑1 (MQ=255)
                  cgcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAA                                                                                                                                  <  2:627611‑M1/150‑1 (MQ=255)
                  cgcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAA                                                                                                                                  >  1:325699‑M1/1‑150 (MQ=255)
                   gcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTAGAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAAC                                                                                                                                 <  2:22143‑M1/150‑1 (MQ=255)
                    ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACC                                                                                                                                <  2:1110001‑M1/150‑1 (MQ=255)
                      ttttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGGTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCg                                                                                                                              <  1:1418677/150‑1 (MQ=255)
                         tGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGt                                                                                                                           >  1:1120711/1‑150 (MQ=255)
                          gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTc                                                                                                                          >  1:820464/1‑150 (MQ=255)
                          gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTc                                                                                                                          <  1:762516/150‑1 (MQ=255)
                          gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTc                                                                                                                          >  1:285687/1‑150 (MQ=255)
                          gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTc                                                                                                                          <  1:994762/150‑1 (MQ=255)
                              ttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCATAGGCGTTATTACCCAACTTGAGTAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGcc                                                                                                                      <  2:1492923/150‑1 (MQ=255)
                              ttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGcc                                                                                                                      <  1:942450/150‑1 (MQ=255)
                              ttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGcc                                                                                                                      <  1:1285759/150‑1 (MQ=255)
                                  gggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTa                                                                                                                  >  2:744429/1‑150 (MQ=255)
                                  gggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTa                                                                                                                  <  2:1161055/150‑1 (MQ=255)
                                   ggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAg                                                                                                                 <  2:973982/150‑1 (MQ=255)
                                   ggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAg                                                                                                                 >  1:764976/1‑150 (MQ=255)
                                     ccTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGt                                                                                                               >  2:1139892/1‑150 (MQ=255)
                                       ttGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTcc                                                                                                             <  1:389951/150‑1 (MQ=255)
                                        tGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCg                                                                                                            <  2:1169743/150‑1 (MQ=255)
                                          tgcccgtaaCGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGc                                                                                                          >  1:324555/8‑150 (MQ=255)
                                            cGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGctc                                                                                                        >  1:1338385/1‑150 (MQ=255)
                                                aaaCGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTggg                                                                                                    >  1:1015416‑M2/1‑147 (MQ=255)
                                                     aTTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGGTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTcagc                                                                                               <  2:1446039/150‑1 (MQ=255)
                                                              atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      >  2:1334693/1‑150 (MQ=255)
                                                              atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      <  1:1410603/150‑1 (MQ=255)
                                                              atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      >  2:220739/1‑150 (MQ=255)
                                                                cACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAAc                                                                                    >  1:129369/1‑150 (MQ=255)
                                                                 aCGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCCGGGCGTAAAca                                                                                   <  2:1022841/150‑1 (MQ=255)
                                                                    gggAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAAcatca                                                                                >  1:1261945/1‑150 (MQ=255)
                                                                     ggAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAAcatcat                                                                               <  2:324555/150‑1 (MQ=255)
                                                                      gAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATg                                                                              <  1:236939/150‑1 (MQ=255)
                                                                      gAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATg                                                                              >  2:278757/1‑150 (MQ=255)
                                                                         ccTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGaa                                                                           >  2:841232/1‑150 (MQ=255)
                                                                           tctcAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTTAGCAGTGCGTGAACATCATGGAAtt                                                                         >  1:879045/1‑150 (MQ=255)
                                                                            ctcAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTc                                                                        >  2:788280/1‑150 (MQ=255)
                                                                              cAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCagagccggacctactggcggactcgggttagccat                                                                                            >  1:120474‑M1/1‑93 (MQ=255)
                                                                              cAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCagagccggacctactggcggactcgggttagccat                                                                                            <  2:120474‑M1/128‑36 (MQ=255)
                                                                               agagGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGc                                                                     <  1:334082/150‑1 (MQ=255)
                                                                               agagGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGc                                                                     >  2:120029/1‑150 (MQ=255)
                                                                                gagGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCa                                                                    <  1:568710/150‑1 (MQ=255)
                                                                                cccgCGTCATTACCCCACTTGAGGAATTTATAATGCCTAAGAAAGTACAAGCCTATGTCAACCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCa                                                                    <  2:668198/147‑1 (MQ=255)
                                                                                 agGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCaa                                                                   <  1:1360463/150‑1 (MQ=255)
                                                                                    cGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGc                                                                <  1:1092298/150‑1 (MQ=255)
                                                                                     gTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCg                                                               >  1:1529487/1‑150 (MQ=255)
                                                                                     gTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCg                                                               >  2:921191/1‑150 (MQ=255)
                                                                                      ttattaCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGt                                                              >  1:1314018/1‑150 (MQ=255)
                                                                                        attaCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATTATGGAATTCTGCAAAGCGTTc                                                            >  2:859203/1‑150 (MQ=255)
                                                                                        attaCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTc                                                            <  1:1334693/150‑1 (MQ=255)
                                                                                          taCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCaa                                                          <  1:152615/150‑1 (MQ=255)
                                                                                           aCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAAc                                                         <  2:525042/150‑1 (MQ=255)
                                                                                           aCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAAc                                                         >  2:823353/1‑150 (MQ=255)
                                                                                           aCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGAAGCTGGTATGGCTAACCCGACTCCGCCAGTAGGTCCGGCTCTGTGTCAGCAGGGCGTAAACATCAAGGAATTCTGCAAAGCGTTCAAc                                                         >  2:1021733/1‑150 (MQ=255)
                                                                                            cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACg                                                        <  2:1534719/150‑1 (MQ=255)
                                                                                            cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACg                                                        <  1:751920/150‑1 (MQ=255)
                                                                                            cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCAGAGTCCGCCAGAACGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAAcc                                                        >  2:1446177/1‑149 (MQ=255)
                                                                                              cAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTca                                                                                                 <  2:751358/107‑1 (MQ=255)
                                                                                              cAACTTGAGGAATTTATAATGGATAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCa                                                      >  2:1375714/1‑150 (MQ=255)
                                                                                                 cTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCaaaa                                                   >  1:1375210/1‑150 (MQ=255)
                                                                                                    gAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACt                                                >  1:73027/1‑150 (MQ=255)
                                                                                                        aaTTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGAtt                                            <  1:1153264/150‑1 (MQ=255)
                                                                                                           ttATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCa                                         <  2:1173815/150‑1 (MQ=255)
                                                                                                            tataATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCAt                                        >  1:153207/1‑150 (MQ=255)
                                                                                                              taATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCg                                      <  1:1397581/150‑1 (MQ=255)
                                                                                                                     tAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAgg                               >  2:100775/1‑150 (MQ=255)
                                                                                                                      aaGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGt                              >  2:575245/1‑150 (MQ=255)
                                                                                                                        gAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCt                            >  1:825424/1‑150 (MQ=255)
                                                                                                                          aaGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGc                          <  1:575245/150‑1 (MQ=255)
                                                                                                                           aGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGcc                         <  1:744429/150‑1 (MQ=255)
                                                                                                                            gTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCg                        >  2:219413/1‑150 (MQ=255)
                                                                                                                               cAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCTAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGAtt                     <  1:315216/150‑1 (MQ=255)
                                                                                                                               cAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGAtt                     >  1:378275/1‑150 (MQ=255)
                                                                                                                                 aGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTcc                   <  2:560498/150‑1 (MQ=255)
                                                                                                                                 aGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTcc                   <  2:1108105/150‑1 (MQ=255)
                                                                                                                                   ccTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCgg                 <  2:1500724/150‑1 (MQ=255)
                                                                                                                                    cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                >  1:1177447/1‑150 (MQ=255)
                                                                                                                                    cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                <  2:1261945/150‑1 (MQ=255)
                                                                                                                                    cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                <  1:1106772/150‑1 (MQ=255)
                                                                                                                                      cataccaGCTGCAACCTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgtag              >  1:1213010/7‑150 (MQ=255)
                                                                                                                                       ataccaGCTGCAACCTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGAt                      <  2:1247734/136‑1 (MQ=255)
                                                                                                                                        gTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgtagta            >  2:1524647/1‑150 (MQ=255)
                                                                                                                                         tCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTaa           <  1:1375714/150‑1 (MQ=255)
                                                                                                                                          cAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAAt          >  2:913379/1‑150 (MQ=255)
                                                                                                                                             gCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCAc       >  2:866803/1‑150 (MQ=255)
                                                                                                                                               tGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCg     >  2:732903/1‑150 (MQ=255)
                                                                                                                                                  aGGTTGCAGCTTGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGttt  <  1:634955/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGTTT  >  NC_000913/4178335‑4178630

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