Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,093 G→A 18.0% 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,0930GA18.0% 263.5 / 45.5 111V563M (GTG→ATG) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (47/44);  new base A (13/7);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 111 total aligned reads displayed.

TCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3542001‑3542189
                                                                                            |                                                                                                
tctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGa                                                                                           <  1:719354/100‑1 (MQ=255)
tctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGa                                                                                           >  1:142412/1‑100 (MQ=255)
 ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                                                          <  2:1985774/100‑1 (MQ=255)
 ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                                                          <  2:2130185/100‑1 (MQ=255)
  tcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                                                         <  1:2093055/100‑1 (MQ=255)
   cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGa                                                                                        >  2:945767/1‑100 (MQ=255)
   cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGa                                                                                        <  2:736244/100‑1 (MQ=255)
     cTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATa                                                                                      >  2:2519708/1‑100 (MQ=255)
       gAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAAc                                                                                    <  1:2274235/100‑1 (MQ=255)
         gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTg                                                                                  >  2:1437473/1‑100 (MQ=255)
          caggAAAATCGTCAATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTgg                                                                                 >  1:1754618/3‑100 (MQ=255)
           gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggc                                                                                <  1:2260546/100‑1 (MQ=255)
            ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggca                                                                               >  2:2455757/1‑100 (MQ=255)
              aaaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcagg                                                                             <  2:189493/100‑1 (MQ=255)
               aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcaggc                                                                            <  2:235735/100‑1 (MQ=255)
               aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcaggc                                                                            <  1:933842/100‑1 (MQ=255)
                   cgtcgATAACATCCACGACTCGGGGCTGGCGTCCGTCAGCCGGGTGATCACACCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACg                                                                        >  2:1882014/1‑100 (MQ=255)
                   cgtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACg                                                                        >  1:756184/1‑100 (MQ=255)
                       gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                    >  1:19960/1‑100 (MQ=255)
                       gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                    >  2:422645/1‑100 (MQ=255)
                         tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                  <  1:1743317/100‑1 (MQ=255)
                           aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                >  1:173365/1‑100 (MQ=255)
                           aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                >  2:2250579/1‑100 (MQ=255)
                           aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                >  2:695684/1‑100 (MQ=255)
                            cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  2:157086/1‑100 (MQ=255)
                            cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               <  1:2163121/100‑1 (MQ=255)
                            cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  2:2137995/1‑100 (MQ=255)
                            cATCAACGACTCGGCGCTGGCGTCCGTCAGCCAGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  2:1456186/1‑100 (MQ=255)
                            cATCAACAACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  1:1370048/1‑100 (MQ=255)
                               cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            <  1:2449620/100‑1 (MQ=255)
                               cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGtct                                                            <  1:1750682/100‑3 (MQ=255)
                               cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            <  1:1543933/100‑1 (MQ=255)
                                aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                           >  1:2552362/1‑100 (MQ=255)
                                 aCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTAc                                                          <  1:1277987/100‑1 (MQ=255)
                                   gACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAg                                                        >  2:2135032/1‑100 (MQ=255)
                                   gACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAg                                                        <  1:1114613/100‑1 (MQ=255)
                                     cTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGt                                                      >  2:1877185/1‑100 (MQ=255)
                                     cTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGt                                                      <  1:163536/100‑1 (MQ=255)
                                      tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     <  1:994299/100‑1 (MQ=255)
                                      tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     <  1:1820663/100‑1 (MQ=255)
                                       cGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGc                                                    >  1:2044155/1‑100 (MQ=255)
                                       cGGCGCTGGCGTCCGTCAGCCGGGAGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGc                                                    >  2:1634390/1‑100 (MQ=255)
                                        ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                   >  1:1323206/1‑100 (MQ=255)
                                        ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                   >  1:2195403/1‑100 (MQ=255)
                                         gcgcTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCa                                                  >  1:2412167/1‑100 (MQ=255)
                                             tGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc                                              <  1:1524526/100‑1 (MQ=255)
                                              ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             >  1:1006628/1‑100 (MQ=255)
                                               gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             <  1:667420/99‑1 (MQ=255)
                                                cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                           <  1:1547676/100‑1 (MQ=255)
                                                  tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCTGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaag                                         <  2:1604315/100‑1 (MQ=255)
                                                   ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaaga                                        <  2:173363/100‑1 (MQ=255)
                                                    cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  1:2197868/1‑100 (MQ=255)
                                                    cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  1:1265788/1‑100 (MQ=255)
                                                     gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                      >  1:330755/1‑100 (MQ=255)
                                                     gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                      <  2:1689935/100‑1 (MQ=255)
                                                     gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                      <  2:770928/100‑1 (MQ=255)
                                                     aagaccggggTGATCACCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGt                                     >  2:1857716/8‑100 (MQ=255)
                                                      agaccggggTGATCACCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                    >  2:344349/7‑100 (MQ=255)
                                                        aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                   >  1:1272996/1‑100 (MQ=255)
                                                          ccGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa                                 >  2:2036670/1‑100 (MQ=255)
                                                            gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGt                               >  1:436938/1‑100 (MQ=255)
                                                            gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGt                               >  2:1870760/1‑100 (MQ=255)
                                                             ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                              <  2:1973196/100‑1 (MQ=255)
                                                             ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                              >  2:2559500/1‑100 (MQ=255)
                                                               tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                            >  1:1491416/1‑100 (MQ=255)
                                                               tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                            <  2:574159/100‑1 (MQ=255)
                                                                gATCACCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                          <  1:344353/100‑1 (MQ=255)
                                                                 aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                          >  1:1362838/1‑100 (MQ=255)
                                                                 aTCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                          <  1:476014/100‑1 (MQ=255)
                                                                  tCACCCCGGTCTTCAAGACAATTGGCATGCATGAAGATAACTGGCAGGCACCGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGtag                             >  2:1407761/1‑94 (MQ=255)
                                                                   cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                        >  2:1436627/1‑100 (MQ=255)
                                                                   cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                        >  2:88531/1‑100 (MQ=255)
                                                                   cACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                        >  1:1401495/1‑100 (MQ=255)
                                                                    aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGc                       <  1:1495799/100‑1 (MQ=255)
                                                                     ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGTAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                      >  2:1082540/1‑100 (MQ=255)
                                                                     ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                      >  2:1766914/1‑100 (MQ=255)
                                                                       ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    >  2:515614/1‑100 (MQ=255)
                                                                       ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    >  1:182416/1‑100 (MQ=255)
                                                                       ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    <  1:2198671/100‑1 (MQ=255)
                                                                         ggTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    <  2:131424/98‑1 (MQ=255)
                                                                         ggTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                  <  2:1002987/100‑1 (MQ=255)
                                                                            cTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTAGGTACGCTCAACAcc               >  1:2567147/1‑100 (MQ=255)
                                                                            cTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc               >  1:2237592/1‑100 (MQ=255)
                                                                               naaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc            >  2:2004076/2‑100 (MQ=255)
                                                                               cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc            >  2:2393802/1‑100 (MQ=255)
                                                                                aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAGGAAGTGGTAGTGGGTACGCTCAACACCct             >  2:2169053/1‑98 (MQ=255)
                                                                                aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct           <  2:196225/100‑1 (MQ=255)
                                                                                aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct           <  2:880609/100‑1 (MQ=255)
                                                                                 aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa          >  1:1696337/1‑100 (MQ=255)
                                                                                 aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa          >  1:66123/1‑100 (MQ=255)
                                                                                 aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAAGCTCAACACCCTCTa          >  1:2561321/1‑100 (MQ=255)
                                                                                    cAATTGGCGTGCATGAAGATAACTGGCAGGCCACGGTTGGCCTGTGTACAGGTGACATGGCGAAAGAAGAGGGAGGGGGTACGATGAACAACCTATacaa       >  2:1294828/1‑99 (MQ=255)
                                                                                    cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac       <  2:625993/100‑1 (MQ=255)
                                                                                    cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac       <  1:2169065/100‑1 (MQ=255)
                                                                                     aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc      >  1:1613730/1‑100 (MQ=255)
                                                                                     aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACCCTCAACACCCTCTACAcc      <  1:2474864/100‑1 (MQ=255)
                                                                                        tGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa   <  2:1763423/100‑1 (MQ=255)
                                                                                         ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  >  1:1815961/1‑100 (MQ=255)
                                                                                         ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  <  1:766842/100‑1 (MQ=255)
                                                                                         ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  <  2:329000/100‑1 (MQ=255)
                                                                                            |                                                                                                
TCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3542001‑3542189

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