Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,914,552 G→A P217S (CCG→TCG)  glmU ← fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,914,5520GA90.5% 323.8 / 2.2 106P217S (CCG→TCG) glmUfused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base A (41/54);  total (50/56)
Fisher's exact test for biased strand distribution p-value = 4.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTG  >  NC_000913/3914406‑3914697
                                                                                                                                                  |                                                                                                                                                 
cGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATg                                                                                                                                                >  1:568850/1‑150 (MQ=255)
cGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAAAGTTGCGGATg                                                                                                                                                >  2:1234511/1‑150 (MQ=255)
   aCCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACg                                                                                                                                            <  2:528492/151‑1 (MQ=255)
   aCCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAAc                                                                                                                                             <  1:1176220/150‑1 (MQ=255)
    ccACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACg                                                                                                                                            <  2:92505/150‑1 (MQ=255)
      aCGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTATTCCCCCCCTCTACCTCCCTTACACCTTGCCGaacaccgg                                                                                                                                           >  1:693243/1‑142 (MQ=255)
      aCGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGATATTCACGCCCTCTACTTCGCTTTAAAGTTGCGAATGAACgg                                                                                                                                           >  2:129161/1‑149 (MQ=255)
                aaaCGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCAc                                                                                                                               <  2:1707392/151‑1 (MQ=255)
                  aCGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGc                                                                                                                             <  2:1333558/151‑1 (MQ=255)
                     cgcTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGGCGATTTCACGccc                                                                                                                           >  1:60365/1‑150 (MQ=255)
                          gATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTg                                                                                                                     <  1:1696811/151‑1 (MQ=255)
                            tCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCCCTTCAACGTTGCGAATGAACGGCGACGATTTTACCCCCTTCCTGAt                                                                                                                   >  1:30126/1‑151 (MQ=255)
                             cgcgCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGAAGAACGGCGACGATTTCACGCCCTTcc                                                                                                                       >  1:487499/1‑146 (MQ=255)
                             cgcgCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGAt                                                                                                                   <  2:1214316/150‑1 (MQ=255)
                             cgcgCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCGCTTAAACGTTGCGAATGAAAGGCGGCGATTTCACGCCCTTCCTgtt                                                                                                                   >  1:110293/1‑148 (MQ=255)
                              gcgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATAc                                                                                                                 <  1:1630189/151‑1 (MQ=255)
                               cgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATAct                                                                                                                <  1:1593596/151‑2 (MQ=255)
                                  agcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTTGCTTAAAAGTTGCGAATGAACGGCGACGATTTTACGCCCTTCCTGATAcccc                                                                                                              >  2:259724/1‑150 (MQ=255)
                                   gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAg                                                                                                            <  2:1059842/151‑1 (MQ=255)
                                    caTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgc                                                                                                           <  1:1619397/151‑1 (MQ=255)
                                     aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcg                                                                                                          <  2:358217/151‑1 (MQ=255)
                                     aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTCCCCCCCTTCCTGATACGCCAgc                                                                                                           >  1:1537851/1‑150 (MQ=255)
                                        aCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcg                                                                                                          >  2:315028/1‑148 (MQ=255)
                                           ccTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataa                                                                                                    <  2:712074/151‑1 (MQ=255)
                                             tGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCCGCGCaataa                                                                                                    >  2:907046/1‑149 (MQ=255)
                                              gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataat                                                                                                   <  1:160561/149‑1 (MQ=255)
                                              gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataat                                                                                                   <  2:1437311/149‑1 (MQ=255)
                                              gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGt                                                                                                 <  1:587899/151‑1 (MQ=255)
                                                 aaCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCAATAATGTCgg                                                                                              >  1:1624415/1‑151 (MQ=255)
                                                 aaCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAAAGTTGCGGATGAACGGCGACGATTTTACGCCCCTCCTGATACGCCAGCGCAATAATGTCgg                                                                                              >  1:472260/1‑151 (MQ=255)
                                                 aaCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCCCTTAAACGTTGCGAATGAACGGGGGCGATTTCACGCCCTTCCTTATACGCCCGCGCAATAATGTCgg                                                                                              >  2:822722/1‑151 (MQ=255)
                                                 aaCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCgg                                                                                              >  2:1115628/1‑151 (MQ=255)
                                                   cagcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTTAACGTTGCGGATGGAAGGCGACGATTTCACCCCCCTCCCTATACCCCAGCGCCATAATGGCgg                                                                                              >  1:590499/1‑149 (MQ=255)
                                                    agcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTACTTTGCTTAAAAGTTGCGAATGAACGGCGGCGATTTTACGCCCTTCCTGATACACCAGCGCATTAATGTCGGTgg                                                                                           >  1:1524054/1‑150 (MQ=255)
                                                     gcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCAGGCCCTCTACATTCCTTTAACGGTGCGAATGAAAGGAGACGATTTTCCGCCCTTTCTGTAACGCATGCGCCAAAATTTCGGTGAt                                                                                          >  1:1213822/1‑151 (MQ=255)
                                                     gcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGAt                                                                                          <  1:1565829/151‑1 (MQ=255)
                                                      cagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATg                                                                                         <  1:53269/151‑1 (MQ=255)
                                                      cagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATg                                                                                         <  1:556081/151‑1 (MQ=255)
                                                            ttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTAt                                                                                   >  2:75456/1‑151 (MQ=255)
                                                             tCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGGTGTAGTAtt                                                                                  >  1:262221/1‑151 (MQ=255)
                                                              cAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGGCGATTTCACGCCCTTCCTGATACGCCCGCGCAATAATGTCGGGGGTGTAGTAtt                                                                                  >  1:947411/1‑150 (MQ=255)
                                                              cAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAtt                                                                                  <  2:19105/150‑1 (MQ=255)
                                                               aGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTAtt                                                                                  >  1:26309/1‑149 (MQ=255)
                                                                  cTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACCGCTTCTAGATCCCTTCAATCTTGGGTAATTACGGCGAGGATTTGACCCCCCTCCTCGTACGCCAGCGCAACAATGACGGGGGTTTAGTATTCGccc                                                                             >  1:184629/1‑151 (MQ=255)
                                                                   tGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTTGCTTAAACGTTGCGAATGAACGGCGACGAATTCACGCCCCTCCTGATACGCCCGCGCAATAATGTCGGGGGTGTAGTAtttcccc                                                                             >  2:1656744/1‑145 (MQ=255)
                                                                     ttCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGGCGGTTTCACGCCCTTCCTGATACGCCCGCACAATAATGTCGGTGGTGTAGTATTTGccc                                                                             >  1:575201/1‑148 (MQ=255)
                                                                      tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGa                                                                          <  2:582049/150‑1 (MQ=255)
                                                                      tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGa                                                                          <  1:783778/150‑1 (MQ=255)
                                                                        ggACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGATATTCACGCCTTCTACTTCCCTTAAACGTTGCGGATGAAACGCGAGCATTTCCCGCCCTTCCTTATACCCCAACGCAATAATTTCGGTGGTGTAGGATTCGCCCTggg                                                                         >  1:1414695/1‑147 (MQ=255)
                                                                          aCTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCatt                                                                     <  1:1235186/151‑1 (MQ=255)
                                                                                aaaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGGGAATGAACGGCGACGATTTTTCGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGGAATCGCCCCTAGCattatt                                                                  >  1:938097/1‑148 (MQ=255)
                                                                                aaaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgt                                                                <  2:1661850/150‑1 (MQ=255)
                                                                                 aaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgttg                                                              <  2:104456/151‑1 (MQ=255)
                                                                                 aaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgttg                                                              <  2:1062596/151‑1 (MQ=255)
                                                                                   cacGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGt                                                            <  1:1565861/151‑1 (MQ=255)
                                                                                     cGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCa                                                          <  1:1413238/151‑1 (MQ=255)
                                                                                      gCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAg                                                         <  2:468949/151‑1 (MQ=255)
                                                                                      gCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAg                                                         <  1:1357148/151‑1 (MQ=255)
                                                                                      gCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAg                                                         <  2:1192142/151‑1 (MQ=255)
                                                                                         ccAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCtt                                                      <  2:845302/151‑1 (MQ=255)
                                                                                         ccAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCt                                                       <  1:295302/150‑1 (MQ=255)
                                                                                          cAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTc                                                     <  2:494329/151‑1 (MQ=255)
                                                                                          cAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTc                                                     <  1:456601/151‑1 (MQ=255)
                                                                                            gACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgc                                                   <  1:1301534/151‑1 (MQ=255)
                                                                                              cGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTTGCCCTGAGCATTATTGTTGGGCAGCTTCgcc                                                  >  1:943322/1‑150 (MQ=255)
                                                                                              cGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgcca                                                 <  1:474069/151‑1 (MQ=255)
                                                                                                ggAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagc                                               <  1:1458203/151‑1 (MQ=255)
                                                                                                  agagTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagcc                                              >  1:1512377/1‑150 (MQ=255)
                                                                                                  agagTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagcca                                             <  2:1162597/151‑1 (MQ=255)
                                                                                                  agagTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagcca                                             <  2:108500/151‑1 (MQ=255)
                                                                                                   gagTTGCAGGCGGTTATTTACGCCCTCTACTTCGCTTAAACGTTGCGAAAGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGCCCTGAGCCTTATTGTTGGTCAGCTTCgccagcc                                              >  1:360578/1‑149 (MQ=255)
                                                                                                   gagTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagcca                                             <  1:599818/150‑1 (MQ=255)
                                                                                                    agTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAAAGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGGATTCCCCCTTAGCATTATTGTTGGTCAGCTTcgcccgcc                                              >  1:1497426/1‑148 (MQ=255)
                                                                                                       tGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGt                                         <  2:377503/150‑1 (MQ=255)
                                                                                                         cAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTc                                      <  2:969503/151‑1 (MQ=255)
                                                                                                               gTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCt                                <  1:1659573/151‑1 (MQ=255)
                                                                                                               gTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCt                                <  1:503433/151‑1 (MQ=255)
                                                                                                                ttATTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTg                               >  1:287974/1‑151 (MQ=255)
                                                                                                                  aTTCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTgcg                             >  2:607394/1‑151 (MQ=255)
                                                                                                                    tCACGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGcc                           >  1:23763/1‑151 (MQ=255)
                                                                                                                      aCGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGt                         >  2:1634641/1‑151 (MQ=255)
                                                                                                                      aCGCCTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGt                         >  2:669370/1‑151 (MQ=255)
                                                                                                                         ctttCTACTTAGCTTAAACGTTGCGAATGAGCGGCGACGATTTCACGCCCTTCCTGATACCCCAGAGCAATAATGTCCGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTACGCCAGCCAGCGTTTCATATCTGCGCCGGTgg                      >  1:1633830/3‑151 (MQ=255)
                                                                                                                         ccTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTgg                      >  2:1333764/1‑151 (MQ=255)
                                                                                                                         ccTTCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTgg                      >  2:629524/1‑151 (MQ=255)
                                                                                                                           ttCTACTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGc                     >  2:1574678/1‑150 (MQ=255)
                                                                                                                               aCTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATc                 >  1:647916/1‑150 (MQ=255)
                                                                                                                                cTTCGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAAt                  <  1:1337561/148‑1 (MQ=255)
                                                                                                                                   cGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt            <  2:1430961/151‑1 (MQ=255)
                                                                                                                                   cGCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt            <  1:1703945/151‑1 (MQ=255)
                                                                                                                                    gCTTAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATg           >  2:440835/1‑151 (MQ=255)
                                                                                                                                      ttAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGcc         >  1:872297/1‑151 (MQ=255)
                                                                                                                                       tAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGcc         >  1:853609/1‑150 (MQ=255)
                                                                                                                                       tAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGc          <  2:60369/149‑1 (MQ=255)
                                                                                                                                       tAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGc          <  1:906978/149‑1 (MQ=255)
                                                                                                                                       tAAACGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGc          <  2:778301/149‑1 (MQ=255)
                                                                                                                                        aaaCGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCgg       >  2:1696027/1‑151 (MQ=255)
                                                                                                                                        aaaCGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCgg       <  1:651944/151‑1 (MQ=255)
                                                                                                                                        aaaCGTTGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCg        <  1:1158289/150‑1 (MQ=255)
                                                                                                                                             ttGCGAATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTg  >  2:1663488/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTG  >  NC_000913/3914406‑3914697

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