Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,646 +T coding (5/1500 nt) pitA → metal phosphate:H(+) symporter PitA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,637,6451.T97.9% 689.7 / 6.2 189coding (4/1500 nt)pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base T (82/103);  total (84/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 196 total aligned reads displayed.

GGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCA  >  NC_000913/3637547‑3637739
                                                                                                    |                                                                                              
ggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCt                                                                                               >  1:2995585/1‑100 (MQ=255)
  cAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCtta                                                                                             >  1:2205912/1‑98 (MQ=255)
   aaaTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCttac                                                                                            >  1:754103/1‑97 (MQ=255)
     aTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCttacat                                                                                          >  1:2821443/1‑95 (MQ=255)
      tCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtt                                                                                          <  1:1100248/100‑1 (MQ=255)
       cAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAttt                                                                                         <  1:1256931/100‑1 (MQ=255)
        aaaaaaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttg                                                                                        <  1:2079877/100‑1 (MQ=255)
            aaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCtt                                                                                              <  2:2512307/89‑2 (MQ=255)
             aaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                   <  1:1219543/100‑1 (MQ=255)
              aGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCt                                                                                  >  1:580292/1‑100 (MQ=255)
                  tataTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGcc                                                                              <  1:3410769/100‑1 (MQ=255)
                   atatTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGCTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttg                                                                                        <  2:1196438/89‑1 (MQ=39)
                    tatTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTg                                                                            <  1:1617521/100‑1 (MQ=255)
                     atTTCACTTTGCCCGCGCCGCGAAAGTCGCTGATAATGCGCCGAATGGATGTCTTCAAAATGGCGTAACGTCCTATGATTACAtttgtt                                                                                      <  2:1261399/89‑1 (MQ=25)
                     atTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtt                                                                                      <  2:3760158/89‑1 (MQ=255)
                     atTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTgg                                                                           >  1:900218/1‑100 (MQ=255)
                      tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGa                                                                          >  1:394695/1‑100 (MQ=255)
                       ttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtttg                                                                                    >  2:1502499/1‑89 (MQ=255)
                        tCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAtt                                                                        <  1:2256812/100‑1 (MQ=255)
                         cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCt                                                                                  <  2:754103/89‑1 (MQ=255)
                         cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCt                                                                                  >  2:341766/1‑89 (MQ=255)
                          aCTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTg                                                                      <  1:2550959/100‑1 (MQ=255)
                            tttGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGc                                                                               <  2:180106/89‑1 (MQ=255)
                               gCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTg                                                                            <  2:453937/89‑1 (MQ=255)
                               gCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTg                                                                            <  2:2897220/89‑1 (MQ=255)
                               gCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAc                                                                 >  1:1735408/1‑100 (MQ=255)
                                  cgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCgg                                                              >  1:150897/1‑100 (MQ=255)
                                    cgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                            <  1:3008767/100‑1 (MQ=255)
                                    cgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                            <  1:3160568/100‑1 (MQ=255)
                                       cgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCa                                                                    <  2:2634961/89‑1 (MQ=255)
                                        gcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCAt                                                                   >  2:2577341/1‑89 (MQ=255)
                                         cgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATa                                                                  >  2:4450382/1‑89 (MQ=255)
                                           aaaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAcc                                                                >  2:3084683/1‑89 (MQ=255)
                                             aGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCgg                                                              <  2:4596457/89‑1 (MQ=255)
                                              gTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTattatt                                                  <  1:4539633/100‑1 (MQ=255)
                                              gTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                            <  2:4318692/89‑1 (MQ=255)
                                               tCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                            <  2:3214063/89‑1 (MQ=255)
                                                cACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgc                                                >  1:3910609/1‑100 (MQ=255)
                                                 aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgct                                               <  1:2649957/100‑1 (MQ=255)
                                                 aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgct                                               <  1:2773913/100‑1 (MQ=255)
                                                  cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGt                                                         >  2:4523037/1‑89 (MQ=255)
                                                  cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgctt                                              <  1:1387755/100‑1 (MQ=255)
                                                   tGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtt                                                        <  2:1652812/89‑1 (MQ=255)
                                                   tGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgcttg                                             <  1:3035850/100‑1 (MQ=255)
                                                    gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtta                                                       >  2:3008523/1‑89 (MQ=255)
                                                     aTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttat                                                      <  2:686673/89‑1 (MQ=255)
                                                     aTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttat                                                      <  2:346583/89‑1 (MQ=255)
                                                     aTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCa                                           <  1:4303773/100‑1 (MQ=255)
                                                      tAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCAc                                          <  1:861492/100‑1 (MQ=255)
                                                       aaTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttatta                                                    >  2:283955/1‑89 (MQ=255)
                                                        aTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTg                                        >  1:1409647/1‑100 (MQ=255)
                                                         tGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTattatt                                                  <  2:3677647/89‑1 (MQ=255)
                                                         tGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTgg                                       >  1:4866186/1‑100 (MQ=255)
                                                          gcgcCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttg                                                 <  2:2771176/89‑1 (MQ=255)
                                                          gcgcCGCGTTAATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGc                                      >  1:2333348/1‑100 (MQ=255)
                                                           cgccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCt                                     <  1:3201158/100‑1 (MQ=255)
                                                            gccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgct                                               <  2:3137817/89‑1 (MQ=255)
                                                             ccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTTCTTGCACTGGCttt                                   >  1:2032742/1‑100 (MQ=255)
                                                               gcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgcttgc                                            <  2:829156/89‑1 (MQ=255)
                                                               gcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                 <  1:4264082/100‑1 (MQ=255)
                                                               gcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                 <  1:957892/100‑1 (MQ=255)
                                                               gcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                 <  1:3598004/100‑1 (MQ=255)
                                                                cgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCa                                           <  2:2259299/89‑1 (MQ=255)
                                                                  ttCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGc                              >  1:2395694/1‑100 (MQ=255)
                                                                   tCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTg                                        >  2:202074/1‑89 (MQ=255)
                                                                    cATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTgg                                       >  2:982666/1‑89 (MQ=255)
                                                                     aTGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGc                                      <  2:1031370/89‑1 (MQ=255)
                                                                      tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCt                                     <  2:4704243/89‑1 (MQ=255)
                                                                      tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGtt                          >  1:1582659/1‑100 (MQ=255)
                                                                      tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGtt                          >  1:888594/1‑100 (MQ=255)
                                                                       gTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtt                                    >  2:2377932/1‑89 (MQ=255)
                                                                         ccTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtttt                                  >  2:2435365/1‑89 (MQ=255)
                                                                         ccTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtttt                                  >  2:549373/1‑89 (MQ=255)
                                                                         ccTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTa                       <  1:2969239/100‑1 (MQ=255)
                                                                          ctaaaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                 >  2:383940/4‑89 (MQ=255)
                                                                          cTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTAc                      >  1:1889803/1‑100 (MQ=255)
                                                                          cTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTAc                      <  1:3826946/100‑1 (MQ=255)
                                                                           tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgt                                <  2:3390892/89‑1 (MQ=255)
                                                                            cAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                               >  2:3434906/1‑89 (MQ=255)
                                                                            cAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                               >  2:357864/1‑89 (MQ=255)
                                                                            cAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGc                              >  2:1412204/1‑89 (MQ=255)
                                                                             agaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTATTATTATTGTGTGCACTGGCTTTTGTGCTGTTCTACGaa                   <  1:1714513/98‑1 (MQ=255)
                                                                               aaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                            >  2:2301990/1‑89 (MQ=255)
                                                                                aTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGt                           <  2:3521588/89‑1 (MQ=255)
                                                                                aTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGcc                <  1:3154538/100‑1 (MQ=255)
                                                                                 tGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTACATACCGGGCTGTTATTATTGCTTGCGCTGGCTTTTGTGCTGTTCTACGAAGCCa               >  1:1285167/1‑100 (MQ=255)
                                                                                   gCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCt                        <  2:325672/89‑1 (MQ=255)
                                                                                    cGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTa                       >  2:392085/1‑89 (MQ=255)
                                                                                    cGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTa                       >  2:2059330/1‑89 (MQ=255)
                                                                                     gTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCaa           <  1:1112534/100‑1 (MQ=255)
                                                                                     gTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCaa           <  1:4557638/100‑1 (MQ=255)
                                                                                     gTAACGTCCTATGCTTACAATTGTTTGCTGGCCTGGATTTGCATGCGTAGCTGTTACTATCGGTTGGCCTGGCTTTGGTGTTGTtataa                      <  2:1925979/89‑2 (MQ=255)
                                                                                      tAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACg                     <  2:894418/89‑1 (MQ=255)
                                                                                       aaCGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATg         >  1:4009992/1‑100 (MQ=255)
                                                                                        aCGTCCTATGCTTACATTTGTTTGCTTGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGaa                   >  2:4013033/1‑89 (MQ=255)
                                                                                         cGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAg                  <  2:1967823/89‑1 (MQ=255)
                                                                                          gTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGc                 <  2:4824677/89‑1 (MQ=255)
                                                                                              tATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATc             >  2:1225113/1‑89 (MQ=255)
                                                                                              tATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa  <  1:983670/100‑1 (MQ=255)
                                                                                              tATGCTTACATTTGTTTGCTGGCCTGGATTGGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATc             <  2:4272521/89‑1 (MQ=255)
                                                                                                    |                                                                                              
GGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCA  >  NC_000913/3637547‑3637739

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: