Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,395 G→T 5.2% C98* (TGC→TGA fucA ← L‑fuculose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,3950GT 5.2% 317.4 / 2.3 94C98* (TGC→TGAfucAL‑fuculose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base T (2/4);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTT  >  NC_000913/2933246‑2933543
                                                                                                                                                     |                                                                                                                                                    
ttATGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTg                                                                                                                                                      >  2:1412415/1‑150 (MQ=255)
 tACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGc                                                                                                                                                     <  2:635384/150‑1 (MQ=255)
  aCGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCa                                                                                                                                                    <  2:858775/150‑1 (MQ=255)
    gATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAAt                                                                                                                                                  <  2:155455/150‑1 (MQ=255)
     aTTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATg                                                                                                                                                 <  1:610962/150‑1 (MQ=255)
      tttttGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGa                                                                                                                                                >  2:212418/1‑150 (MQ=255)
        tttGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCAc                                                                                                                                              >  1:775737/1‑150 (MQ=255)
        tttGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCAc                                                                                                                                              >  2:923384/1‑150 (MQ=255)
        tttGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                                                              <  1:937302/150‑1 (MQ=255)
           gagaGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGc                                                                                                                                           >  2:982901/1‑150 (MQ=255)
             gagCCAGCGCAACATGTTCATAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGATTTACCGCCAGCCGCCGAAATCATGTTGTGAATAGCGGTGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTTCACGGCag                                                                                                                                         <  2:285227/150‑2 (MQ=255)
             gagCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGAGCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                                                                         <  2:519590/150‑2 (MQ=255)
              agCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATg                                                                                                                                        <  2:784491/150‑1 (MQ=255)
               gCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtt                                                                                                                                       <  1:141417/150‑4 (MQ=255)
                ccAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGGTCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                                                                      <  1:941339/150‑1 (MQ=255)
                ccAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                                                                      >  2:680138/1‑150 (MQ=255)
                 cAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAtt                                                                                                                                     >  1:1595130/1‑150 (MQ=255)
                  aGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcc                                                                                                                                    <  1:923384/150‑7 (MQ=255)
                  aGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                                                                    <  1:232886/150‑1 (MQ=255)
                     gcAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGACCGGCATGATTGTGa                                                                                                                                 >  2:1583037/1‑150 (MQ=255)
                            gTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                                                                          >  1:1423750/1‑150 (MQ=255)
                            gTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                                                                          <  1:219229/150‑1 (MQ=255)
                              tCAGAAAGTTGGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        >  2:1327380/1‑150 (MQ=255)
                              tCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        <  1:242973/150‑1 (MQ=255)
                              tCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        <  2:998642/150‑1 (MQ=255)
                              tCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                                                                        >  1:1169568/1‑150 (MQ=255)
                               cAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                                                                       >  1:751343/1‑150 (MQ=255)
                                 gAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                                                                     <  1:63562/150‑1 (MQ=255)
                                  aaaGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                                                                    <  1:1067331/150‑1 (MQ=255)
                                      ttCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATc                                                                                                                >  1:887351/1‑150 (MQ=255)
                                      ttCGCGTGTTCAAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATc                                                                                                                <  1:934516/150‑1 (MQ=255)
                                         gcgTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCgg                                                                                                             <  1:374430/150‑1 (MQ=255)
                                              ttCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                                                                        >  1:1526023/1‑150 (MQ=255)
                                              ttCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                                                                        <  2:711930/150‑1 (MQ=255)
                                               tCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCt                                                                                                       >  1:1166465/1‑150 (MQ=255)
                                                 cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCttt                                                                                                     >  2:40204/1‑150 (MQ=255)
                                                 cAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCttt                                                                                                     >  2:844212/1‑150 (MQ=255)
                                                  aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTg                                                                                                    >  2:49780/1‑150 (MQ=255)
                                                  aaaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTg                                                                                                    <  1:111374/150‑1 (MQ=255)
                                                    aGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAt                                                                                                  <  1:982901/150‑1 (MQ=255)
                                                     ggTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATa                                                                                                 <  2:1423750/150‑1 (MQ=255)
                                                      gTCGCATAAGGCGCGCACGGAATACAATTACCGCCACCCGCCGCAATCATGTAGTGAATCGCGGGTATCGATCGGTTAAGAATGTAAACTGCCGTGCAATGAACGGCATGATTGTGAACACCCGCGTTGGCATCCGGTCTGCTTTGATAg                                                                                                <  2:1135505/150‑1 (MQ=255)
                                                        cGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGc                                                                                              <  1:1012003/150‑1 (MQ=255)
                                                        cGCATAAGGCGCGCAAGGAACAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGc                                                                                              <  2:964121/150‑1 (MQ=255)
                                                         gCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCt                                                                                             >  1:1343101/1‑150 (MQ=255)
                                                          cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                                                            >  1:1600004/1‑150 (MQ=255)
                                                                  gcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                                                                    >  2:367097/1‑150 (MQ=255)
                                                                  gcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                                                                    >  2:1320212/1‑150 (MQ=255)
                                                                   cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGa                                                                                   >  2:811357/1‑150 (MQ=255)
                                                                        ggAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGc                                                                              >  2:1165415/1‑150 (MQ=255)
                                                                         gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                                                             >  1:1338178/1‑150 (MQ=255)
                                                                            tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCAtt                                                                          >  2:1375804/1‑150 (MQ=255)
                                                                              gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                                                              >  1:10715/1‑80 (MQ=255)
                                                                              gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                                                              <  2:10715/80‑1 (MQ=255)
                                                                                aTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                      <  1:844212/150‑1 (MQ=255)
                                                                                aTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                                                                      >  2:151757/1‑150 (MQ=255)
                                                                                 ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                                                                     >  1:1094498/1‑150 (MQ=255)
                                                                                   aCCGCCAGCCGCCGCAATCATGTAGTGAATCGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGGGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGa                                                                   <  2:359279/150‑1 (MQ=255)
                                                                                   aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGa                                                                   >  1:523598/1‑150 (MQ=255)
                                                                                   aCCGCCACCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGCCTGCCATATGCAAACGCCATTCGCTTGa                                                                   >  2:394121/1‑150 (MQ=255)
                                                                                     cGCCAGCCTCAGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                                                                 <  1:606432/150‑1 (MQ=255)
                                                                                     cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                                                                 >  2:1254225/1‑150 (MQ=255)
                                                                                     cGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                                                                 <  1:49780/150‑1 (MQ=255)
                                                                                      gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                                                                >  1:677889/1‑150 (MQ=255)
                                                                                      gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                                                                >  2:1154324/1‑150 (MQ=255)
                                                                                       ccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                                                               <  2:677889/150‑1 (MQ=255)
                                                                                          gccgccGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGc                                                            >  1:685926/1‑150 (MQ=255)
                                                                                          gccgccGCAATCATGTAGTGAATAGCGGGGATCGAGCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTGGAGGGGAGc                                                            <  2:1253632/150‑1 (MQ=255)
                                                                                           ccgccgCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAATGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCt                                                           >  2:809456/1‑150 (MQ=255)
                                                                                                  aaTCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc                                                    <  2:1321311/150‑1 (MQ=255)
                                                                                                   aTCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc                                                   <  1:367097/150‑1 (MQ=255)
                                                                                                    tCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct                                                  <  2:1166465/150‑1 (MQ=255)
                                                                                                    acATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGTGGAGCTTtccttcct                                                  <  2:337807/149‑1 (MQ=255)
                                                                                                         tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt                                             >  1:1387902/1‑150 (MQ=255)
                                                                                                         tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt                                             >  1:1004381/1‑150 (MQ=255)
                                                                                                         tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt                                             <  2:1595130/150‑1 (MQ=255)
                                                                                                          cgTGAATAGCGGGGATCGCTCGGTTAAGACTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGtt                                            <  2:1286300/149‑1 (MQ=255)
                                                                                                                tAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGt                                      <  1:1165415/150‑1 (MQ=255)
                                                                                                                  gCGGGGATCGATCGGTTAAGAATGGAAACTGGCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTg                                    <  1:1618122/150‑1 (MQ=255)
                                                                                                                       gatcgatcGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATc                               <  2:461950/150‑1 (MQ=255)
                                                                                                                          cgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGc                                   <  1:1507805/143‑1 (MQ=255)
                                                                                                                          cgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGc                                   >  2:1507805/1‑143 (MQ=255)
                                                                                                                          cgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAAt                            >  1:1468535/1‑150 (MQ=255)
                                                                                                                           gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATa                           <  2:1526023/150‑1 (MQ=255)
                                                                                                                                gTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGAc                      >  1:241492/1‑150 (MQ=255)
                                                                                                                                gTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCCTTGCCATCAATAAAGAc                      >  1:591933/1‑150 (MQ=255)
                                                                                                                                    aGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAata                  <  2:1169568/150‑1 (MQ=255)
                                                                                                                                     caaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAatat                 <  2:333307/149‑1 (MQ=255)
                                                                                                                                          gAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGAc            >  1:659458/1‑150 (MQ=255)
                                                                                                                                              cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCg        <  1:1327380/150‑1 (MQ=255)
                                                                                                                                              cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCg        >  1:1543713/1‑150 (MQ=255)
                                                                                                                                                  cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAg    <  1:1412415/150‑1 (MQ=255)
                                                                                                                                                  cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAg    >  2:981067/1‑150 (MQ=255)
                                                                                                                                                   gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGt   <  1:151757/150‑1 (MQ=255)
                                                                                                                                                    tGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGtt  <  2:1169326/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTACGATTTTTGAGAGCCAGCGCAACATGTTCAGAAAGTTCGCGTGTTCCAAAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTT  >  NC_000913/2933246‑2933543

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