Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,910,662 A→T 100% 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,6620AT94.5% 293.7 / 5.3 92D288E (GAT→GAApstSphosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (38/48);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGT  >  NC_000913/3910567‑3910756
                                                                                               |                                                                                              
cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCtt                                                                                            >  2:449574/1‑100 (MQ=255)
 nggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCttt                                                                                           >  2:165291/2‑100 (MQ=255)
  ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtg                                                                                          >  2:1804462/1‑100 (MQ=255)
   gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                         <  2:143578/100‑1 (MQ=255)
   gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                         >  1:1465907/1‑100 (MQ=255)
   gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                         >  1:2197286/1‑100 (MQ=255)
   gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                         >  1:167620/1‑100 (MQ=255)
     gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                       >  1:1809404/1‑100 (MQ=255)
      tAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                      >  2:1250802/1‑100 (MQ=255)
      tAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                      <  2:1142964/100‑1 (MQ=255)
       aaTCCAGGTCGGTCGCCTGTTTAGCCCCGGTTTTGTACGGCCAGTCGGAGAATTTCAGCACTTCTGTGCCTTGTTTTGGTCTCTTTTGTTCTTtttgta                                                                                      >  2:1174974/1‑97 (MQ=255)
          ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGCGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                  <  1:1804476/100‑1 (MQ=255)
          ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCGGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                   <  1:1080589/99‑1 (MQ=255)
           cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                 >  1:117611/1‑100 (MQ=255)
            aGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGaa                                                                                <  1:1962568/100‑1 (MQ=255)
             ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                               <  1:1851476/100‑1 (MQ=255)
                cggtCGCCTGTTTAGCCCCGGTTTTTTACGCCCAGTCGTAGAGTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                            <  2:794880/97‑1 (MQ=255)
                 gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAAc                                                                           >  1:2076447/1‑100 (MQ=255)
                 gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAAc                                                                           >  2:2509884/1‑100 (MQ=255)
                 gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAAc                                                                           >  1:1665390/1‑100 (MQ=255)
                  ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                          >  1:2139683/1‑100 (MQ=255)
                  ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                          <  2:678992/100‑1 (MQ=255)
                   tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGt                                                                         >  1:380701/1‑100 (MQ=255)
                     gCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTgg                                                                       >  1:1085722/1‑100 (MQ=255)
                      ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                      <  2:2310256/100‑1 (MQ=255)
                      ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                      <  2:892431/100‑1 (MQ=255)
                      ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGt                                                                      <  2:1292263/100‑1 (MQ=255)
                       cTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTa                                                                     <  1:1893340/100‑1 (MQ=255)
                       cTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTa                                                                     <  1:1893511/100‑1 (MQ=255)
                        tGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                    >  2:354532/1‑100 (MQ=255)
                        tGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTag                                                                    <  2:2019465/100‑1 (MQ=255)
                          tgtAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTagag                                                                  <  1:189168/98‑1 (MQ=255)
                           ttAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                 <  1:354536/100‑1 (MQ=255)
                            tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                <  2:1741093/100‑1 (MQ=255)
                            tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATGTCAGCACTTCTGTGCCTTGTTCTGGGTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                <  1:2481149/100‑1 (MQ=255)
                             aGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTa                                                               <  1:2283408/100‑1 (MQ=255)
                             aGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTa                                                               <  1:465367/100‑1 (MQ=255)
                              gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGGGCCTCGTACTGGCTTCTTCTGCTCTTTGTGGATCAGAATGAACGTGGTAGTGGt                                                                >  2:549129/1‑98 (MQ=255)
                               ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                             >  1:1026139/1‑100 (MQ=255)
                               ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                             >  1:1468183/1‑100 (MQ=255)
                               ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGGTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                 >  2:2375974/1‑96 (MQ=255)
                                cccGGTTTTGTACGCCCAGTCGAAGAATTTCATCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                            <  2:2378189/100‑1 (MQ=255)
                                 ccGGTTTTGTACGCCCAGTCGAAGAATCTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAg                                                           <  1:2086729/100‑1 (MQ=255)
                                  ccgTTTTGGACGCGCAGTCGTAGAATTTCAGCACTTCTGCGCCTTGTACTGGTTTCTTCTGTTCTTTGTGGAGCAGAATGAACGTGGTAGAGGTa                                                               <  2:1358098/93‑1 (MQ=255)
                                  cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                          <  2:1825479/100‑1 (MQ=255)
                                   ggTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                         >  1:741797/1‑100 (MQ=255)
                                   ggTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                         <  2:1924902/100‑1 (MQ=255)
                                    gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGcc                                                        >  2:2120056/1‑100 (MQ=255)
                                    gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGc                                                         <  1:222325/99‑1 (MQ=255)
                                     ttttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCa                                                       <  2:1860382/100‑1 (MQ=255)
                                       ttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATg                                                     <  1:2328458/100‑1 (MQ=255)
                                       ttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATg                                                     <  2:350580/100‑1 (MQ=255)
                                       ttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATg                                                     <  2:1360519/100‑1 (MQ=255)
                                        tGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGc                                                    >  1:617348/1‑100 (MQ=255)
                                          tatgctcattCGAAGAATTTCATCATTTCTGTGCTTTGTTCTGGTTTCTTTTGTTTTTTGTGGATCAGAATGAATGTGGTAGAGGTAATAGGCCATGCa                                                   <  2:663255/90‑1 (MQ=255)
                                          tACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCAt                                                  >  1:955561/1‑100 (MQ=255)
                                             gCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCtt                                               >  1:2372128/1‑100 (MQ=255)
                                              cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTCTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTc                                              >  2:1379677/1‑100 (MQ=255)
                                                cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGc                                            >  2:309904/1‑100 (MQ=255)
                                                  gTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCt                                          <  1:1595556/100‑1 (MQ=255)
                                                     gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTCGAGGTAATAGGCCATGCATCTTCGCCTTTc                                       >  2:2295151/1‑100 (MQ=255)
                                                      aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                      <  2:895502/100‑1 (MQ=255)
                                                       agaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtg                                     <  2:1475219/100‑1 (MQ=255)
                                                        gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtgg                                    >  2:2077254/1‑100 (MQ=255)
                                                        gaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtgg                                    <  1:449578/100‑1 (MQ=255)
                                                         aaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                   <  1:439588/100‑1 (MQ=255)
                                                          aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                  >  2:475122/1‑100 (MQ=255)
                                                          aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                  <  1:567825/100‑1 (MQ=255)
                                                            ttCAGCCCTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGTATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgc                                <  2:1144790/100‑2 (MQ=255)
                                                            ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                >  2:31600/1‑100 (MQ=255)
                                                               aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCa                             <  1:2195544/100‑1 (MQ=255)
                                                                gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAg                            >  1:613745/1‑100 (MQ=255)
                                                                 cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                           >  1:1191292/1‑100 (MQ=255)
                                                                  aCTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGAt                          <  2:751778/100‑1 (MQ=255)
                                                                     tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                       >  2:159769/1‑100 (MQ=255)
                                                                        gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                    <  1:1863582/100‑1 (MQ=255)
                                                                        gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGACTTTCTGGTTGGTCAGATCCTGAg                    >  2:1795312/1‑100 (MQ=255)
                                                                         tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                   >  1:1769487/1‑100 (MQ=255)
                                                                            cTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGaa                <  2:1809390/100‑1 (MQ=255)
                                                                            cTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGaa                >  1:2158440/1‑100 (MQ=255)
                                                                              tGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAgg              <  1:1283400/100‑1 (MQ=255)
                                                                              tGTTCTGGTTTCTTCTGTTCTTTGTGAATCAGAACGACCGTGGTAGAGGTAATAGGCCATGCAACTTCGCCTTTCTGGTTGATCAGATCCTGAGCGAAgg              <  2:2430971/100‑1 (MQ=255)
                                                                                ttCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtt            >  1:437007/1‑100 (MQ=255)
                                                                                 tCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGttt           <  1:874183/100‑1 (MQ=255)
                                                                                   tGGTTTCTTCTGTTCTTTGTGGAGCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtttt          <  2:2099560/99‑1 (MQ=255)
                                                                                    ngTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc        >  2:2253664/2‑100 (MQ=255)
                                                                                    ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc        >  1:575804/1‑100 (MQ=255)
                                                                                      tttCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTc      <  1:2253678/100‑1 (MQ=255)
                                                                                         cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg   <  1:920255/100‑1 (MQ=255)
                                                                                         cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg   >  1:514660/1‑100 (MQ=255)
                                                                                          ttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcc     <  2:1552689/97‑1 (MQ=255)
                                                                                          ttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGt  <  2:1744897/100‑1 (MQ=255)
                                                                                               |                                                                                              
CTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGT  >  NC_000913/3910567‑3910756

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: