Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,179,383 G→A K423K (AAG→AAA tolC → outer membrane channel TolC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,179,3830GA98.5% 227.9 / NA 68K423K (AAG→AAAtolCouter membrane channel TolC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/42);  minor base C (1/0);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTC  >  NC_000913/3179145‑3179623
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
gCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGt                                                                                                                                                                                                                                      <  1:121659/251‑1 (MQ=255)
  gtcgtGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTAc                                                                                                                                                                                                                                    <  1:408439/251‑1 (MQ=255)
       cccaaccgttgttccctcctTCAACAACATAAATGCATCTAGCAGTAGCATTACCGCCTCCAAACAAGCCGTAGTTTCCTCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGGTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGa                                                                                                                                                                                                                               <  2:16023/237‑1 (MQ=255)
        cAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGaa                                                                                                                                                                                                                              <  1:158207/251‑1 (MQ=255)
           aCCGTGCGTTCCTCCTTCAACAAAATTAAGGCACCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTTTACAACGCCAAGCAAGAGCTGGCTAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTTGGTACGTTTAACGa                                                                                                                                                                                                                           <  2:393697/251‑1 (MQ=255)
              gTGCGTTCCTCCTTAAACAACATTAATGCACCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCa                                                                                                                                                                                                                        <  2:169885/251‑1 (MQ=255)
                gCGTTCCCCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAg                                                                                                                                                                                                                       <  1:212124/250‑1 (MQ=255)
                        cctTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATctgctg                                                                                                                                                                                                              <  1:391549/251‑1 (MQ=255)
                        cctTCAACAACATTAATGCATCTATCAGTAGCATTAAAGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATctgctg                                                                                                                                                                                                              >  2:235693/1‑251 (MQ=255)
                             aaaaacaTTAATGAATCTATCAGTAGCATTAACGCCTCCAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTAGATGTGTTGGATGCGACCACCACTTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTGATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACt                                                                                                                                                                                                         <  2:3669/248‑1 (MQ=255)
                                 acaTTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAAc                                                                                                                                                                                                     >  1:167280/1‑251 (MQ=255)
                                  caTTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATGAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACa                                                                                                                                                                                                    <  2:292483/251‑1 (MQ=255)
                                   aTTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACa                                                                                                                                                                                                    >  2:133712/1‑250 (MQ=255)
                                     tAATGCATCCATCAGTAGCATAAACGCCTACAAACAAGCCATAGATTAAGCTACAAGCTCATTAGACCCGATGGAGAAGAGCTACTCGGACGGTACGCATACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGCGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCAGAATACTAAATCAGCTCTGGGTACGTTGAAAGAGCCGGATCTGCTGGAACTGAACAATg                                                                                                                                                                                                 >  2:41126/1‑251 (MQ=255)
                                         gCATCTAACATTAGCATTAACGCCGACAAACAAGCCGTAGTTGTCTCGCAAAGCTCATTAGACGCGATGGAAGCGGGCTACGCGGGCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACTCCAAGCAAGAGGTGGCGAATGCGCGTTATAACTACCTGATGAATCAGCTGAATATGAAATCAGCGCTTGGAACGTTGAACTAGCAGGATCTCATGGCAGGGAACAATGCGCt                                                                                                                                                                                             <  2:341810/251‑1 (MQ=255)
                                            tCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCCTTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAACTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTg                                                                                                                                                                                            >  2:326684/1‑249 (MQ=255)
                                               aTCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCa                                                                                                                                                                                        >  2:28919/1‑250 (MQ=255)
                                                acAGTATCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCaaa                                                                                                                                                                                      <  1:282088/250‑1 (MQ=255)
                                                  aGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTTCGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAAcc                                                                                                                                                                                    <  2:483842/251‑1 (MQ=255)
                                                     aGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGt                                                                                                                                                                                 >  2:34640/1‑251 (MQ=255)
                                                       cattcaaaccttcaaacaaGCCGTAGTTTAAACTCCACGCTCATTAGTCGCGATCAAATCGGGCTCCTCGGTCGGTACGCGTACCATTGTTCATGTGCTGGATGCGCCCACCACGTTGTACTACGCCAAGCAAGAGCTGGCTACTGCGCGTTATACCTACCTGATTAATCAGCTGAATATTAACTCCGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGCGCAAACCGGttt                                                                                                                                                                               >  2:80409/13‑251 (MQ=255)
                                                        aTTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTc                                                                                                                                                                              >  1:207496/1‑251 (MQ=255)
                                                                cTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATc                                                                                                                                                                      >  1:104813/1‑251 (MQ=255)
                                                                  acaaacaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATc                                                                                                                                                                      <  2:167280/249‑1 (MQ=255)
                                                                     aacaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGGACGCGTACCAATGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTCGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTTGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGaa                                                                                                                                                                 <  2:416603/251‑1 (MQ=255)
                                                                      acaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTGTCCACTAATCCGGaaa                                                                                                                                                                <  1:431955/251‑1 (MQ=255)
                                                                        aaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCAACACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGATCAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAAc                                                                                                                                                              <  1:132146/251‑1 (MQ=255)
                                                                         aGCCGTCGTTTCCGCTCAAAGCTCATTCGACGCGATGGCAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACg                                                                                                                                                             >  2:246919/1‑251 (MQ=255)
                                                                            cGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTg                                                                                                                                                          <  1:409216/251‑1 (MQ=255)
                                                                               aGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCa                                                                                                                                                        <  1:560799/250‑1 (MQ=255)
                                                                                   tCCGCTCAAAGCTCATTAGACGCGATGGAATCTGGCTACTCGGTCGGTTCGCGTACCATTTGTGATGTGTGGGATGCTACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCGGATTAATCAGCTGAATATGAAATCAGCGCTGGTTACGTTGAACGGGCAGGATCTGTGGGTACTGTACAATGGTCTGAGGAAACGTGATTCTACTAATCCGGAAAACGTTGTACCGc                                                                                                                                                    <  2:571027/250‑1 (MQ=255)
                                                                                        tCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGc                                                                                                                                              >  2:311810/1‑251 (MQ=255)
                                                                                                    aGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGCATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGAACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGAACCGCAAACGCCGGAACAGAATg                                                                                                                                  >  2:537083/1‑251 (MQ=255)
                                                                                                         cgATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTAtt                                                                                                                             >  2:454083/1‑251 (MQ=255)
                                                                                                          gATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTg                                                                                                                            <  1:16373/251‑1 (MQ=255)
                                                                                                                   gggCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGtt                                                                                                                   <  2:207496/251‑1 (MQ=255)
                                                                                                                        aCTCGGTCCGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATg                                                                                                                <  2:13352/249‑1 (MQ=255)
                                                                                                                               cggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGAt                                                                                                        <  2:82168/250‑1 (MQ=255)
                                                                                                                               cggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTGTCCACTAATCCGGAAAACGGTGCACCGCAAACGCCGGAACAGAATGCTATTGTTGGTGGTTATGCTCCTGATa                                                                                                       <  2:286509/251‑1 (MQ=255)
                                                                                                                                ggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGAt                                                                                                        >  2:414281/1‑249 (MQ=255)
                                                                                                                                     gcgTACCATTGTTGATGTGTTGGATGCGTCCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATTCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCgg                                                                                                 <  1:396983/251‑1 (MQ=255)
                                                                                                                                          ccATTGTTGATGTGTTTGATGCGAGCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCCCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCGGATGGTTATGCGCCTGATATCCCGGCACCa                                                                                            <  1:137529/251‑1 (MQ=255)
                                                                                                                                          acATTGTTGATGTGTTGGATGCGACCACCACGGTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGATCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTTCTGATGGTTATGCGCCTGATAGCCCGGCACCa                                                                                            <  2:510535/250‑1 (MQ=255)
                                                                                                                                           cATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGCTGGTTATGAGCCTGATAGCCCGGCACCa                                                                                            >  1:415476/1‑250 (MQ=255)
                                                                                                                                              tgttgATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCTGAAAACGGTGCACCGCAAACGCCGGAAGAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTc                                                                                         <  1:80409/250‑1 (MQ=255)
                                                                                                                                                        ttGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAAc                                                                              >  1:16666/1‑251 (MQ=255)
                                                                                                                                                          ggATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCTGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACa                                                                             <  2:508272/250‑1 (MQ=255)
                                                                                                                                                          ggATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACa                                                                             <  1:28919/250‑1 (MQ=255)
                                                                                                                                                             tGCGACAACCACGTTGTACAACGCCAAGCAAGAGCTGGCTAATGCGCGTTATAACTACCTGATTACTCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAACACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCcg                                                                         >  2:310119/1‑251 (MQ=255)
                                                                                                                                                                                    ccAAGAAAGAGCTGGCGACTGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGCAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGt                                                  >  2:555226/1‑251 (MQ=255)
                                                                                                                                                                                        gcaagAGCTGGCGAATGCGCGTTATAACTACCTGAGTAATCAGCTGAATATTAAATCAGCTCGGGGTACGTTGAACGAGCAGGATCTGCGGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATATCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATa                                              <  2:572465/251‑1 (MQ=255)
                                                                                                                                                                                              gCTGGAGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAACCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGCATGCTATTGCTGATGGTTATGCGCCTGATAGACAGGCACCAGTCGTTCAGCAAACCTCCGCACGCACTACCCCCAGTAACGGTCATAACCCtt                                        >  2:576359/1‑251 (MQ=255)
                                                                                                                                                                                                         cgcgTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACtga                             <  1:311810/251‑1 (MQ=255)
                                                                                                                                                                                                          gcgTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCTCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTCACtgat                            >  2:134767/1‑251 (MQ=255)
                                                                                                                                                                                                               ataACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATgacg                        <  1:321109/250‑1 (MQ=255)
                                                                                                                                                                                                                     aCCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTTCTGGCACTGAACAATGCGCTGAGCAAACCTGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCGATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAATCATCCGCACGCACTTCCACCAGTAACGGTCATAACCCTTCCCTTATCTGATGacgacgacg                  <  2:394158/250‑1 (MQ=255)
                                                                                                                                                                                                                     aCCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgg                 >  2:87620/1‑251 (MQ=255)
                                                                                                                                                                                                                     aCCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgg                 >  1:454646/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACggg                >  1:23505/1‑251 (MQ=255)
                                                                                                                                                                                                                         gATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgggg               <  1:133712/249‑1 (MQ=255)
                                                                                                                                                                                                                            tAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGACAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGATATGCGCCTGATAGCCCGGCACCCGTCGTTCAGCAAACATCcgca                                                                       >  1:566493/1‑190 (MQ=255)
                                                                                                                                                                                                                            tAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTc           <  1:34640/250‑1 (MQ=255)
                                                                                                                                                                                                                            tAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTc           <  1:585836/250‑1 (MQ=255)
                                                                                                                                                                                                                            tAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTGCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTTGTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCATACATCcgca                                                                       <  2:566493/190‑1 (MQ=255)
                                                                                                                                                                                                                                 agcttaagatTAAATCAGCTCTGGGTACGTTTAACGAGCAGGATCTGCTGTCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAAGCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGcccc     <  1:361245/243‑1 (MQ=255)
                                                                                                                                                                                                                                     taaTATTAAATCAGCTCTGGGGACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTGTCCACTAATCCTGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCTCCATTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTc  <  2:504582/249‑1 (MQ=255)
                                                                                                                                                                                                                                      aaTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg                                                                                                  >  1:202144/1‑153 (MQ=255)
                                                                                                                                                                                                                                      aaTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg                                                                                                  <  2:202144/153‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
GCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTC  >  NC_000913/3179145‑3179623

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