Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,910,662 A→T D288E (GAT→GAA)  pstS ← phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,910,6620AT97.7% 303.6 / ‑1.3 86D288E (GAT→GAA) pstSphosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (43/41);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

GGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910565‑3910758
                                                                                                 |                                                                                                
ggCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTc                                                                                                >  2:1220960/1‑100 (MQ=255)
  cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  1:1644470/1‑100 (MQ=255)
     gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  1:222377/100‑1 (MQ=255)
     gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           >  2:1525439/1‑100 (MQ=255)
      cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          >  1:1495126/1‑100 (MQ=255)
       gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGCAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                         <  2:1733744/100‑1 (MQ=255)
       gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                         <  2:525111/100‑1 (MQ=255)
        tAATCCAGGTCGTTCGCCGGTTTTGTCCCCGTTTTGTACGCCCAGTCGAAGTATTTCAGCACTCCTGTGGCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          <  2:1633967/98‑1 (MQ=255)
          aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  2:1967273/100‑1 (MQ=255)
          aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  1:1654710/100‑1 (MQ=255)
           tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     <  1:18737/100‑1 (MQ=255)
           tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     <  1:464665/100‑1 (MQ=255)
           tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCCGTTCTTTGTGGATCa                                                                                     >  1:234314/1‑100 (MQ=255)
            ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    >  2:226918/1‑100 (MQ=255)
             cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                   >  1:1351240/1‑100 (MQ=255)
             cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                   >  1:214119/1‑100 (MQ=255)
              aGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCAAGTGGAAGTATTGCAGCTCTTCGGTGCCTTGGTCTGGTTTCTTCTGTTCTTTGTGGATCAGaa                                                                                  <  2:703231/100‑1 (MQ=255)
                gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATg                                                                                >  1:526319/1‑100 (MQ=255)
                   gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAAc                                                                             >  1:868393/1‑100 (MQ=255)
                      cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          <  1:226919/100‑1 (MQ=255)
                      cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          <  2:952620/100‑1 (MQ=255)
                       gCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTgg                                                                         >  1:762961/1‑100 (MQ=255)
                        ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        <  1:1220966/100‑1 (MQ=255)
                          tGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                      <  1:193320/100‑1 (MQ=255)
                          tGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                      >  1:1414285/1‑100 (MQ=255)
                           gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     >  2:576765/1‑100 (MQ=255)
                            gatctccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTagag                                                                    <  1:248942/95‑1 (MQ=255)
                              tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  1:576767/100‑1 (MQ=255)
                              tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  2:1495117/100‑1 (MQ=255)
                               aGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTa                                                                 <  1:686318/100‑1 (MQ=255)
                                gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTaa                                                                >  2:1914811/1‑100 (MQ=255)
                                gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTaa                                                                >  2:297125/1‑100 (MQ=255)
                                 ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTagagagatc                                                               >  2:248941/1‑95 (MQ=255)
                                 ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                               >  2:1661260/1‑100 (MQ=255)
                                  cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              <  1:1680784/100‑1 (MQ=255)
                                  cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  1:1227089/1‑100 (MQ=255)
                                      gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGcc                                                          >  1:837997/1‑100 (MQ=255)
                                        tttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGCTCAGAATGAACGTGGTAGCGGTACTAGGCCAt                                                        <  2:1706465/100‑1 (MQ=255)
                                          tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGc                                                      <  1:934422/100‑1 (MQ=255)
                                               gCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCtt                                                 <  1:591959/100‑1 (MQ=255)
                                               gCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCtt                                                 >  1:299123/1‑100 (MQ=255)
                                                cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTc                                                >  1:45642/1‑100 (MQ=255)
                                                cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTc                                                >  2:836083/1‑100 (MQ=255)
                                                    gTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCt                                            >  1:505379/1‑100 (MQ=255)
                                                     tCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAGTAGGCCATGCATCTTCGCCtt                                           >  1:86052/1‑100 (MQ=255)
                                                      cGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                          >  2:771414/1‑100 (MQ=255)
                                                      cGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                          <  2:517911/100‑1 (MQ=255)
                                                      cGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCtt                                           <  2:483494/99‑1 (MQ=255)
                                                       gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  2:609660/1‑100 (MQ=255)
                                                       gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         <  2:526318/100‑1 (MQ=255)
                                                       gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:587472/1‑100 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    >  2:2017742/1‑100 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  2:219262/100‑1 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  1:1666616/100‑1 (MQ=255)
                                                             tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                   <  1:1039614/100‑1 (MQ=255)
                                                             tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  2:937061/99‑1 (MQ=255)
                                                              ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  >  2:1180510/1‑100 (MQ=255)
                                                                cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTATGTGGATCAGAATGAACGTGGAAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                >  1:721498/1‑100 (MQ=255)
                                                                 aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCa                               >  2:1513225/1‑100 (MQ=255)
                                                                   cACTTCTGTGCCTTGTTCTGGTTTCTTCCGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                             <  2:234313/100‑1 (MQ=255)
                                                                      ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATcc                          <  1:1419897/100‑1 (MQ=255)
                                                                      ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATcc                          <  1:1407710/100‑1 (MQ=255)
                                                                       tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGCGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                         >  1:925119/1‑100 (MQ=255)
                                                                        cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        <  1:836086/100‑1 (MQ=255)
                                                                        cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        <  1:283706/100‑1 (MQ=255)
                                                                         tgtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGa                       <  2:1644457/100‑1 (MQ=255)
                                                                          gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:570152/1‑100 (MQ=255)
                                                                          gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:1633215/1‑100 (MQ=255)
                                                                          gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:1643500/1‑100 (MQ=255)
                                                                            gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                     <  2:914349/99‑1 (MQ=255)
                                                                            gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    >  1:1082710/1‑100 (MQ=255)
                                                                            gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTAGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    >  1:1181590/1‑100 (MQ=255)
                                                                               ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAg                 >  1:1656697/1‑100 (MQ=255)
                                                                               ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAg                 <  1:869412/100‑1 (MQ=255)
                                                                                tGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAgg                <  2:1181584/100‑1 (MQ=255)
                                                                                  ttCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAACAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtt              >  2:1128022/1‑100 (MQ=255)
                                                                                      ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          >  1:418165/1‑100 (MQ=255)
                                                                                      ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          >  1:938666/1‑100 (MQ=255)
                                                                                      ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:1902261/100‑1 (MQ=255)
                                                                                      ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGAGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  2:1082705/100‑1 (MQ=255)
                                                                                         ttcttcTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcc       <  2:1351231/100‑1 (MQ=255)
                                                                                           cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg     >  2:48334/1‑100 (MQ=255)
                                                                                           cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg     <  1:1794057/100‑1 (MQ=255)
                                                                                           cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg     >  2:1959477/1‑100 (MQ=255)
                                                                                             nctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc   >  2:1631790/2‑100 (MQ=255)
                                                                                              ctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  2:629148/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910565‑3910758

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