Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,401 A→C 15.9% V96V (GTT→GTG fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,4010AC15.9% 391.6 / 60.0 156V96V (GTT→GTGfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (56/73);  new base C (10/15);  total (66/90)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATG  >  NC_000913/2933304‑2933499
                                                                                                 |                                                                                                  
cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                  >  1:840321/1‑100 (MQ=255)
  tAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACgg                                                                                                <  2:1383037/100‑1 (MQ=255)
  tAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACgg                                                                                                >  2:1520045/1‑100 (MQ=255)
    aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                              >  1:2384138/1‑100 (MQ=255)
     ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                             <  2:670633/100‑2 (MQ=255)
     ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                             <  1:1340649/100‑2 (MQ=255)
      gcgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagt                                                                                            <  2:2583747/100‑3 (MQ=255)
      gcgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATg                                                                                            <  2:2455638/100‑1 (MQ=255)
       ctcgcAAGGCATAGAATTACCGCCAGTCGCCGCAATCAAGTAGTGAATAGCGGGGATCGTTCGGTTAAGAATGGAAACTGCCGTGCAATGTACGACAt                                                                                             <  2:1988747/96‑1 (MQ=255)
       cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGa                                                                                           <  2:23209/100‑1 (MQ=255)
         cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAtt                                                                                         >  1:1426196/1‑100 (MQ=255)
          gtaaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                        <  2:469400/98‑1 (MQ=255)
          gcagataaTAGAATCACCGCCCGCCGCTCCCTTCATGTGGTGAAGAGCGGCGATGGATGGGTTACGAATCGAGGCTTCCGTGCAATGAACGGTAt                                                                                             <  2:3034759/89‑1 (MQ=255)
          gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcc                                                                                        <  1:67344/100‑7 (MQ=255)
          gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcc                                                                                        <  1:2810848/100‑7 (MQ=255)
          gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAtcg                                                                                        <  1:3542115/100‑3 (MQ=255)
          gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                        <  2:1153707/100‑1 (MQ=255)
           cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcca                                                                                       <  1:44513/100‑8 (MQ=38)
           cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                       <  1:510063/100‑1 (MQ=255)
           cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                       >  1:1757301/1‑100 (MQ=255)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                      <  1:2025616/100‑9 (MQ=38)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                      <  1:3503360/100‑9 (MQ=38)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                      <  1:978088/100‑9 (MQ=38)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                      <  2:2753627/100‑1 (MQ=255)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                      <  2:3774439/100‑1 (MQ=255)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                      <  2:536364/100‑1 (MQ=255)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                      <  2:3774585/100‑1 (MQ=255)
            aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                      <  2:3256892/100‑1 (MQ=255)
            aaGGAATAGAATCACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATTGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                      <  1:1335622/100‑9 (MQ=37)
             aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                     <  2:2429005/100‑1 (MQ=255)
              ggAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaa                                                                                    <  1:1813616/100‑1 (MQ=255)
               gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac                                                                                   <  2:3401189/100‑1 (MQ=255)
               gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac                                                                                   <  1:2343025/100‑1 (MQ=255)
                   aGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAcc                                                                               <  1:3307717/100‑1 (MQ=255)
                    gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                              >  2:3723088/1‑100 (MQ=255)
                      aTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCAGGGATCGATCGGTTAAGAATGGAAACAGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                            >  2:848543/1‑100 (MQ=255)
                      aTTACCGCCAGCCGCCGCAAACATGTAGTGAATAGCGGGGATCGTTCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                            >  1:680061/1‑100 (MQ=255)
                       ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                           >  1:544150/1‑100 (MQ=255)
                         aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                         <  1:1055799/100‑1 (MQ=255)
                          ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                        >  2:1905503/1‑100 (MQ=255)
                          ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                        <  2:1511217/100‑1 (MQ=255)
                            gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                      >  2:2638349/1‑100 (MQ=255)
                            gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                      >  2:3687965/1‑100 (MQ=255)
                            gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                      >  1:2416153/1‑100 (MQ=255)
                            gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                      <  1:39973/100‑1 (MQ=255)
                                 ccgccgCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAGCAACCGCGTTGGCATCggg                                                                 >  2:2755171/1‑100 (MQ=255)
                                 ccgccgCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCgg                                                                 <  1:3412621/100‑1 (MQ=255)
                                      gcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                            <  2:1764397/100‑1 (MQ=255)
                                      gcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                            <  2:532813/100‑1 (MQ=255)
                                      gcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                            >  2:2565280/1‑100 (MQ=255)
                                           naTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                       >  2:71737/2‑100 (MQ=255)
                                           cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                       <  1:65110/100‑1 (MQ=255)
                                            aTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                       <  1:2319941/99‑1 (MQ=255)
                                             tctggTGAGTAGGGGGGATCGATCGGTTGAGACTGCTTGCTGCCGTGGTATGTCCGGCATGATTGTCAACAACCGCGCTGTCAGCTGGTCTGCTTTtata                                                     <  2:2215627/96‑1 (MQ=255)
                                             tGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                       >  2:3052404/1‑98 (MQ=255)
                                             tGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATa                                                     <  2:3514083/100‑1 (MQ=255)
                                             tGTAGTGAATAGCGGGGACCGATCGGTTAAGAATGGCGACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTgct                                                      >  2:2044617/1‑97 (MQ=255)
                                              gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                                    >  1:1700349/1‑100 (MQ=255)
                                               tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                                    >  2:2071299/1‑99 (MQ=255)
                                                aGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGc                                                  <  1:2215793/100‑1 (MQ=255)
                                                aGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGc                                                  <  1:3733003/100‑1 (MQ=255)
                                                  tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                >  1:3255851/1‑100 (MQ=255)
                                                   gAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGc                                               >  1:2738016/1‑100 (MQ=255)
                                                     aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCa                                             >  2:66355/1‑100 (MQ=255)
                                                       aGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCata                                           <  1:1183319/100‑1 (MQ=255)
                                                        gCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCatat                                          >  1:612710/1‑100 (MQ=255)
                                                          ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                        >  1:3000098/1‑100 (MQ=255)
                                                              atcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                    <  2:2165504/100‑1 (MQ=255)
                                                                 gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                 <  1:1241072/100‑1 (MQ=255)
                                                                 gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                 <  2:1568800/100‑1 (MQ=255)
                                                                 gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCAATATGGAAACGcc                                 <  1:2565297/100‑1 (MQ=255)
                                                                     ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                             <  1:1236457/100‑1 (MQ=255)
                                                                     ggTTAAGAATGGAAACTGCCGTGCAATGAACGACATGATGGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                             <  2:2675685/100‑1 (MQ=255)
                                                                      gTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGCTAGGCTGCCATATGGAAACGCCATTCg                            >  2:2268134/1‑100 (MQ=255)
                                                                       ttAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGc                           <  1:1570918/100‑1 (MQ=255)
                                                                        tAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                          <  1:1905514/100‑1 (MQ=255)
                                                                         aaGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                         <  1:1026207/100‑1 (MQ=255)
                                                                          ngAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTg                        >  2:1199742/2‑100 (MQ=255)
                                                                           gAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGa                       >  1:1929016/1‑100 (MQ=255)
                                                                            aaTGGAAACTGCCGTGCAATGAACGGCATGTTTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                      <  2:624427/100‑1 (MQ=255)
                                                                            aaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                      <  2:1907929/100‑1 (MQ=255)
                                                                            aaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                      <  1:571012/100‑1 (MQ=255)
                                                                             aTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                     <  2:1795635/100‑1 (MQ=255)
                                                                              acaatcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                    >  1:400752/8‑100 (MQ=38)
                                                                               ggAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                   >  1:3708755/1‑100 (MQ=255)
                                                                               caatcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                   >  2:1384386/7‑100 (MQ=255)
                                                                                gAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                  >  2:3517306/1‑100 (MQ=255)
                                                                                gAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                  >  2:2521575/1‑100 (MQ=255)
                                                                                 aaaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg                 <  2:3415371/100‑1 (MQ=255)
                                                                                   catGCCGTTCATTGCACGGCATGATTGTGAATAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCt               >  1:997859/3‑100 (MQ=255)
                                                                                   catGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCt               >  1:2061566/3‑100 (MQ=255)
                                                                                    atGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt              >  1:2838602/2‑100 (MQ=255)
                                                                                      gCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttc            >  1:2225363/1‑100 (MQ=255)
                                                                                       ccGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcc           <  1:3340233/100‑1 (MQ=255)
                                                                                         gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt         <  1:2728452/100‑1 (MQ=255)
                                                                                           gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc       <  2:1972401/100‑1 (MQ=255)
                                                                                           gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc       <  2:2009359/100‑1 (MQ=255)
                                                                                            cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct      >  1:1191261/1‑100 (MQ=255)
                                                                                                gAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATg  <  1:1347311/100‑1 (MQ=255)
                                                                                                gAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATg  <  1:715093/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATG  >  NC_000913/2933304‑2933499

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