Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,279,858 A→C 6.1% intergenic (+30/‑122) ghxP → / → yjcE guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,8580AC 6.1% 374.9 / 9.6 127intergenic (+30/‑122)ghxP/yjcEguanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (59/57);  new base C (3/5);  total (63/64)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

GGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGA  >  NC_000913/4279709‑4280004
                                                                                                                                                     |                                                                                                                                                  
ggCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCa                                                                                                                                                    >  1:2210846/1‑150 (MQ=255)
ggCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCa                                                                                                                                                    >  2:62878/1‑150 (MQ=255)
   ttCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATg                                                                                                                                                 <  1:1684935/150‑1 (MQ=255)
      gcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCt                                                                                                                                              >  1:314303/1‑150 (MQ=255)
        gACTCTGGGGATTGGTCGTCTGGTTTCCGGCGAATGGCGAAAGTTGAACATCTGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGTGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGc                                                                                                                                            <  1:1957927/150‑1 (MQ=255)
          ctctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGccc                                                                                                                                          >  2:1173410/1‑150 (MQ=255)
             tGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGtt                                                                                                                                       >  2:939925/1‑150 (MQ=255)
              ggTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttt                                                                                                                                      >  2:2106658/1‑150 (MQ=255)
                tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                    <  2:85863/150‑2 (MQ=255)
                tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                    <  2:2117023/150‑2 (MQ=255)
                tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                    <  1:14820/150‑2 (MQ=255)
                  cttGGTCGGCTGGTATCCGGCGAATGGCGCAAGGTGAAGAGCGGTACGGTCGTTATCGGCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcag                                                                                                                                  <  1:127183/149‑4 (MQ=255)
                  aTTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                                                  <  2:1482438/150‑1 (MQ=255)
                     ggTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTTGCAATGGCTGCCCGTTTTTAtttt                                                                                                                               >  2:1852695/1‑150 (MQ=255)
                        cgtcTGGTTTCCGGAGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtctc                                                                                                                            >  1:1686805/1‑150 (MQ=255)
                         gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCGTCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTcc                                                                                                                           <  2:630984/150‑1 (MQ=255)
                         gacTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTCGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCACGTTTTTATTTTATcc                                                                                                                           >  2:462757/3‑150 (MQ=255)
                           cTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGc                                                                                                                         >  1:1824394/1‑150 (MQ=255)
                            tGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGca                                                                                                                        >  1:1608992/1‑150 (MQ=255)
                             ggTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcac                                                                                                                       <  2:563277/150‑1 (MQ=255)
                                  ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgt                                                                                                                  >  2:1643836/1‑150 (MQ=255)
                                  ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGAGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCAAGTTTTTATTTTCTCCGCACATtgt                                                                                                                  >  2:131590/1‑150 (MQ=255)
                                    ggCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTTGTGACATTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                                                <  2:1082040/150‑1 (MQ=255)
                                    ggCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                                                >  1:1434509/1‑150 (MQ=255)
                                    ggCGAATGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                               <  1:1050058/150‑1 (MQ=255)
                                      cGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGa                                                                                                              <  1:1171281/150‑1 (MQ=255)
                                       gAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAg                                                                                                             <  1:422287/150‑1 (MQ=255)
                                             cgcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCACCGCACATTTTGTGAGCTTTTg                                                                                                       >  1:215699/1‑150 (MQ=255)
                                              gcAAGTTGAACATCGGTACGGTCGTTATCGCCGGGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                      <  2:219298/150‑1 (MQ=255)
                                               cAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                                                     <  2:1119052/150‑1 (MQ=255)
                                                  gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttt                                                                                                                                      <  2:116586/114‑1 (MQ=255)
                                                  gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttt                                                                                                                                      >  1:116586/1‑114 (MQ=255)
                                                    tGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCTATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATAtt                                                                                                >  1:51867/1‑150 (MQ=255)
                                                     gAACATCGGTTCGGTCGTTATCGCCGTGGCGCGGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                               <  2:1396370/150‑1 (MQ=255)
                                                     gAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                               <  1:100858/150‑1 (MQ=255)
                                                     gAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                               <  1:1490164/150‑1 (MQ=255)
                                                        cATGGGTACGGTCGTAATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGTTTCTGAAAACGGGTGGCAATGGCTTCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                            <  2:1073785/150‑1 (MQ=255)
                                                        cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                            >  1:1506978/1‑150 (MQ=255)
                                                        cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                            >  2:626966/1‑150 (MQ=255)
                                                        cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCAATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                            <  2:1989416/150‑1 (MQ=255)
                                                         aTCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaa                                                                                           <  2:14291/150‑1 (MQ=255)
                                                         aTCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaa                                                                                           >  1:1707890/1‑150 (MQ=255)
                                                            ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCAGTTTTTATTTTCTCAGCACATTGTGTGAGCTTTTGCGATATTCTGaaaaa                                                                                        >  2:1319844/1‑150 (MQ=255)
                                                               aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                     <  2:1180372/150‑1 (MQ=255)
                                                               aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                     <  2:1172636/150‑1 (MQ=255)
                                                                    cGTTATCGCCGTGGCGCTGGTGTACTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTTCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGaa                                                                                <  1:416900/150‑1 (MQ=255)
                                                                    cGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGaa                                                                                <  2:314303/150‑1 (MQ=255)
                                                                     gTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAt                                                                               >  1:2177585/1‑150 (MQ=255)
                                                                        aTCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGCGCTTTTGCGATATTCTGAAAAAAATGAGAATTCa                                                                            >  2:637232/1‑150 (MQ=255)
                                                                          cGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCGGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAgg                                                                          <  1:928967/150‑1 (MQ=255)
                                                                           gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                         >  1:15327/1‑150 (MQ=255)
                                                                           gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGCGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                         >  2:210531/1‑150 (MQ=255)
                                                                            ccGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCa                                                                        <  2:1894457/150‑1 (MQ=255)
                                                                                   gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCTGATTCTGATAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                 <  2:2115592/150‑1 (MQ=255)
                                                                                   gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                 <  2:1012145/150‑1 (MQ=255)
                                                                                   gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                 >  1:1770916/1‑150 (MQ=255)
                                                                                      ggTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCttt                                                              >  2:1230406/1‑150 (MQ=255)
                                                                                      ggTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGAAATGGCTGCCCGTTTTTATTTTCTCCGAACATTGTTTGAGCTTTTGCGATATTCTGAATAAAATGAGAATTCAGGCATAACGTCttt                                                              >  1:362957/1‑150 (MQ=255)
                                                                                         gACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctc                                                           <  2:2201031/150‑1 (MQ=255)
                                                                                           ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                                                         <  1:576491/150‑1 (MQ=255)
                                                                                             ttCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGAt                                                       >  1:1222187/1‑150 (MQ=255)
                                                                                               cTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCa                                                     <  2:499114/150‑1 (MQ=255)
                                                                                                 aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATa                                                   <  2:1682668/150‑1 (MQ=255)
                                                                                                  tGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaa                                                  <  1:1288592/150‑1 (MQ=255)
                                                                                                   gCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaa                                                 <  1:965906/150‑1 (MQ=255)
                                                                                                    cGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGTCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaa                                                <  1:62878/150‑1 (MQ=255)
                                                                                                      ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTATGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaaaa                                              >  2:369279/1‑150 (MQ=255)
                                                                                                          ggtGGGCTATCTAATGTTTCATCGGATTGAGTAAACGGGTGGCTATGGCTGCCCGTTTTTATTTTCTCCGCAATTTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGTGAATTGAGGCATAACGTCTTTCTCAGATCATAAAAAACATc                                          <  2:1390590/148‑1 (MQ=255)
                                                                                                             gggCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCa                                       >  1:1819689/1‑150 (MQ=255)
                                                                                                               gCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCCCATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCaaa                                     >  1:1186147/1‑150 (MQ=255)
                                                                                                               gCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCaaa                                     >  2:1658837/1‑150 (MQ=255)
                                                                                                                cTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAAc                                    >  1:1023948/1‑150 (MQ=255)
                                                                                                                  aTCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg                                  >  2:2025186/1‑150 (MQ=255)
                                                                                                                     tAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGa                               <  2:913494/150‑1 (MQ=255)
                                                                                                                      aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGaa                              <  2:1237188/150‑1 (MQ=255)
                                                                                                                      aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATCCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGaa                              >  2:532400/1‑150 (MQ=255)
                                                                                                                      aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCAGAAAAAAATGAGAATTCAGGAATAACGTCTTTCTAAGATCATAAAAAACATAGCAAACAGGGaa                              >  2:1278/1‑150 (MQ=255)
                                                                                                                       aTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAAc                             >  1:416545/1‑150 (MQ=255)
                                                                                                                         cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                           <  1:552916/150‑1 (MQ=255)
                                                                                                                         cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                           <  1:1592683/150‑1 (MQ=255)
                                                                                                                             cATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATgg                       >  1:2115853/1‑150 (MQ=255)
                                                                                                                             aataggATTCTGAAAACGGGTGGCAATGGCTGCCAGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATgg                       <  2:1712108/146‑1 (MQ=255)
                                                                                                                              aTCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                      <  2:1770916/150‑1 (MQ=255)
                                                                                                                                 ggATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAATAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAt                   >  2:423920/1‑150 (MQ=255)
                                                                                                                                  gaaaataaaaaCGGGCAGCCATTGCCACCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtc                  >  2:1286847/8‑150 (MQ=255)
                                                                                                                                  gATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtc                  <  2:552334/150‑1 (MQ=255)
                                                                                                                                  gATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtc                  <  1:756395/150‑1 (MQ=255)
                                                                                                                                   agtCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCCGATCATAAAAACCATCGCAAACAGGGAACGCATGGAActct                 <  1:1547378/148‑1 (MQ=255)
                                                                                                                                   aTTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtct                 <  1:1547373/150‑1 (MQ=255)
                                                                                                                                     aataaaaaCGGGCAGCCATTGCCACCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTATCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttc               >  2:881417/5‑150 (MQ=255)
                                                                                                                                      cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct              <  1:939925/150‑1 (MQ=255)
                                                                                                                                      cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct              >  2:889158/1‑150 (MQ=255)
                                                                                                                                       tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt             >  2:1953323/1‑150 (MQ=255)
                                                                                                                                        gAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCttcttc            <  1:1006940/150‑1 (MQ=255)
                                                                                                                                         aaaaCGGGTGGCAATGGCTGCCCGTTTTTAATTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTACACAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCa           >  1:1332853/1‑150 (MQ=255)
                                                                                                                                          aaaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCAc          >  1:1858471/1‑150 (MQ=255)
                                                                                                                                             cGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCAt       <  1:563791/150‑1 (MQ=255)
                                                                                                                                                 gggCAATGGCTGCCCGTTTTTATTTTCTACGTACATTGTGTGAGCTTTTGAGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACt  <  1:1582595/149‑2 (MQ=255)
                                                                                                                                                  ggCAATGGCTGCCCGTTTTTATTTTCTCCTCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGa  <  1:648744/150‑1 (MQ=255)
                                                                                                                                                  ggCAATGGCTGCCCGTTTTTATTTTCTCAGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGATTTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAAAAGGGAACGCATGGAAATCTTCTTCACCTTCCTGa  >  2:944038/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGA  >  NC_000913/4279709‑4280004

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