Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,093 G→A 69.2% V563M (GTG→ATG)  feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,0930GA69.2% 142.8 / 85.3 112V563M (GTG→ATG) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base A (28/49);  total (46/66)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 115 total aligned reads displayed.

AACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3541994‑3542189
                                                                                                   |                                                                                                
aaCAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCg                                                                                                  >  1:891945/1‑100 (MQ=255)
  cAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATg                                                                                                >  1:2329377/1‑100 (MQ=255)
   aGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGc                                                                                               <  2:1829151/100‑1 (MQ=255)
   aGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGc                                                                                               <  1:146967/100‑1 (MQ=255)
   aGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGc                                                                                               >  2:1335992/1‑100 (MQ=255)
    gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCa                                                                                              >  1:802858/1‑100 (MQ=255)
    gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                              >  1:2622544/1‑100 (MQ=255)
    gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                              <  2:1323785/100‑1 (MQ=255)
    gCTTCTCGCTGAGCGGGAAAATAGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCa                                                                                              <  1:2028563/100‑1 (MQ=255)
     cTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGgcatgcat                                                                                             <  1:1944445/100‑1 (MQ=255)
     cTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGgcatgcat                                                                                             >  2:1515984/1‑100 (MQ=255)
        ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                                                          >  2:307214/1‑100 (MQ=255)
        ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGaa                                                                                          <  2:840473/100‑1 (MQ=255)
          cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                        >  1:2324722/1‑100 (MQ=255)
           gCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGAt                                                                                       <  2:2077398/100‑1 (MQ=255)
           gCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGAt                                                                                       <  1:45531/100‑1 (MQ=255)
            cTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGACGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATa                                                                                      >  1:106228/1‑100 (MQ=255)
              gAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAAc                                                                                    >  2:2445696/1‑100 (MQ=255)
              gAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAAc                                                                                    <  1:445084/100‑1 (MQ=255)
              gAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGACTTCAAGCCAATTGGCATGCATGAAGATAAc                                                                                    >  1:341490/1‑100 (MQ=255)
               aGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                   <  1:1056445/100‑1 (MQ=255)
                gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAACCAATTGGCATGCATGAAGATAACTg                                                                                  >  1:1734979/1‑100 (MQ=255)
                 cGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTgg                                                                                 >  2:1459397/1‑100 (MQ=255)
                  gggAAAATCGTCGATAACATCAACGTCTCGGCGCTGGCGTCCGTCAGTCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggc                                                                                <  1:528002/100‑1 (MQ=255)
                   ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggca                                                                               >  1:1718906/1‑100 (MQ=255)
                   ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggca                                                                               <  2:1031947/100‑1 (MQ=255)
                   ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggca                                                                               <  1:88106/100‑1 (MQ=255)
                      aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcaggc                                                                            <  1:2030927/100‑1 (MQ=255)
                         tcgtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAAc                                                                         <  2:2446970/100‑1 (MQ=255)
                             cgATAACATCAACGACTCGGCGCTGGCGTCCGTCCGCCGGGTGATAACCCCAGTCTTCACGCCAATTGGCGTCCCAGAAGCTACCTggca                                                                               >  2:186900/1‑90 (MQ=255)
                              gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGtt                                                                     >  1:294321/1‑99 (MQ=255)
                              gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                    <  2:2630723/100‑1 (MQ=255)
                              gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTg                                                                    <  2:817286/100‑1 (MQ=255)
                               aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                   <  1:2336804/100‑1 (MQ=255)
                                tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                  <  2:611383/100‑1 (MQ=255)
                                tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                  <  1:2311421/100‑1 (MQ=255)
                                 aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                 <  2:288473/100‑1 (MQ=255)
                                 aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                 <  2:1204322/100‑1 (MQ=255)
                                 aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                 <  1:2629668/100‑1 (MQ=255)
                                  aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                >  1:1082860/1‑100 (MQ=255)
                                   cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               <  2:448691/100‑1 (MQ=255)
                                   cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  2:739950/1‑100 (MQ=255)
                                      cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTTTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            <  1:421701/100‑1 (MQ=255)
                                      cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            >  1:1345312/1‑100 (MQ=255)
                                      cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            <  1:1143205/100‑1 (MQ=255)
                                             tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     <  2:568222/100‑1 (MQ=255)
                                             tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     <  2:611603/100‑1 (MQ=255)
                                              cGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGc                                                    <  1:1977992/100‑1 (MQ=255)
                                               ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                   <  2:2239388/100‑1 (MQ=255)
                                                   cTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                               <  2:926812/100‑1 (MQ=255)
                                                     ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             <  2:2127670/100‑1 (MQ=255)
                                                     ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             <  2:753585/100‑1 (MQ=255)
                                                      gCGTCTGTCAGCCGGGTGATCACCCCGGTTTTCAAGCTTATTTGCTTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             <  1:877192/99‑1 (MQ=255)
                                                       cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                           >  1:1096251/1‑100 (MQ=255)
                                                        gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGTATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                          <  1:1734122/100‑1 (MQ=255)
                                                         tccgtcAGCCGGGTGATCACTCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaag                                         <  1:1973011/100‑1 (MQ=255)
                                                         tccgtcAGCCGGGTGATCACTCCGGTCTTCAAGCCAATTAGAGTGCAGGAAGATAACTGGCAGGCAACAGTTGGGCAGTCTACAGCTGACGTggagcaa                                          >  2:2425202/1‑94 (MQ=255)
                                                         tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaag                                         >  1:1583238/1‑100 (MQ=255)
                                                            gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                      >  2:245791/1‑100 (MQ=255)
                                                             tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATTCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGt                                     <  2:2281669/100‑1 (MQ=255)
                                                              cAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                    <  1:1694702/100‑1 (MQ=255)
                                                              cAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                    <  2:327547/100‑1 (MQ=255)
                                                               aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTCACAGGTGCCATGGCGAAAGAAGTgg                                   <  1:835468/100‑1 (MQ=255)
                                                                gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCCGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  2:135442/1‑95 (MQ=255)
                                                                gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                  >  1:1043906/1‑100 (MQ=255)
                                                                gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAATGGTTGGACTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                  >  2:1859578/1‑100 (MQ=255)
                                                                gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                  <  2:1584216/100‑1 (MQ=255)
                                                                   gggTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTggaa                                 >  2:656827/1‑96 (MQ=255)
                                                                    ggTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                              <  1:1591283/100‑1 (MQ=255)
                                                                      tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                            <  2:1096241/100‑1 (MQ=255)
                                                                      tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                            <  1:1336005/100‑1 (MQ=255)
                                                                       gATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                           >  2:42320/1‑100 (MQ=255)
                                                                        aTCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                          <  2:2050982/100‑1 (MQ=255)
                                                                           aCCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGc                       >  1:2401629/1‑100 (MQ=255)
                                                                             cccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                     <  1:2353392/100‑1 (MQ=255)
                                                                              ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGTCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    <  1:769524/100‑1 (MQ=255)
                                                                              ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    <  2:1206151/100‑1 (MQ=255)
                                                                                ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                  <  2:549426/100‑1 (MQ=255)
                                                                                ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                  <  2:2324693/100‑1 (MQ=255)
                                                                                ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                  <  2:611528/100‑1 (MQ=255)
                                                                                 gTCTTCAAGCCAGTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                 >  1:1737290/1‑100 (MQ=255)
                                                                                   cTTCAAGCCAATTGGCATGCATGATGATAACTGGCAGGCAACGGTTGGTCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc               <  1:202651/100‑1 (MQ=255)
                                                                                   cTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAAGCTCAaca                 >  1:1185494/1‑98 (MQ=255)
                                                                                    ttaacgCCCATCGGCATGCATGAAGAGAACTGGCAGGAAACGGTTGGGATGTTGAGAGGTGATATGGCGAATGACGTGGTAGTGGGTACGCTGTacac                <  2:2214669/93‑1 (MQ=255)
                                                                                    ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc              <  1:903933/100‑1 (MQ=255)
                                                                                     tCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct             <  1:307216/100‑1 (MQ=255)
                                                                                      cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCACCGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc            >  2:1411688/1‑100 (MQ=255)
                                                                                      cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc            <  2:2329348/100‑1 (MQ=255)
                                                                                       aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct           >  2:658354/1‑100 (MQ=255)
                                                                                       aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct           <  2:741325/100‑1 (MQ=255)
                                                                                       aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct           <  1:1431038/100‑1 (MQ=255)
                                                                                          ccAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGAATACAGGTGCCATGGCGAACGGAGTGGAAGTGGGTACGCTCAAAACCCTCTaca        >  2:425888/1‑100 (MQ=255)
                                                                                            aaTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc      >  1:463473/1‑100 (MQ=255)
                                                                                             aTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg     <  1:1927991/100‑1 (MQ=255)
                                                                                              ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCTTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc    <  1:1578068/100‑1 (MQ=255)
                                                                                              ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc    >  1:2392134/1‑100 (MQ=255)
                                                                                              ttGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc    <  1:691532/100‑1 (MQ=255)
                                                                                                ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  >  2:146123/1‑100 (MQ=255)
                                                                                                ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGATACGCTCAACATCCTCTACACCGCAg  >  2:1318048/1‑100 (MQ=255)
                                                                                                ggCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  >  1:1674195/1‑100 (MQ=255)
                                                                                                   |                                                                                                
AACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3541994‑3542189

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