Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,279,865 G→A 6.1% intergenic (+37/‑115) 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,8650GA 6.1% 413.3 / 9.8 131intergenic (+37/‑115)ghxP/yjcEguanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (63/58);  new base A (3/5);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 131 total aligned reads displayed.

GACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACC  >  NC_000913/4279717‑4280012
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                            <  1:14820/150‑2 (MQ=255)
        tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                            <  2:2117023/150‑2 (MQ=255)
          cttGGTCGGCTGGTATCCGGCGAATGGCGCAAGGTGAAGAGCGGTACGGTCGTTATCGGCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcag                                                                                                                                          <  1:127183/149‑4 (MQ=255)
          aTTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                                                          <  2:1482438/150‑1 (MQ=255)
             ggTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAtttt                                                                                                                                       >  2:1917905/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)
                   cTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTAATCGGATTCTGAAAACGGTTGGAAATGGCTGCCAGTTTTTAattttcagca                                                                                                                                 >  2:214856/1‑144 (MQ=255)
                     ggTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcac                                                                                                                               <  2:563277/150‑1 (MQ=255)
                          ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgt                                                                                                                          >  2:1643836/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)
                             gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTTGGCTATCTAATCTGTCATCGGAGTCTGAAAACGGGTGGCAATGGCTGCCCGTGTTTATTTTCTCCGCACATtgggtg                                                                                                                       >  2:333986/1‑150 (MQ=255)
                             gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                                       <  2:1283013/150‑1 (MQ=255)
                              cGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGa                                                                                                                      <  1:1171281/150‑1 (MQ=255)
                                 aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCt                                                                                                                   <  1:2188534/150‑1 (MQ=255)
                                      gcAAGTTGAACATCGGTACGGTCGTTATCGCCGGGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGc                                                                                                              <  2:219298/150‑1 (MQ=255)
                                       cAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                                                             <  2:1119052/150‑1 (MQ=255)
                                          gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttt                                                                                                                                              >  1:116586/1‑114 (MQ=255)
                                          gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGttt                                                                                                                                              <  2:116586/114‑1 (MQ=255)
                                            tGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCTATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATAtt                                                                                                        >  1:51867/1‑150 (MQ=255)
                                             gAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                                       <  1:100858/150‑1 (MQ=255)
                                                cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                                    >  2:626966/1‑150 (MQ=255)
                                                cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                                    >  1:1506978/1‑150 (MQ=255)
                                                cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCAATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                                    <  2:1989416/150‑1 (MQ=255)
                                                 aTCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaa                                                                                                   <  2:14291/150‑1 (MQ=255)
                                                    ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCAGTTTTTATTTTCTCAGCACATTGTGTGAGCTTTTGCGATATTCTGaaaaa                                                                                                >  2:1319844/1‑150 (MQ=255)
                                                       aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                             <  2:1180372/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)
                                                                           gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                         <  2:1012145/150‑1 (MQ=255)
                                                                           gcTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                                                         >  1:1770916/1‑150 (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)
                                                                                           gCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAATCAGGCATAACGTCTTTCTCAGATCATaaa                                                         >  1:1342616/1‑150 (MQ=255)
                                                                                            cGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGTCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaa                                                        <  1:62878/150‑1 (MQ=255)
                                                                                             gggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaaa                                                       <  1:578024/150‑1 (MQ=255)
                                                                                                  ggtGGGCTATCTAATGTTTCATCGGATTGAGTAAACGGGTGGCTATGGCTGCCCGTTTTTATTTTCTCCGCAATTTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGTGAATTGAGGCATAACGTCTTTCTCAGATCATAAAAAACATc                                                  <  2:1390590/148‑1 (MQ=255)
                                                                                                  gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACAta                                                  >  1:390395/1‑149 (MQ=255)
                                                                                                    tGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGc                                                >  2:2073398/1‑150 (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)
                                                                                                             tAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGa                                       <  2:913494/150‑1 (MQ=255)
                                                                                                              aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGaa                                      <  2:1237188/150‑1 (MQ=255)
                                                                                                              aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCAGAAAAAAATGAGAATTCAGGAATAACGTCTTTCTAAGATCATAAAAAACATAGCAAACAGGGaa                                      >  2:1278/1‑150 (MQ=255)
                                                                                                               aTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAAc                                     >  1:416545/1‑150 (MQ=255)
                                                                                                                 cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                                   >  2:1504236/1‑150 (MQ=255)
                                                                                                                 cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                                   <  1:1592683/150‑1 (MQ=255)
                                                                                                                 cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                                   <  1:552916/150‑1 (MQ=255)
                                                                                                                 cTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGc                                   >  2:105179/1‑150 (MQ=255)
                                                                                                                     cATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATgg                               >  1:2115853/1‑150 (MQ=255)
                                                                                                                     cATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAGTTTCTCCGCACATTTTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTATCAGATCAGAAAAAACAGCGCAAACAGTGAACGCAggg                               >  1:1501405/1‑150 (MQ=255)
                                                                                                                     aataggATTCTGAAAACGGGTGGCAATGGCTGCCAGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATgg                               <  2:1712108/146‑1 (MQ=255)
                                                                                                                      aTCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                              >  1:1789918/1‑150 (MQ=255)
                                                                                                                      aTCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                              <  2:1770916/150‑1 (MQ=255)
                                                                                                                      aTCGGATTCTGAAAACGGGTGGCAATGGATGCCCATTTTTATTTTATCCGCAAATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTATTTCTCAGATCATAAAAAACATCGCAAACAGGTAACGCATGGa                              >  2:1210543/1‑150 (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)
                                                                                                                           aTTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtct                         <  1:1547373/150‑1 (MQ=255)
                                                                                                                           aTTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtct                         <  1:1555195/150‑1 (MQ=255)
                                                                                                                             aataaaaaCGGGCAGCCATTGCCACCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTATCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttc                       >  2:881417/5‑150 (MQ=255)
                                                                                                                             aataaaaaCGGGCAGCCATTGCCACCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAGGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttc                       >  1:1675787/5‑150 (MQ=255)
                                                                                                                              cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct                      >  2:889158/1‑150 (MQ=255)
                                                                                                                              cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct                      <  1:939925/150‑1 (MQ=255)
                                                                                                                               tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGCGCTTTTGCGATATTCTGAAAAAAATTAGAATTCAGGCATTACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt                     >  1:434746/1‑150 (MQ=255)
                                                                                                                               tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt                     <  2:1395401/150‑1 (MQ=255)
                                                                                                                               tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt                     >  2:1953323/1‑150 (MQ=255)
                                                                                                                                gAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCttcttc                    <  1:1006940/150‑1 (MQ=255)
                                                                                                                                gAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGAAAAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCttattc                    >  1:1860600/1‑150 (MQ=255)
                                                                                                                                  aaaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCAc                  >  1:1858471/1‑150 (MQ=255)
                                                                                                                                   aaCGGGTGGCTATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCAcc                 <  2:51867/150‑1 (MQ=255)
                                                                                                                                         gggCAATGGCTGCCCGTTTTTATTTTCTACGTACATTGTGTGAGCTTTTGAGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACtt           <  1:1582595/149‑2 (MQ=255)
                                                                                                                                          ggCAATGGCTGCCCGTTTTTATTTTCTCAGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGATTTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAAAAGGGAACGCATGGAAATCTTCTTCACCTTCCTGa          >  2:944038/1‑150 (MQ=255)
                                                                                                                                           gCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGaa         >  1:812135/1‑149 (MQ=255)
                                                                                                                                              aTGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAAt      <  1:353193/150‑1 (MQ=255)
                                                                                                                                                ggCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGa    >  1:1631245/1‑150 (MQ=255)
                                                                                                                                                 gCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAACCATCGCAACCAGGGCCCGCATGGAATTCTGCTTCACCATACTGCTaaagcc   >  2:317529/1‑146 (MQ=255)
                                                                                                                                                  cTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGAcc  >  2:564441/1‑150 (MQ=255)
                                                                                                                                                  cTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGAcc  <  1:977337/150‑1 (MQ=255)
                                                                                                                                                  cTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGAcc  >  2:1739247/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATAATGACC  >  NC_000913/4279717‑4280012

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: