Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,523 A→G L239P (CTG→CCG)  yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,5230AG98.8% 513.2 / ‑2.8 161L239P (CTG→CCG) yieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (71/88);  total (72/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 168 total aligned reads displayed.

TACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGA  >  NC_000913/3937432‑3937612
                                                                                           |                                                                                         
tACGATACCGATCGAGAAAGGGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTTGCGTGTGTCGTGCATGGAGAATATAGAGAAGCAg                                                                                          >  2:2215589/1‑93 (MQ=255)
 aCGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGt                                                                                         <  1:338875/93‑1 (MQ=255)
 aCGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGt                                                                                         <  1:3741645/93‑1 (MQ=255)
  cGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTa                                                                                        >  2:1244602/1‑93 (MQ=255)
  cGATACAGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTa                                                                                        <  1:2865468/93‑1 (MQ=255)
   gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTaa                                                                                       >  1:1298312/1‑93 (MQ=255)
   gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTaa                                                                                       >  1:3292072/1‑93 (MQ=255)
   gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTaa                                                                                       >  1:953917/1‑93 (MQ=255)
    aTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAc                                                                                      <  2:3451872/93‑1 (MQ=255)
    aTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAc                                                                                      >  2:874220/1‑93 (MQ=255)
     tACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAccc                                                                                     >  1:1256990/1‑93 (MQ=255)
        cGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGa                                                                                  <  2:945665/93‑1 (MQ=255)
        cGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGa                                                                                  <  1:1281851/93‑1 (MQ=255)
        cGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGa                                                                                  >  1:2165039/1‑93 (MQ=255)
          aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGagc                                                                                >  1:1280367/1‑91 (MQ=255)
          aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATg                                                                                >  1:375643/1‑93 (MQ=255)
          aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATg                                                                                >  1:2018139/1‑93 (MQ=255)
          aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATg                                                                                >  1:657120/1‑93 (MQ=255)
           tCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               <  2:3611560/93‑1 (MQ=255)
           tCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               <  2:462737/93‑1 (MQ=255)
           tCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               <  2:2893037/93‑1 (MQ=255)
           tCGAGAAAGTGCGGGTTTTAAATAAATCTAAAGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               <  1:4083819/93‑1 (MQ=255)
           acGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               <  2:1818140/92‑1 (MQ=255)
            cGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCt                                                                              <  2:874813/93‑1 (MQ=255)
             gagaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTg                                                                             <  2:4380526/93‑1 (MQ=255)
             gagaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTg                                                                             >  1:4739631/1‑93 (MQ=255)
              agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTgg                                                                            >  2:436420/1‑93 (MQ=255)
              agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTgg                                                                            <  1:1691853/93‑1 (MQ=255)
                aaaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTg                                                                          <  1:1638935/93‑1 (MQ=255)
                 aaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGa                                                                         >  2:649479/1‑93 (MQ=255)
                 aaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGa                                                                         <  2:3933358/93‑1 (MQ=255)
                   gTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACa                                                                       >  1:2226896/1‑93 (MQ=255)
                   gTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACa                                                                       >  2:1825850/1‑93 (MQ=255)
                   gTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACa                                                                       >  1:4599862/1‑93 (MQ=255)
                      cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAAtt                                                                    <  1:3026163/93‑1 (MQ=255)
                      cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAAtt                                                                    <  2:729898/93‑1 (MQ=255)
                       gggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTa                                                                   >  2:3905874/1‑93 (MQ=255)
                       gggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTa                                                                   >  2:1527120/1‑93 (MQ=255)
                        ggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTAc                                                                  <  2:3304813/93‑1 (MQ=255)
                        ggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTAc                                                                  <  2:3146806/93‑1 (MQ=255)
                          ttttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCg                                                                <  2:4601548/93‑1 (MQ=255)
                           tttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGt                                                               <  1:4224592/93‑1 (MQ=255)
                               aataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                           <  2:2174683/93‑1 (MQ=255)
                               aataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                           <  1:2376362/93‑1 (MQ=255)
                                 taaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGc                                                         >  1:3414171/1‑93 (MQ=255)
                                 taaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGc                                                         <  1:1758560/93‑1 (MQ=255)
                                   aatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCaa                                                       >  1:1998247/1‑93 (MQ=255)
                                    atCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAAt                                                      >  1:499493/1‑93 (MQ=255)
                                    atCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAAt                                                      >  1:627339/1‑93 (MQ=255)
                                     tCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATc                                                     >  1:3219730/1‑93 (MQ=255)
                                      cTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATcc                                                    >  2:3779184/1‑93 (MQ=255)
                                      cTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATcc                                                    >  2:2708408/1‑93 (MQ=255)
                                       tAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCa                                                   >  2:1957729/1‑93 (MQ=255)
                                       tAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCa                                                   >  1:2154581/1‑93 (MQ=255)
                                          tGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCt                                                <  1:728393/93‑1 (MQ=255)
                                          tGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCt                                                <  2:2358183/93‑1 (MQ=255)
                                           gAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTg                                               >  2:1071115/1‑93 (MQ=255)
                                           gAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTg                                               <  1:4608311/93‑1 (MQ=255)
                                           gAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTg                                               <  1:3967239/93‑1 (MQ=255)
                                              aTTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCg                                            >  1:694027/1‑93 (MQ=255)
                                               ttAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGa                                           >  2:3202214/1‑93 (MQ=255)
                                               ttAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGCCAAGGCAATCCAGCTGGCGa                                           >  2:1762642/1‑93 (MQ=255)
                                                tAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGAc                                          >  2:1389182/1‑93 (MQ=255)
                                                 aaTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACa                                         >  2:4239826/1‑93 (MQ=255)
                                                  aTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACaa                                        >  1:138223/1‑93 (MQ=255)
                                                    gggTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATc                                      >  1:1094007/1‑93 (MQ=255)
                                                      gTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCtt                                    <  2:4656918/93‑1 (MQ=255)
                                                       tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  2:911329/93‑1 (MQ=255)
                                                       tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  1:1095730/93‑1 (MQ=255)
                                                       tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  2:2213743/93‑1 (MQ=255)
                                                        ttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttct                                  <  2:3292934/93‑3 (MQ=255)
                                                            cgtgcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTTACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcccc                              <  2:4148207/93‑1 (MQ=255)
                                                            cgtgcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcccc                              <  2:2038566/93‑1 (MQ=255)
                                                             gtgcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCg                             >  2:1722573/1‑93 (MQ=255)
                                                             gtgcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCg                             >  1:513242/1‑93 (MQ=255)
                                                               gcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGaa                           >  2:2647165/1‑93 (MQ=255)
                                                               gcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGaa                           <  1:1769580/93‑1 (MQ=255)
                                                               gcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGaa                           <  1:1937613/93‑1 (MQ=255)
                                                                ggtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAAt                          <  1:2870054/92‑1 (MQ=255)
                                                                cgtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAAt                          <  2:2621534/93‑1 (MQ=255)
                                                                 gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         <  1:270523/93‑1 (MQ=255)
                                                                 gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         >  2:1660560/1‑93 (MQ=255)
                                                                 gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         <  2:1426049/93‑1 (MQ=255)
                                                                      gCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctct                    <  1:3851105/93‑1 (MQ=255)
                                                                        aTGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAt                  >  1:1456053/1‑93 (MQ=255)
                                                                         tGGAGAATATAGAGAAGCGGTAACCGGATGCTAGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 <  2:2592694/93‑1 (MQ=255)
                                                                         tGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 <  2:4066974/93‑1 (MQ=255)
                                                                         tGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 >  1:4482710/1‑93 (MQ=255)
                                                                         tGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 >  2:784776/1‑93 (MQ=255)
                                                                         tGGAGAATATAGAGAAGCGGCAACCCGATGCTGGGGACAATTACCGTCAAGCCAATCCAGCTGGCGACACTCTTTTCACCGAATATCTCTAtt                 >  2:4475784/1‑93 (MQ=255)
                                                                          ggAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTc                <  2:2046535/93‑1 (MQ=255)
                                                                          ggAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTc                <  2:2411129/93‑1 (MQ=255)
                                                                            agaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCg              <  2:2430637/93‑1 (MQ=255)
                                                                            agaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCg              <  2:2980539/93‑1 (MQ=255)
                                                                             gaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGc             <  1:2974710/93‑1 (MQ=255)
                                                                              aaTATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCt            <  2:4096397/93‑1 (MQ=255)
                                                                              aaTATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCt            <  1:3111950/93‑1 (MQ=255)
                                                                                    gagaAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaa      >  2:1718460/1‑93 (MQ=255)
                                                                                      gaAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaaga    >  2:3555082/1‑93 (MQ=255)
                                                                                        aGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaagaga  >  1:2492324/1‑93 (MQ=255)
                                                                                           |                                                                                         
TACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGA  >  NC_000913/3937432‑3937612

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