Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,653 A→G 5.4% L250L (CTA→CTG yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6530AG 5.4% 240.5 / 3.9 92L250L (CTA→CTGyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (40/45);  new base G (2/3);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

AGAGGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAT  >  NC_000913/2357504‑2357800
                                                                                                                                                     |                                                                                                                                                   
agagGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCAGCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTa                                                                                                                                                     >  1:434160/1‑150 (MQ=255)
    gCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGc                                                                                                                                                 <  2:357926/150‑1 (MQ=255)
     gctCACTAATGACAGCAATCTAGAAACACTGTTTAAGTATCTGGTGCTATGCGGGGATAGCTCGCGTTCCACCCGCTTACTTCCTGATGCTGCCGAACCGTCACCTCTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCa                                                                                                                                                <  2:761038/148‑1 (MQ=16)
        cACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGa                                                                                                                                             >  2:670371/1‑150 (MQ=255)
        cACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCAGCAATGCGGCGATACCTCCCGTTTCACCCGTTATATTGAGGACATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGa                                                                                                                                             >  2:444074/1‑150 (MQ=255)
          cTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTGCCCCCGTTTTATTGAGGAGATTCCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGaat                                                                                                                                           >  2:235340/1‑149 (MQ=255)
              tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCa                                                                                                                                       >  1:1228143/1‑150 (MQ=255)
              tGACAGCAAGCTACAAACACTGTTTAATTATCTGCGGCAATTCGGCGATACCTCCCGTTCGACCAGTTTTCTTGAGGAGATTGCCGAACGTTCACCACTAAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACa                                                                                                                                       <  2:879627/150‑1 (MQ=37)
                aCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATc                                                                                                                                     >  1:1070702/1‑150 (MQ=255)
                   gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                  >  1:442823/1‑150 (MQ=255)
                   gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTCCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                  >  1:904954/1‑150 (MQ=255)
                     cagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGGTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAtt                                                                                                                                <  2:519003/150‑1 (MQ=255)
                       gcTACAAACACTGTTTAATTAACTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCCCCACTACAAAAGGAGAGATTAAAGACTATAGCTGAACGGCTCCCGCAGGAAGTACATCAAATtgg                                                                                                                              >  2:308873/1‑150 (MQ=255)
                         tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                            >  1:69096/1‑150 (MQ=255)
                         tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTATATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                            <  1:38286/150‑1 (MQ=255)
                          acacacacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCCTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGCGATTAATGACTATTCCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                           <  2:185454/150‑1 (MQ=255)
                          aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                           <  2:972801/150‑1 (MQ=255)
                           cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCGCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTGCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                          <  1:637894/150‑1 (MQ=255)
                            aaaCACTGTTTAATTATCTGCTGCAATGCGTCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGTCTACGGCAGGAAGGGCATCAAATttgctggc                                                                                                                         <  1:198491/150‑1 (MQ=255)
                            aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGTCAGGAAGGGCATCAAATtggctggc                                                                                                                         <  1:735443/150‑1 (MQ=255)
                                   gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                  >  1:1038800/1‑150 (MQ=255)
                                   gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                  >  1:20326/1‑150 (MQ=255)
                                   gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                  >  2:833678/1‑150 (MQ=255)
                                    tttatttaTCAGCTGCAATGAGGCGATACCTCGCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGt                                                                                                                 <  2:20326/150‑1 (MQ=25)
                                      tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAt                                                                                                               >  2:1293934/1‑150 (MQ=255)
                                        aTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGc                                                                                                             <  2:345680/150‑1 (MQ=255)
                                         ttATCTGCGGCAATGCGGCGATACCTCCCGGTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCa                                                                                                            <  1:679509/150‑1 (MQ=255)
                                               gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGCCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACa                                                                                                      >  1:840311/1‑150 (MQ=255)
                                                 tgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAg                                                                                                    <  2:904954/150‑1 (MQ=255)
                                                  gcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGc                                                                                                   <  2:1188287/150‑1 (MQ=255)
                                                   cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGcc                                                                                                  <  1:271788/150‑1 (MQ=255)
                                                   cAATGCGGCGATACCTCCAGTTTCACCCGTATTATTGATGAGATTGCAGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAGGAACAAGcc                                                                                                  <  1:194122/150‑1 (MQ=37)
                                                    aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACTTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATTAACAAGCCa                                                                                                 <  2:785519/150‑1 (MQ=255)
                                                    aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCa                                                                                                 >  2:854431/1‑150 (MQ=255)
                                                    aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCa                                                                                                 <  1:443912/150‑1 (MQ=255)
                                                       gcggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTa                                                                                              >  1:699203/1‑150 (MQ=255)
                                                        cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaa                                                                                             >  1:1005782/1‑150 (MQ=255)
                                                          gcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                           >  2:791933/1‑150 (MQ=255)
                                                            taTACCTCCCGTTTCACCCGTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTg                                                                                        <  2:745474/149‑1 (MQ=255)
                                                              tACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGc                                                                                       <  2:351953/150‑1 (MQ=255)
                                                                ccTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCtt                                                                                     >  1:199443/1‑150 (MQ=255)
                                                                 cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                    >  1:990591/1‑150 (MQ=255)
                                                                 cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGACATTAAAATTGCttt                                                                                    >  1:294164/1‑150 (MQ=255)
                                                                 cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGATACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATAAATGAACAAGCCATTAAAATTGcgat                                                                                    >  1:1124443/1‑147 (MQ=255)
                                                                 cTCCCGTGTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                    <  2:679745/150‑1 (MQ=255)
                                                                 cTCCCGATTCACCCGTTTTATTGAGGAGATTGCCGAACGTACACCACTACAAAACGAGAGACTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGTCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                    <  1:381945/150‑1 (MQ=37)
                                                                   cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCGAATAAACTTGCGTTgc                                                                                  >  2:1183613/1‑150 (MQ=255)
                                                                   cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                  <  2:699203/150‑1 (MQ=255)
                                                                      gTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATTAATGAACAATACATTAAAATAGCTTTgagaa                                                                               >  2:376194/1‑146 (MQ=25)
                                                                      gTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATACATGAACAAGAGATTAAAATTGCTTTgataa                                                                               >  2:929076/1‑146 (MQ=255)
                                                                        ttCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATg                                                                             <  1:882578/150‑1 (MQ=255)
                                                                            cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTgg                                                                         >  1:454458/1‑150 (MQ=255)
                                                                            cccGTTTTATTGAGGAGATTGCCGAACGGTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCGACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTgg                                                                         <  2:710475/150‑1 (MQ=255)
                                                                             ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTTAAATTGCGTTGCGCATGCTGGa                                                                        >  1:626184/1‑150 (MQ=255)
                                                                             ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGa                                                                        <  1:667825/150‑1 (MQ=255)
                                                                               gctttCTTGAGGAGATAGCCGACCGTTCCCCACTACAAAAGCAGAGATTAATGACTATTGCTGACCGGCTACGGCAGGAACGGCAACAAATTTGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAACTTGCTTTGCGCATGCTGGAGc                                                                      <  1:1103907/148‑1 (MQ=25)
                                                                                  ttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAgg                                                                   <  2:703047/150‑1 (MQ=255)
                                                                                  ttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAgg                                                                   >  2:1086090/1‑150 (MQ=255)
                                                                                        cggagATTCCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTAGTGCAGGAAGTGCATCATATTGGCTTGCAGGAAGGTTTGGATGAACAAGCCATTATAATTGCTTTGCGCATGCTGGAGCAGGGCTTTg                                                             <  1:623857/149‑1 (MQ=25)
                                                                                        aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTg                                                             <  2:1038800/150‑1 (MQ=255)
                                                                                        aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATAATTGAACGATGCATTAAAAATGCATTGAAAATTCTTGAGCAGAGCATTg                                                             >  2:235348/1‑150 (MQ=18)
                                                                                           agaTTGCCGAACGTTCACCACTACAAAAGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                          <  2:436957/150‑1 (MQ=255)
                                                                                            gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                         <  2:454458/150‑1 (MQ=255)
                                                                                             aTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGt                                                        <  2:426239/150‑1 (MQ=255)
                                                                                                gCCGAACGTTCACCACTACAAAAGGAGAGAATAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTgag                                                     <  1:1159234/150‑1 (MQ=255)
                                                                                                       gTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAATACATTAAAAATGCTTTGAGCATGCTGGAGCAGGGCTTTGATCGTGACAATGTGc                                              >  1:501347/1‑150 (MQ=21)
                                                                                                         tCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCATTGAGCATGCAGGTGCAGGGCTTTGAACGTGACCTTGAGCTg                                            >  2:1305986/1‑150 (MQ=25)
                                                                                                           accacTACAAAAGGAGAGATTAATGACTATTGTTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                          <  1:914841/150‑1 (MQ=255)
                                                                                                            ccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCg                                         <  1:854431/150‑1 (MQ=255)
                                                                                                              acTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGATTGAACAATACATTAAAATTGCTATGAGAATGCTGTAGCACGGCATTGAACGTGACTTTGAGCtagagaa                                       >  2:979331/1‑144 (MQ=18)
                                                                                                              acTACAAAAGGAGAGATTAATGACTAGTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGCTATGCTTGAACAAGCCATTAAAATAGCTTTACGCATGCGTGAGCAAGGCATTGATCCTGACATTGTGcgcgcccc                                       >  2:1292282/1‑143 (MQ=21)
                                                                                                                 aCAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCGATTAAACTAGCGTTGCGAATTCAGTAGCAGAGCTTTGAACGTGACCTTGAGCTGGCGaacaa                                    >  2:719322/1‑146 (MQ=18)
                                                                                                                  cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTATGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGACATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGACCTTGTGCTCGCGAccccc                                   >  2:1107233/1‑146 (MQ=21)
                                                                                                                   aaaaGGAGAGATTAATGACTATTGCTGAACGTCTACGGCAGGAAGGGCATCAAATTGGCTTGCAGGTAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACTTGATATTGTGCTGGCGACAAccc                                  <  2:557498/150‑1 (MQ=25)
                                                                                                                   aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAccc                                  <  2:990591/150‑1 (MQ=255)
                                                                                                                    aaaGGAGAGATTAATGACTATTGCTGAACGGCTACTGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCa                                 <  2:958737/150‑1 (MQ=255)
                                                                                                                      tggAGAGATTAATGACTATTGCTGATCGTCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAATGTATGCATGAACAAGCCATTAAAATTGGTTTGCGCATGCTGGAGCTGGGCTTTGAACGTTAGATTGTGCTGGCGACAACCCaaa                               <  2:429803/149‑2 (MQ=37)
                                                                                                                      aGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCTTGAACAAGACAATAACATTGCTTTGAGCATGGGGGTGAAGGGCTTTGACCGTGACATGGTGCGGGCGGCagctctcg                               >  2:641696/1‑143 (MQ=17)
                                                                                                                        gagagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGAATGCTGGAGCAGGGCTTTGAACTTGACATTGTGCTAGCGACCACCCACCta                             >  1:289824/1‑149 (MQ=21)
                                                                                                                            gaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCTTGAACAAGCCATTAAAATTGCTTTTCTCATGCTGGAGCAGGGCTTTGAACGTGACAATGTGCTAGCGACAACCCACCTAAcaa                         >  1:249812/1‑148 (MQ=21)
                                                                                                                             aTTTATGACTATTGCTGATCGGCTACGGCAGGAAGGGCATCATATTGGCTGGCAGGAAGGTATGCATGAACAAGACATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGATCGTGTGATTGTGCTGGCGACAAGCCAACTCACTGa                        <  1:205321/150‑1 (MQ=37)
                                                                                                                             aTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGa                        <  1:1086090/150‑1 (MQ=255)
                                                                                                                                aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACGAGCCATTAAAATTGCTTTTAGCATGCTCGAGCAGGGCTTTGAACGTGACAATGTTCTGGAGACAACCCATCTAACTGAAGc                     >  1:556937/1‑150 (MQ=25)
                                                                                                                                 aTGACTATTTCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAGGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt                    <  1:1142674/150‑1 (MQ=255)
                                                                                                                                        ttccTGAACGGCTACGGCAGTAATGGCAGCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGTAGCAGGGCGTTGAACGTGAGATTGTGCTGGCGACAACCCACCTCACTGATGCGGATATTc             <  2:967855/147‑1 (MQ=37)
                                                                                                                                          gCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGAATGCTGGAACAGGGCTTTGAACGTGACAATGTGCTAGCGCCCACCCAGCTAACCGAAGCCGatctcgca           >  2:598965/1‑144 (MQ=17)
                                                                                                                                           cTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                         >  2:14775/1‑103 (MQ=255)
                                                                                                                                           cTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                         <  1:14775/103‑1 (MQ=255)
                                                                                                                                           cTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGACAATGTGCTCGCGGCAACCCAACTCAGAGAAGCAGATATTccac          >  1:1311901/1‑148 (MQ=21)
                                                                                                                                             gatcGGCTACGGCTGGATGTGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAAc        <  2:468979/147‑1 (MQ=255)
                                                                                                                                             gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAAc        <  1:686772/150‑1 (MQ=255)
                                                                                                                                                 ggCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGATCGTGACATTGTGCTGGCGCCCACCCACCTAACAGAAGCGGAAATACCGAACTgac    >  1:847136/1‑148 (MQ=21)
                                                                                                                                                   cTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAt  <  1:791933/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGAGGCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAT  >  NC_000913/2357504‑2357800

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