Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 15.8% intergenic (+43/‑187) tufB → / → secE translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA15.8% 210.0 / 30.1 88intergenic (+43/‑187)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (30/44);  new base A (8/6);  total (38/50)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

ATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTAT  >  NC_000913/4177072‑4177268
                                                                                                   |                                                                                                 
aTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGcatcat                                                                                                   <  1:203578/100‑1 (MQ=255)
 tCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAtt                                                                                                  <  2:2215972/100‑1 (MQ=255)
   cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                >  2:1120134/1‑100 (MQ=255)
   cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                >  1:334162/1‑100 (MQ=255)
    gTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                               <  2:2542563/100‑1 (MQ=255)
    gTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                               <  1:1153751/100‑1 (MQ=255)
     tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                              <  1:116109/100‑1 (MQ=255)
        aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                           <  2:1897069/100‑1 (MQ=255)
        aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                           >  1:1894194/1‑100 (MQ=255)
        aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                           >  2:2432040/1‑100 (MQ=255)
         ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAACAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                          >  2:1531160/1‑100 (MQ=255)
           cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                        <  1:680691/100‑1 (MQ=255)
            ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                       <  1:489113/100‑1 (MQ=255)
             gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTTTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                      <  2:2022555/100‑1 (MQ=255)
             gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                      <  1:576619/100‑1 (MQ=255)
               cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                    >  2:2501989/1‑100 (MQ=255)
               cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                    <  1:881228/100‑1 (MQ=255)
               cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                    <  2:245108/100‑3 (MQ=255)
                 tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                  <  1:1703122/100‑1 (MQ=255)
                 tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                  <  1:2072972/100‑1 (MQ=255)
                 tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                  <  1:222884/100‑5 (MQ=38)
                       tGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCttt                                                                            <  2:483870/100‑1 (MQ=255)
                        ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                           >  2:540666/1‑100 (MQ=255)
                         gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                          <  2:2329790/100‑1 (MQ=255)
                          cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                         <  2:2495220/100‑1 (MQ=255)
                             gggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  2:1560869/100‑1 (MQ=255)
                             gggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  2:398479/100‑1 (MQ=255)
                             gggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  1:1282850/100‑1 (MQ=255)
                             gggCGTTGTAGCAAAAGTTCTGAGCGAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  1:1753805/100‑1 (MQ=255)
                              ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                     >  1:1876706/1‑100 (MQ=255)
                                  ttGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAc                                                                 <  2:843266/100‑1 (MQ=255)
                                   tGTAGCAAAAGTTCTGGGCTAATTGGCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAc                                                                 <  2:884486/99‑1 (MQ=255)
                                   tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                <  1:1715180/100‑1 (MQ=255)
                                     tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                              >  1:1627447/1‑100 (MQ=255)
                                      aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                             <  1:2437075/100‑1 (MQ=255)
                                      aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                             >  2:722350/1‑100 (MQ=255)
                                      aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                             >  2:1119746/1‑100 (MQ=255)
                                       gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            >  2:3936/1‑100 (MQ=255)
                                       gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            <  1:60633/100‑1 (MQ=255)
                                        cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                           >  1:2380850/1‑100 (MQ=255)
                                        cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                           <  2:1939563/100‑1 (MQ=255)
                                        cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                           <  2:337166/100‑1 (MQ=255)
                                              ttCTGAGCTAATTGCCGATAACATTGGACGCAATGCGCACTAAAAGGTCATCATTTGATGCTCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                     <  1:608848/100‑1 (MQ=255)
                                               tCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATAATTTGTTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                      >  2:151009/1‑98 (MQ=255)
                                                cTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                   >  2:80627/1‑100 (MQ=255)
                                                 tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                  <  1:870634/100‑1 (MQ=255)
                                                  gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                 <  1:99528/100‑1 (MQ=255)
                                                   aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                 >  1:1426501/1‑99 (MQ=255)
                                                    gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGCttt                                               >  2:1423380/1‑100 (MQ=255)
                                                       aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                            >  2:980424/1‑100 (MQ=255)
                                                       aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                            <  1:1590260/100‑1 (MQ=255)
                                                        aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                           <  1:503162/100‑1 (MQ=255)
                                                          tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                          <  2:1453939/99‑1 (MQ=255)
                                                          tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                         <  2:1824670/100‑1 (MQ=255)
                                                           gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                        >  1:2306530/1‑100 (MQ=255)
                                                             ccaTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                      >  1:1733683/3‑100 (MQ=255)
                                                             cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                      <  1:2502004/100‑1 (MQ=255)
                                                              gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                     <  1:1000411/100‑1 (MQ=255)
                                                                 aaCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAAt                                  <  1:1119754/100‑1 (MQ=255)
                                                                   cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                >  2:305870/1‑100 (MQ=255)
                                                                    aTTTGACGTAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                               >  2:414989/1‑100 (MQ=255)
                                                                      ttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTg                             >  1:1831154/1‑100 (MQ=255)
                                                                        gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                           <  1:333116/100‑1 (MQ=255)
                                                                         aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                          <  1:1650925/100‑1 (MQ=255)
                                                                         aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                          >  2:2516468/1‑100 (MQ=255)
                                                                          cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                         >  2:2008691/1‑100 (MQ=255)
                                                                          cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                         <  1:129508/100‑1 (MQ=255)
                                                                          cGCAATGCGCACCAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                          <  1:1126330/99‑1 (MQ=255)
                                                                                gcgtgcaaaaaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGACTTTCTTTGTCATAATCATTGCTGAGACAgg                   >  2:2551141/8‑100 (MQ=21)
                                                                                 cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                  <  1:850094/100‑1 (MQ=255)
                                                                                 cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                  <  1:602506/100‑1 (MQ=255)
                                                                                 cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                  <  2:1050256/100‑1 (MQ=255)
                                                                                 cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                  <  1:1292554/100‑1 (MQ=255)
                                                                                  gcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                 >  2:1845824/1‑100 (MQ=255)
                                                                                   tgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                >  1:1053268/5‑100 (MQ=38)
                                                                                   tgcaaaaaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                >  2:868895/5‑100 (MQ=25)
                                                                                   cACTAGAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCCGTGATTTTCTTTCTCATAATCATTGCTGAGACAGGct                 >  2:2237067/1‑99 (MQ=255)
                                                                                   cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                >  2:438277/1‑100 (MQ=255)
                                                                                    aCTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct               <  1:1570506/100‑1 (MQ=255)
                                                                                     caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg              >  1:417327/3‑100 (MQ=255)
                                                                                      aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATTAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt             >  2:1803244/2‑100 (MQ=39)
                                                                                      aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt             >  1:388390/2‑100 (MQ=255)
                                                                                      aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCCCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt             >  2:1267519/2‑100 (MQ=39)
                                                                                        aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg           <  1:1211146/100‑1 (MQ=255)
                                                                                         aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGCAACTCCATCATTAGTGATTTTATATGTCCAAATCCTTGCTGAGCCAGCCACTGt             >  2:1135275/1‑97 (MQ=255)
                                                                                          aGGGCATCATTTGATGCCCTTTTTTTTTTTTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTCAg         >  2:382467/1‑100 (MQ=255)
                                                                                          aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTCTCATAATCATTGCTGAGACAGGCTCTGTTGAg         <  1:2432054/100‑1 (MQ=255)
                                                                                           gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg        <  1:1962267/100‑1 (MQ=255)
                                                                                           gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg        >  1:574304/1‑100 (MQ=255)
                                                                                           gggCATCATTTGATGCCCTGTGTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTGCTTTGTCATGATCATTGCTGAGACAGGCTCTGTTGAgg        <  2:219712/100‑1 (MQ=255)
                                                                                              catcatTTGATGTCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg     <  1:405477/100‑1 (MQ=255)
                                                                                               atcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt    >  1:1097358/1‑100 (MQ=255)
                                                                                               atcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt    <  2:1067207/100‑1 (MQ=255)
                                                                                                tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta   <  2:1638055/100‑1 (MQ=255)
                                                                                                 catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat  >  1:477584/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTAT  >  NC_000913/4177072‑4177268

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: