Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,178,482 G→A 6.0% Q12Q (CAG→CAA 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,4820GA 6.0% 425.4 / 9.2 132Q12Q (CAG→CAArplK50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base G (60/64);  new base A (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 = 9.91e-01

Only 100 of 133 total aligned reads displayed.

GAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGTTT  >  NC_000913/4178333‑4178630
                                                                                                                                                     |                                                                                                                                                    
gagaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATAACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAg                                                                                                                                                      >  2:634955/1‑150 (MQ=255)
  gaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGt                                                                                                                                                    <  1:1190172/150‑1 (MQ=255)
  gaTTGGAATACAATTTCGCCCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGt                                                                                                                                                    >  2:1384145/1‑150 (MQ=255)
    ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTg                                                                                                                                                  >  2:1394671/1‑150 (MQ=255)
     tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGc                                                                                                                                                 <  2:272988/150‑1 (MQ=255)
      ggAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCa                                                                                                                                                <  2:768845/150‑1 (MQ=255)
            cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCCAGCTGCAACCTGCAGCttga                                                                                                                                          <  2:540364/150‑4 (MQ=255)
            cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGGAggtggt                                                                                                                                          >  1:668198/1‑150 (MQ=255)
            cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGt                                                                                                                                          >  1:3170/1‑150 (MQ=255)
            cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGt                                                                                                                                          >  2:1285759/1‑150 (MQ=255)
            cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGt                                                                                                                                          >  2:1397581/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                                                                                                                                     <  2:1226264/150‑1 (MQ=255)
                 tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGc                                                                                                                                     <  1:265945/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)
                        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:994762/150‑1 (MQ=255)
                            gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTc                                                                                                                          >  1:285687/1‑150 (MQ=255)
                                ttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCATAGGCGTTATTACCCAACTTGAGTAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGcc                                                                                                                      <  2:1492923/150‑1 (MQ=255)
                                ttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGcc                                                                                                                      <  1:1285759/150‑1 (MQ=255)
                                    gggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTa                                                                                                                  <  2:1161055/150‑1 (MQ=255)
                                    gggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTa                                                                                                                  >  2:744429/1‑150 (MQ=255)
                                     ggCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAg                                                                                                                 <  2:973982/150‑1 (MQ=255)
                                       ccTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGt                                                                                                               >  2:1139892/1‑150 (MQ=255)
                                         ttGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTcc                                                                                                             <  1:389951/150‑1 (MQ=255)
                                            tgcccgtaaCGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGc                                                                                                          >  1:324555/8‑150 (MQ=255)
                                              cGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGctc                                                                                                        >  1:1338385/1‑150 (MQ=255)
                                                  aaaCGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTggg                                                                                                    >  2:561534‑M2/1‑147 (MQ=255)
                                                  aaaCGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTggg                                                                                                    >  1:1015416‑M2/1‑147 (MQ=255)
                                                          ttttATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTGACCCagagccggacctactggcggactcgggttagccat                                                                                            <  2:469310‑M1/150‑36 (MQ=255)
                                                                atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      >  2:220739/1‑150 (MQ=255)
                                                                atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      <  1:1410603/150‑1 (MQ=255)
                                                                atCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTaa                                                                                      >  2:1334693/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)
                                                                          gCCTCTCAGAGGCGTTATTACCCAACTCGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGa                                                                            >  2:1153264/1‑150 (MQ=255)
                                                                           ccTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGaa                                                                           >  2:841232/1‑150 (MQ=255)
                                                                            ctctCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCTGTTGGTTTGGTTATCCGGAGTCCGCCTGTAGGTCCGGGTGTTGGTCAGCTGGGCGTAAACTTCTTGGAAt                                                                          <  2:1066620/150‑1 (MQ=255)
                                                                             tctcAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTTAGCAGTGCGTGAACATCATGGAAtt                                                                         >  1:879045/1‑150 (MQ=255)
                                                                              ctcAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTc                                                                        >  2:788280/1‑150 (MQ=255)
                                                                                cAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAATCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTg                                                                      >  2:129578/1‑150 (MQ=255)
                                                                                cAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCagagccggacctactggcggactcgggttagccat                                                                                            >  1:120474‑M1/1‑93 (MQ=255)
                                                                                cAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTg                                                                      >  1:560498/1‑150 (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                                                               >  2:921191/1‑150 (MQ=255)
                                                                                       gTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCg                                                               >  1:1529487/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)
                                                                                          attaCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTc                                                            >  2:1247517/1‑150 (MQ=255)
                                                                                            taCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCaa                                                          <  1:152615/150‑1 (MQ=255)
                                                                                             aCCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAAc                                                         >  2:823353/1‑150 (MQ=255)
                                                                                              cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACg                                                        <  1:751920/150‑1 (MQ=255)
                                                                                              cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACg                                                        <  2:1534719/150‑1 (MQ=255)
                                                                                              cccAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCAGAGTCCGCCAGAACGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAAcc                                                        >  2:1446177/1‑149 (MQ=255)
                                                                                                cAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTca                                                                                                 >  1:751358/1‑107 (MQ=255)
                                                                                                cAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTca                                                                                                 <  2:751358/107‑1 (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)
                                                                                                                   tGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGaaa                                   <  1:942841/150‑1 (MQ=255)
                                                                                                                    ggCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGaaaa                                  >  2:1461983/1‑150 (MQ=255)
                                                                                                                       tAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAgg                               >  2:100775/1‑150 (MQ=255)
                                                                                                                          gAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCt                            >  1:825424/1‑150 (MQ=255)
                                                                                                                            aaGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGc                          <  1:575245/150‑1 (MQ=255)
                                                                                                                                 cAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGAtt                     >  1:378275/1‑150 (MQ=255)
                                                                                                                                   aGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTcc                   <  2:1108105/150‑1 (MQ=255)
                                                                                                                                   aGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTcc                   <  2:560498/150‑1 (MQ=255)
                                                                                                                                     ccTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCgg                 <  1:340033/150‑1 (MQ=255)
                                                                                                                                      cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                >  1:1177447/1‑150 (MQ=255)
                                                                                                                                      cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                <  1:1106772/150‑1 (MQ=255)
                                                                                                                                      cTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGgt                <  2:1179560/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)
                                                                                                                                           accaGCTGCAACCTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTaa           >  1:860460/4‑150 (MQ=255)
                                                                                                                                            cAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAAt          >  2:913379/1‑150 (MQ=255)
                                                                                                                                              aGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAACAGGTCTGCCGATTCCGGTAGTAATCa        >  1:664425/1‑150 (MQ=255)
                                                                                                                                               gCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCAc       >  2:866803/1‑150 (MQ=255)
                                                                                                                                                 tGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCg     >  2:732903/1‑150 (MQ=255)
                                                                                                                                                    aGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGttt  <  1:788280/150‑1 (MQ=255)
                                                                                                                                                    aGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGttt  <  2:1529487/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGAATTTATAATGGCTAAGAAAGTACAAGCCTATGTCAAGCTGCAGGTTGCAGCTGGTATGGCTAACCCGAGTCCGCCAGTAGGTCCGGCTCTGGGTCAGCAGGGCGTAAACATCATGGAATTCTGCAAAGCGTTCAACGCAAAAACTGATTCCATCGAAAAAGGTCTGCCGATTCCGGTAGTAATCACCGTTT  >  NC_000913/4178333‑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: