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,6620AT96.3% 270.5 / ‑0.8 81D288E (GAT→GAApstSphosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (37/41);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910565‑3910758
                                                                                                 |                                                                                                
ggCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTc                                                                                                <  1:10171/100‑1 (MQ=255)
  cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCGCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTATGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  2:1410716/1‑100 (MQ=255)
  cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  2:438420/1‑100 (MQ=255)
    ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtg                                                                                            <  2:2193740/100‑1 (MQ=255)
    ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtg                                                                                            >  2:752933/1‑100 (MQ=255)
     gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATccacaa                                                                                           >  2:1072937/1‑95 (MQ=255)
      cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          >  2:1476684/1‑100 (MQ=255)
       gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                         <  2:2212257/100‑1 (MQ=255)
        tAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                        >  2:421695/1‑100 (MQ=255)
         aaTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGAt                                                                                       <  2:1771629/100‑1 (MQ=255)
         aaTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGAt                                                                                       <  2:841717/100‑1 (MQ=255)
          aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  2:1880318/100‑1 (MQ=255)
           tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     <  2:1501692/100‑1 (MQ=255)
           tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     >  1:69692/1‑100 (MQ=255)
            ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      >  2:465666/1‑98 (MQ=255)
            ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    >  2:483078/1‑100 (MQ=255)
              aGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGaa                                                                                  <  2:794074/100‑1 (MQ=255)
               ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  2:2313039/100‑1 (MQ=255)
               ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  1:1334579/100‑1 (MQ=255)
               ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  2:1634075/100‑1 (MQ=255)
                gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  1:1137479/99‑1 (MQ=255)
                  cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGaa                                                                              >  2:1413867/1‑100 (MQ=255)
                   gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAAc                                                                             >  2:1239869/1‑100 (MQ=255)
                    ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  1:702542/100‑1 (MQ=255)
                      cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          >  2:1497551/1‑100 (MQ=255)
                        ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        >  2:1204388/1‑100 (MQ=255)
                          tGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                      >  1:723128/1‑100 (MQ=255)
                           gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     >  1:1691816/1‑100 (MQ=255)
                           gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     >  2:35565/1‑100 (MQ=255)
                                gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTATTTGTGGATCAGAATGAACGTGGTAGAGGTaa                                                                >  2:947700/1‑100 (MQ=255)
                                 ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCGCTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                               >  2:1122758/1‑100 (MQ=255)
                                  cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  2:1616749/1‑100 (MQ=255)
                                  cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  1:704624/1‑100 (MQ=255)
                                     ggTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                           >  2:1467420/1‑100 (MQ=255)
                                      gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGcc                                                          <  1:1498123/100‑1 (MQ=255)
                                        tttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCAt                                                        <  1:35566/100‑1 (MQ=255)
                                           gTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                     >  1:1614203/1‑100 (MQ=255)
                                           gTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                     <  1:1239872/100‑1 (MQ=255)
                                           gTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                     >  2:1284407/1‑100 (MQ=255)
                                            tACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCAt                                                    >  1:1530628/1‑100 (MQ=255)
                                            tACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCAt                                                    >  1:478627/1‑100 (MQ=255)
                                                 ccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCg                                               <  2:1775515/100‑1 (MQ=255)
                                                    gTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCt                                            <  1:2053165/100‑1 (MQ=255)
                                                       gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:886463/1‑100 (MQ=255)
                                                         agaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtg                                       <  2:2020892/100‑1 (MQ=255)
                                                          gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtgg                                      <  2:1194403/100‑1 (MQ=255)
                                                          gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtgg                                      <  1:360016/100‑1 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  1:437198/100‑1 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  2:478625/100‑1 (MQ=255)
                                                            aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    >  2:2091947/1‑100 (MQ=255)
                                                            aTTTCAGCACCTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGCGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  2:1530624/100‑1 (MQ=255)
                                                             tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                   <  2:2226844/100‑1 (MQ=255)
                                                             tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                   <  1:870760/100‑1 (MQ=255)
                                                              ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  1:1467424/100‑1 (MQ=255)
                                                              ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  2:723127/100‑1 (MQ=255)
                                                                cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                <  1:162886/100‑1 (MQ=255)
                                                                cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                <  2:1543315/100‑1 (MQ=255)
                                                                   cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTCGATGTAATATGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                             >  2:1640892/1‑100 (MQ=255)
                                                                   cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                             <  1:1476688/100‑1 (MQ=255)
                                                                   cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                             >  1:716585/1‑100 (MQ=255)
                                                                     cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATc                           <  2:1690040/100‑1 (MQ=255)
                                                                       tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                         >  1:2021544/1‑100 (MQ=255)
                                                                       tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                         <  2:1319317/100‑1 (MQ=255)
                                                                        cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        <  2:2021537/100‑1 (MQ=255)
                                                                        cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        <  2:2258965/100‑1 (MQ=255)
                                                                        cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        <  2:886461/100‑1 (MQ=255)
                                                                        cTGTGCCTTGTTCTGAATTCTTCAGTTGTAAGTGGATCAGCATGCACGTGGTAGATGTAATATGCCATGCATCTTCGCCTTTCTGTTTGGTCAGATCCTg                        >  2:479748/1‑100 (MQ=255)
                                                                          gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:2148101/1‑100 (MQ=255)
                                                                          gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:113497/1‑100 (MQ=255)
                                                                           tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCCTGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                     <  1:1647880/100‑1 (MQ=255)
                                                                                tGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAgg                <  1:408579/100‑1 (MQ=255)
                                                                                 gTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGt               >  1:320028/1‑100 (MQ=255)
                                                                                  ttCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtt              <  2:1707171/100‑1 (MQ=255)
                                                                                   tCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGttt             <  2:704623/100‑1 (MQ=255)
                                                                                    cTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtttt            >  2:1812539/1‑100 (MQ=255)
                                                                                         ttcttcTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcc       >  1:1032600/1‑100 (MQ=255)
                                                                                           cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg     >  1:1176521/1‑100 (MQ=255)
                                                                                           cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg     >  1:477445/1‑100 (MQ=255)
                                                                                             tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGt    >  2:1901045/1‑99 (MQ=255)
                                                                                             tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc   <  1:465668/100‑1 (MQ=255)
                                                                                              ctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  2:163728/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: