Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,422,395 A→G S45S (AGT→AGC)  hisP ← ATP‑binding component of histidine transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,422,3950AG87.2% 570.2 / 64.0 219S45S (AGT→AGC) hisPATP‑binding component of histidine transport
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (81/110);  total (94/125)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 226 total aligned reads displayed.

TCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTA  >  NC_000913/2422303‑2422487
                                                                                            |                                                                                            
tCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTa                                                                                              >  2:3437953/1‑93 (MQ=255)
 cTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTAc                                                                                             <  1:614759/93‑1 (MQ=255)
  tttGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCt                                                                                            <  1:2804300/93‑1 (MQ=255)
   ttGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACtt                                                                                           <  1:5471846/93‑1 (MQ=255)
    tGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCttt                                                                                          <  2:814836/93‑1 (MQ=255)
    tGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCttt                                                                                          <  1:3134243/93‑1 (MQ=255)
    tGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACttt                                                                                          <  2:5395936/93‑1 (MQ=255)
     gTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCtttt                                                                                         >  2:4419086/1‑93 (MQ=255)
     gTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACtttt                                                                                         >  1:177134/1‑93 (MQ=255)
      tCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTc                                                                                        <  2:2520375/93‑1 (MQ=255)
        gcgcACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTccc                                                                                      >  1:3310146/1‑93 (MQ=255)
         cgcACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTcccc                                                                                     <  1:88045/93‑1 (MQ=255)
          gcACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCg                                                                                    >  2:6599564/1‑93 (MQ=255)
            aCCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGAt                                                                                  >  1:343533/1‑93 (MQ=255)
            aCCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGAt                                                                                  >  2:83751/1‑93 (MQ=255)
             ccAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATc                                                                                 >  2:6254959/1‑93 (MQ=255)
              cAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATcc                                                                                <  1:823182/93‑1 (MQ=255)
               aGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATccc                                                                               >  2:1022321/1‑93 (MQ=255)
               aGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATccc                                                                               <  2:2037063/93‑1 (MQ=255)
                gATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCcg                                                                              <  1:4255908/93‑1 (MQ=255)
                gATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTATTCCCCGATCCcg                                                                              >  1:4226404/1‑93 (MQ=255)
                 aTTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCcga                                                                             <  1:4758255/93‑1 (MQ=255)
                  ttGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCcgac                                                                            <  2:3983964/93‑1 (MQ=255)
                    gATCGTCTTGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCcgacga                                                                          <  2:428486/93‑1 (MQ=255)
                    gATCGTCTGGCCATTGACTACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGAGCTTTTCCCCGATCCcgacga                                                                          >  2:1753079/1‑93 (MQ=255)
                     aTCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGAt                                                                         <  1:5556440/93‑1 (MQ=255)
                      tCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATc                                                                        >  1:2488484/1‑93 (MQ=255)
                      tCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATc                                                                        <  1:4715153/93‑1 (MQ=255)
                       cGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATcc                                                                       <  1:1574706/93‑1 (MQ=255)
                         tCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCga                                                                     >  2:1647289/1‑93 (MQ=255)
                           tGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCgatg                                                                   >  2:3223562/1‑93 (MQ=255)
                           tGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCgatg                                                                   <  2:292254/93‑1 (MQ=255)
                            ggCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCgatga                                                                  <  2:5107479/93‑1 (MQ=255)
                            ggCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCgatga                                                                  >  1:1725650/1‑93 (MQ=255)
                            ggCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCgatga                                                                  >  1:4018743/1‑93 (MQ=255)
                             gCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCgatgat                                                                 <  2:351972/93‑1 (MQ=255)
                               cATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGc                                                               >  1:2713390/1‑93 (MQ=255)
                                aTTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCt                                                              <  2:3625482/93‑1 (MQ=255)
                                aTTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCt                                                              <  1:2210381/93‑1 (MQ=255)
                                aTTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCt                                                              >  2:4355123/1‑93 (MQ=255)
                                 ttGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCtt                                                             <  1:3884520/93‑1 (MQ=255)
                                  tGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCtta                                                            <  1:2988211/93‑1 (MQ=255)
                                   gACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCttat                                                           >  2:2790321/1‑93 (MQ=255)
                                    aCCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCttatt                                                          >  1:2952153/1‑93 (MQ=255)
                                     ccACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCttatta                                                         <  1:2736016/93‑1 (MQ=255)
                                       aCGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACa                                                       >  2:2176968/1‑93 (MQ=255)
                                       aCGATCCACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACa                                                       >  1:4559046/1‑93 (MQ=255)
                                         gATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACAtc                                                     >  2:2704237/1‑93 (MQ=255)
                                         gATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACAtc                                                     <  1:1250643/93‑1 (MQ=255)
                                          aTCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACAtct                                                    <  1:1149619/93‑1 (MQ=255)
                                          aTCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACAtct                                                    <  1:6306114/93‑1 (MQ=255)
                                             gACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCg                                                 >  2:4082639/1‑93 (MQ=255)
                                              aCCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCgg                                                >  1:920927/1‑93 (MQ=255)
                                               cccTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGc                                               >  1:714811/1‑93 (MQ=255)
                                               cccTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGGGCTTTTCTCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGc                                               >  1:2340146/1‑93 (MQ=255)
                                                 cTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATCACATCTCCGGCAt                                             <  1:3179450/93‑1 (MQ=255)
                                                    cACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCg                                          <  1:1928595/93‑1 (MQ=255)
                                                     aCTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGc                                         <  2:755683/93‑1 (MQ=255)
                                                     aCTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGc                                         >  1:1060993/1‑93 (MQ=255)
                                                       tCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                       >  1:750023/1‑93 (MQ=255)
                                                       tCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                       >  2:640315/1‑93 (MQ=255)
                                                       tCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                       <  1:5741830/93‑1 (MQ=255)
                                                       tCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                       >  2:4253431/1‑93 (MQ=255)
                                                       tCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                       <  2:4225116/93‑1 (MQ=255)
                                                          gTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCa                                    >  1:3700514/1‑93 (MQ=255)
                                                           tttttCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAg                                   >  1:1118720/1‑93 (MQ=255)
                                                           tttttCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAg                                   <  1:2593343/93‑1 (MQ=255)
                                                           tttttCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAg                                   <  1:99184/93‑1 (MQ=255)
                                                            ttttCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGt                                  >  2:1947521/1‑93 (MQ=255)
                                                              ttCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGa                                >  1:1410135/1‑93 (MQ=255)
                                                              ttCGAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGAGCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGa                                <  2:273822/93‑1 (MQ=255)
                                                                 gAGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATAc                             >  2:1922472/1‑93 (MQ=255)
                                                                  aGGAAGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATAcc                            >  1:1000664/1‑93 (MQ=255)
                                                                  aGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATAcc                            >  1:355372/1‑93 (MQ=255)
                                                                     aaGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATTATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCt                         <  2:5312064/93‑1 (MQ=255)
                                                                     aaGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCt                         >  2:5652465/1‑93 (MQ=255)
                                                                      aGTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCtt                        <  1:4577133/93‑1 (MQ=255)
                                                                       gTTAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCttt                       >  1:5975909/1‑93 (MQ=255)
                                                                        ttAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTc                      <  1:6226056/93‑1 (MQ=255)
                                                                        ttAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTc                      <  1:5161174/93‑1 (MQ=255)
                                                                        ttAATGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTc                      <  1:266714/93‑1 (MQ=255)
                                                                          aaTGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAg                    >  1:5144974/1‑93 (MQ=255)
                                                                          aaTGCAGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAg                    >  1:5498285/1‑93 (MQ=255)
                                                                          aaTGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAg                    <  1:3651681/93‑1 (MQ=255)
                                                                            aggaGCGCAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCa                  <  2:1382532/90‑1 (MQ=255)
                                                                                  gaaGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCAt            <  1:1363058/91‑1 (MQ=255)
                                                                                   cAGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATg           <  1:644554/93‑1 (MQ=255)
                                                                                    aGAAAGTTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGt          <  1:5629346/93‑1 (MQ=255)
                                                                                    aGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGt          >  2:4795684/1‑93 (MQ=255)
                                                                                    aGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGt          <  1:3756770/93‑1 (MQ=255)
                                                                                    aGAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGt          <  1:5215989/93‑1 (MQ=255)
                                                                                     gAAAGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGtt         >  2:1232205/1‑93 (MQ=255)
                                                                                      aaaGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTc        <  1:5707268/93‑1 (MQ=255)
                                                                                      aaaGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTc        >  2:5144799/1‑93 (MQ=255)
                                                                                      aaaGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTc        >  2:3651557/1‑93 (MQ=255)
                                                                                       aaGGTGCTTTTCGCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCg       >  1:5310027/1‑93 (MQ=255)
                                                                                       aaGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCg       >  1:1825798/1‑93 (MQ=255)
                                                                                        aGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGc      <  1:5167036/93‑1 (MQ=255)
                                                                                        aGGTGCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGc      <  2:3131611/93‑1 (MQ=255)
                                                                                            aCTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTa  <  1:4363531/93‑1 (MQ=255)
                                                                                            |                                                                                            
TCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTA  >  NC_000913/2422303‑2422487

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