Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,093 G→A 69.9% 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.9% 199.8 / 126.9 156V563M (GTG→ATG) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (45/64);  total (65/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 157 total aligned reads displayed.

CAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAA  >  NC_000913/3541996‑3542191
                                                                                                 |                                                                                                  
cAGCTTCTCGCTGAGCGGGAAAATCGTCGCTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATg                                                                                                  <  1:62269/100‑1 (MQ=255)
 aGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGc                                                                                                 <  2:3651323/100‑1 (MQ=255)
  gattCTCGCTGAGTGGGAAAATCGTCGATAAGATCAACGACTAGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                                >  2:2983537/3‑100 (MQ=255)
    ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGcatgcatg                                                                                              >  2:4315828/1‑100 (MQ=255)
      ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                                                            <  2:1185666/100‑1 (MQ=255)
      ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                                                            >  1:2596477/1‑100 (MQ=255)
       tcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAg                                                                                           >  1:703616/1‑100 (MQ=255)
        cGCTGAGCGGGAAACTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAg                                                                                           <  2:4202552/99‑1 (MQ=255)
         gCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCTCTGGCGTCCGTCAGCCGGGTGATCACACCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                          >  2:2076665/1‑99 (MQ=255)
             aGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                     >  2:2514756/1‑100 (MQ=255)
              gCGGGACAGGCGGCGATAACATCAACGACTGGGCTCTGGCGTCTGTCAGCCGGGTGATCACTTCGGTCTTCAAGCCGATTGGCATGCATGCAGATAACt                                                                                     <  1:490870/99‑1 (MQ=255)
                 ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATACCTggca                                                                                 >  1:47057/1‑100 (MQ=255)
                   aaaaTCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcagg                                                                               <  1:220917/100‑1 (MQ=255)
                   aaaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcagg                                                                               <  2:2620101/100‑1 (MQ=255)
                   aaaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcagg                                                                               <  1:1476225/100‑1 (MQ=255)
                      aTCGTCGATAACATCAACGTCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcaggc                                                                              <  1:1014658/98‑1 (MQ=255)
                      aTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCaa                                                                            >  2:4245202/1‑100 (MQ=255)
                      aTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCaa                                                                            >  1:1644667/1‑100 (MQ=255)
                       tcgtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAAc                                                                           <  2:3997926/100‑1 (MQ=255)
                       tcgtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAAc                                                                           <  1:4062804/100‑1 (MQ=255)
                        cgtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACg                                                                          >  1:1790760/1‑100 (MQ=255)
                          tctaTAACATCAACGTCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGt                                                                        <  2:502229/97‑1 (MQ=255)
                           cgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGtt                                                                       >  1:1223145/1‑100 (MQ=255)
                           cgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGtt                                                                       >  2:1509743/1‑100 (MQ=255)
                            gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                      <  2:2563516/100‑1 (MQ=255)
                              tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                    <  1:4151810/100‑1 (MQ=255)
                              tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                    <  2:2214272/100‑1 (MQ=255)
                               aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                   >  2:772192/1‑100 (MQ=255)
                                aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                  >  2:2399330/1‑100 (MQ=255)
                                aCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                  >  2:272965/1‑100 (MQ=255)
                                 cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGTTCACGCCGGTCTGCAAGCCAATTGGCATGCATTGAGATGACTGGGAGGCAACGGTTGGcc                                                                   <  1:3061552/98‑1 (MQ=255)
                                 cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                                 >  2:1288638/1‑100 (MQ=255)
                                 cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                                 >  2:1943164/1‑100 (MQ=255)
                                 cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                                 <  2:907334/100‑1 (MQ=255)
                                  aTCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGt                                                                >  2:1475266/1‑100 (MQ=255)
                                   tCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGtt                                                               <  2:1473341/100‑1 (MQ=255)
                                    cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                              <  2:2556131/100‑1 (MQ=255)
                                    cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                              <  1:436149/100‑1 (MQ=255)
                                    cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                              <  2:2462110/100‑1 (MQ=255)
                                     aaCGACTCGGCGCTTGCGTCCGTCAGCCGGGTGATCACCCCGGACTTCAAGCCAATTGGCATGCATGCAGATAACTGGCCGGCAACGGTTGGCCTGCTTa                                                             >  1:3399523/1‑100 (MQ=255)
                                     aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                             >  1:3020955/1‑100 (MQ=255)
                                       cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                           <  2:1588761/100‑1 (MQ=255)
                                         aCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                                                         <  1:3200374/100‑1 (MQ=255)
                                           tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                       <  2:3020973/100‑1 (MQ=255)
                                             ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                     >  2:2114130/1‑100 (MQ=255)
                                                 cTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                                 <  2:782621/100‑1 (MQ=255)
                                                    gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGa                                              >  1:801296/1‑100 (MQ=255)
                                                      gtccgtcAGCCGTGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                            <  2:1424568/100‑1 (MQ=255)
                                                      gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGACAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                            <  2:2632149/100‑1 (MQ=255)
                                                      gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                            <  2:3746649/100‑1 (MQ=255)
                                                      gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                            <  2:3971099/100‑1 (MQ=255)
                                                      gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                            <  2:950835/100‑1 (MQ=255)
                                                        ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATAAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaaga                                          <  2:2341445/100‑1 (MQ=255)
                                                         cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                         >  1:3539674/1‑100 (MQ=255)
                                                         cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaaga                                          >  1:1221919/1‑99 (MQ=255)
                                                            cAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                      <  1:891262/100‑1 (MQ=255)
                                                              gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                    <  2:589102/100‑1 (MQ=255)
                                                              gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                    >  2:1018529/1‑100 (MQ=255)
                                                                cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAg                                  <  2:2864495/100‑1 (MQ=255)
                                                                  ggTGATCACCCCGGTCTTCAAGCCAATTGTCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                                <  2:1415724/100‑1 (MQ=255)
                                                                  ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                                <  2:2295079/100‑1 (MQ=255)
                                                                   gTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTgg                               >  2:3801240/1‑100 (MQ=255)
                                                                    tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                              <  2:3826161/100‑1 (MQ=255)
                                                                    tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAGGTGGTAGTggg                              <  2:212706/100‑1 (MQ=255)
                                                                    tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                              <  2:703617/100‑1 (MQ=255)
                                                                     gATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGCTACCTGGCAGGCAACGGTTGGCCTGTTTACAGCTGCCATGGCGAAAGAAGTGGTAGTggc                              <  1:3510585/99‑2 (MQ=255)
                                                                      aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                            <  1:1375736/100‑1 (MQ=255)
                                                                       tCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAc                           >  1:374379/1‑100 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                        >  2:907222/1‑100 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                        >  1:2093963/1‑100 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                        >  1:2097705/1‑100 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                        <  2:3959050/100‑1 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                        >  1:2757373/1‑100 (MQ=255)
                                                                             cGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                      <  2:3485774/99‑1 (MQ=255)
                                                                             cGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                      <  2:375267/99‑1 (MQ=255)
                                                                              ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                    >  2:1335319/1‑100 (MQ=255)
                                                                               gTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                   <  1:3223517/100‑1 (MQ=255)
                                                                                tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                  <  1:882131/100‑1 (MQ=255)
                                                                                tCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                  <  1:1181684/100‑1 (MQ=255)
                                                                                 cTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc                 <  2:2596491/100‑1 (MQ=255)
                                                                                  ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                <  2:1644674/100‑1 (MQ=255)
                                                                                  ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGACAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                >  1:155114/1‑100 (MQ=255)
                                                                                   tCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct               <  2:2093970/100‑1 (MQ=255)
                                                                                    cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc              >  1:1892425/1‑100 (MQ=255)
                                                                                    cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGTTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc              >  1:2829270/1‑100 (MQ=255)
                                                                                    cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc              <  2:2859457/100‑1 (MQ=255)
                                                                                     aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct             <  2:3282449/100‑1 (MQ=255)
                                                                                     aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGATGGGCTGTTTACAGTTGCCATGGCGAAAGACGGGGCAGTGGGTACGCTCAACAcc                 >  1:3685182/1‑96 (MQ=255)
                                                                                      aGCCAATTGGCGTGCATGAAGATAACCGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa            >  2:2680492/1‑100 (MQ=255)
                                                                                      aGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa            <  2:310844/100‑1 (MQ=255)
                                                                                      aGCCAACTGGCATGCATGAAGCTAACTGGCACGCAACGGCTGGCCTGTTTCCAGGTACCACGGCGAAAGCAGTGGCAGCGGGCCCGCTCAACACCctct             >  1:3851358/1‑99 (MQ=255)
                                                                                       gCCACTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac           >  2:3070237/1‑100 (MQ=255)
                                                                                       gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac           >  2:1948746/1‑100 (MQ=255)
                                                                                       ccgatctGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac           <  2:4181559/94‑1 (MQ=255)
                                                                                          aaTTAGCGTGGGTGAAGGTAACTGGGAGGCAACGGTGGGCCCGTTTACAGGTGCCATGGCGAGAGGTGGGGTAGTGGGTGCGCTCGAGGCTCTCTacaca        <  1:2110511/100‑2 (MQ=255)
                                                                                           aTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg       >  1:7604/1‑100 (MQ=255)
                                                                                           aTTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCAAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg       >  1:477655/1‑100 (MQ=255)
                                                                                            ttGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc      >  1:2342889/1‑100 (MQ=255)
                                                                                              ggCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg    >  1:2689795/1‑100 (MQ=255)
                                                                                                cGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAAGCTCAACACCCTCTACACCGCAGaa  >  1:1987331/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAA  >  NC_000913/3541996‑3542191

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