Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,798,659 (T)8→7 coding (654/1074 nt) waaU ← lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,6520T.98.2% 508.3 / 3.2 113coding (661/1074 nt)waaUlipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (50/60);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 125 total aligned reads displayed.

GATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGTAACGCTTATACCAATGA  >  NC_000913/3798421‑3798872
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            
gATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGTTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       <  1:800632/239‑2 (MQ=255)
    aGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgca                                                                                                                                                                                                                     >  1:1183996/1‑234 (MQ=255)
    aGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       <  2:2255105/235‑2 (MQ=255)
    aGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       <  2:623879/235‑2 (MQ=255)
    aGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       >  1:2255105/1‑234 (MQ=255)
    aGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       >  1:623879/1‑234 (MQ=255)
       cGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  2:1467257/239‑1 (MQ=255)
       cGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  1:1397610/239‑1 (MQ=255)
       cGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAAAGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  1:1017962/239‑1 (MQ=255)
           gCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGa                                                                                                                                                                                                           <  1:1273091/239‑1 (MQ=255)
               tatGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaa                                                                                                                                                                                                       >  1:1306660/1‑239 (MQ=255)
                atGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaat                                                                                                                                                                                                      <  2:61970/239‑1 (MQ=255)
                    aCAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTaataat                                                                                                                                                                                                   <  2:136307/238‑1 (MQ=255)
                      aaGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGAc                                                                                                                                                                                                <  1:1147384/239‑1 (MQ=255)
                       agagCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACt                                                                                                                                                                                               <  2:2290695/239‑1 (MQ=255)
                        gagCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTa                                                                                                                                                                                                                >  2:646338/1‑221 (MQ=255)
                        gagCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATg                                                                                                                                                                                                  >  2:1147384/1‑235 (MQ=255)
                         agCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTaataat                                                                                                                                                                                                   <  2:2144544/233‑1 (MQ=255)
                         agCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTaataat                                                                                                                                                                                                   >  1:2144544/1‑233 (MQ=255)
                             gTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTCACTTAATGGATTAATAATGACTAt                                                                                                                                                                                             >  1:659982/1‑235 (MQ=255)
                                      gAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTAt                                                                                                                                                                                             >  2:1149831/1‑226 (MQ=255)
                                      gAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTAt                                                                                                                                                                                             <  1:1149831/226‑1 (MQ=255)
                                            aCAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                                                                                                                                                           <  1:546840/238‑1 (MQ=255)
                                            aCAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                                                                                                                                                           >  2:546840/1‑238 (MQ=255)
                                                 aTCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTCACTTAATGGCTTAATAATGACTATTCTTGTATCTCCAATAAATTAttt                                                                                                                                                                     >  1:2268835/1‑239 (MQ=255)
                                                     cTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatt                                                                                                                                                                 <  2:1414882/239‑1 (MQ=255)
                                                             ggTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  2:2294700/185‑1 (MQ=255)
                                                                caacgcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAAc                                                                                                                                                      <  1:821323/239‑1 (MQ=255)
                                                                  acgcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACa                                                                                                                                                     >  1:2305357/1‑238 (MQ=255)
                                                                  acgcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATCTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATACATTCTTTTATTTTATCTTCAACAt                                                                                                                                                    >  2:706852/1‑239 (MQ=255)
                                                                   cgcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGa                                                                                                                                                                                                 >  2:2175555/1‑193 (MQ=255)
                                                                                    aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGa                                                                                                                                                                                                 <  2:1844061/176‑1 (MQ=255)
                                                                                    aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGa                                                                                                                                                                                                 >  1:1844061/1‑176 (MQ=255)
                                                                                     cTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgc                                                                                                                                                                                                                      <  2:2187087/155‑3 (MQ=255)
                                                                                     cTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgc                                                                                                                                                                                                                      >  1:2187087/1‑153 (MQ=255)
                                                                                       catcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       >  1:1547761/1‑151 (MQ=255)
                                                                                       catcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAttcttttg                                                                                                                                                                                                                       <  2:1547761/152‑2 (MQ=255)
                                                                                         tcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAg                                                                                                                             <  2:1306660/239‑1 (MQ=255)
                                                                                                          gATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTAcc                                                                                                            <  1:755319/239‑1 (MQ=255)
                                                                                                           aTTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGa                                                                                                                                                                                                           <  1:66625/143‑1 (MQ=255)
                                                                                                                          attgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               >  2:1723468/1‑124 (MQ=255)
                                                                                                                             gTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  2:328511/121‑1 (MQ=255)
                                                                                                                             gTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               >  1:328511/1‑121 (MQ=255)
                                                                                                                                        ttGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTAt                                                                              <  1:474828/239‑1 (MQ=255)
                                                                                                                                        ttGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTAt                                                                              <  1:527057/239‑1 (MQ=255)
                                                                                                                                          gCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTAtt                                                                             >  2:327919/1‑238 (MQ=255)
                                                                                                                                            ggTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGc                                                                           >  2:819940/1‑238 (MQ=255)
                                                                                                                                                aCCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCgtg                                                                      <  1:1082367/239‑1 (MQ=255)
                                                                                                                                                              aTCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCt                                                        >  2:1175850/1‑239 (MQ=255)
                                                                                                                                                               tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               >  2:437779/1‑87 (MQ=255)
                                                                                                                                                               tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGa                                                                                                                                                                                                           <  1:791427/91‑1 (MQ=255)
                                                                                                                                                               tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGa                                                                                                                                                                                                           >  2:791427/1‑91 (MQ=255)
                                                                                                                                                               tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAtc                                                                                                                                                                                    <  1:311482/114‑1 (MQ=255)
                                                                                                                                                               tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAtc                                                                                                                                                                                    >  2:311482/1‑114 (MQ=255)
                                                                                                                                                                  gATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtt                                                                                                                                                                      <  1:1903655/125‑1 (MQ=255)
                                                                                                                                                                    tAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               <  2:108327/82‑1 (MQ=255)
                                                                                                                                                                    tAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                                                                                                                                             >  2:723265/1‑132 (MQ=255)
                                                                                                                                                                    tAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                                                                                                                                             <  1:723265/132‑1 (MQ=255)
                                                                                                                                                                         ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaat                                                                                                                                                                                                      >  1:614551/1‑86 (MQ=255)
                                                                                                                                                                         ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaat                                                                                                                                                                                                      <  2:614551/86‑1 (MQ=255)
                                                                                                                                                                         ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCt                                                                                                          >  2:1045502/1‑178 (MQ=255)
                                                                                                                                                                         ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCt                                                                                                          <  1:1045502/178‑1 (MQ=255)
                                                                                                                                                                          tCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTCGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTtggtgg                                                                                                                      >  1:2079112/1‑165 (MQ=255)
                                                                                                                                                                                  tgtTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                                                                                                                                               >  1:575935/1‑68 (MQ=255)
                                                                                                                                                                                       tCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGt                               <  1:1636096/239‑1 (MQ=255)
                                                                                                                                                                                          cttcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       >  1:1189291/1‑52 (MQ=255)
                                                                                                                                                                                          cttcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttg                                                                                                                                                                                                                       <  2:1189291/53‑2 (MQ=255)
                                                                                                                                                                                          cttcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTa                                                                               <  1:1520075/188‑1 (MQ=255)
                                                                                                                                                                                          cttcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTaaa                                                  <  2:1004776/217‑1 (MQ=255)
                                                                                                                                                                                          cttcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTaaa                                                  >  1:1004776/1‑217 (MQ=255)
                                                                                                                                                                                            tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTa                                                                    <  2:198905/197‑1 (MQ=255)
                                                                                                                                                                                             cttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATg                         <  1:327919/239‑1 (MQ=255)
                                                                                                                                                                                              ttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATg                                                            >  2:1512228/1‑203 (MQ=255)
                                                                                                                                                                                                   atataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGt                   <  2:2305357/239‑1 (MQ=255)
                                                                                                                                                                                                     ataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGt                   <  1:19794/237‑1 (MQ=255)
                                                                                                                                                                                                     ataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGt                   >  2:19794/1‑237 (MQ=255)
                                                                                                                                                                                                       aaCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTa                                                                                                                                                                                              >  2:2114555/1‑64 (MQ=255)
                                                                                                                                                                                                       aaCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTa                                                                                                                                                                                              >  2:93980/1‑64 (MQ=255)
                                                                                                                                                                                                       aaCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTa                                                                                                                                                                                              <  1:93980/64‑1 (MQ=255)
                                                                                                                                                                                                       aaCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTa                                                                                                                                                                                              <  1:2114555/64‑1 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGa                                                                                                                                                                                                 >  1:2190394/1‑56 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAAt                                                                                                                                        <  2:196408/113‑1 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAAt                                                                                                                                        >  1:196408/1‑113 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTc                                                                                         <  2:906685/160‑1 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTa                                                                                     >  1:846091/1‑164 (MQ=255)
                                                                                                                                                                                                            ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTa                                                                                     <  2:846091/164‑1 (MQ=255)
                                                                                                                                                                                                              aTTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTa                                                                                                                                                 <  2:1668010/102‑1 (MQ=255)
                                                                                                                                                                                                              aTTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTa                                                                                                                                                 >  1:1668010/1‑102 (MQ=255)
                                                                                                                                                                                                                   cTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCt                                                                                                          >  1:1375554/1‑136 (MQ=255)
                                                                                                                                                                                                                    tCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTtggtgg                                                                                                                      >  2:1269993/1‑123 (MQ=255)
                                                                                                                                                                                                                    tCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGTAACGCTTATACCAATGa  >  2:1384503/1‑239 (MQ=255)
                                                                                                                                                                                                                      aaaCGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTAt                                                                                                                                                                                             >  2:942258/1‑50 (MQ=255)
                                                                                                                                                                                                                      aaaCGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTAt                                                                                                                                                                                             <  1:942258/50‑1 (MQ=255)
                                                                                                                                                                                                                                cGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGt                               <  2:555535/198‑1 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTAcc                                                                                                            >  1:2064886/1‑119 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTAcc                                                                                                            <  2:2064886/119‑1 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTAt                                                                              <  2:556307/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTAt                                                                              >  1:556307/1‑149 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGTaa                 >  1:524337/1‑210 (MQ=255)
                                                                                                                                                                                                                                  gCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGTaa                 <  2:524337/210‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            
GATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTACTCATATGTTCTTTTAAACATTCATCGTGGGGATGGTAAAATGAATAGTAACGCTTATACCAATGA  >  NC_000913/3798421‑3798872

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