Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,080 G→C 8.0% intergenic (‑223/+272) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,0800GC 8.0% 423.2 / 24.8 125intergenic (‑223/+272)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base G (40/75);  new base C (3/7);  total (43/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

Only 100 of 128 total aligned reads displayed.

TTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT  >  NC_000913/914932‑915226
                                                                                                                                                    |                                                                                                                                                  
    aaTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtg                                                                                                                                                >  3:982141/1‑149 (MQ=255)
    aaTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtg                                                                                                                                                >  2:177565/1‑149 (MQ=255)
    aaTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTGAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCTAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtg                                                                                                                                                >  2:96125/1‑149 (MQ=255)
     aTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtgt                                                                                                                                               <  3:747180/149‑1 (MQ=255)
     aTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtgt                                                                                                                                               >  3:1370028/1‑149 (MQ=255)
     aTTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtgt                                                                                                                                               <  2:1452528/149‑1 (MQ=255)
       tAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTgg                                                                                                                                             <  2:1023533/149‑1 (MQ=255)
       tAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTgg                                                                                                                                             <  4:1447512/149‑1 (MQ=255)
         aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGc                                                                                                                                           <  4:1562693/149‑1 (MQ=255)
         aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGc                                                                                                                                           <  3:396476/149‑1 (MQ=255)
         aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGc                                                                                                                                           <  3:1264215/149‑1 (MQ=255)
           aGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTg                                                                                                                                         >  1:531877/1‑149 (MQ=255)
            ggAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGa                                                                                                                                        <  1:229895/149‑1 (MQ=255)
             gAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                       >  1:366071/1‑149 (MQ=255)
             gAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                          <  2:684331/146‑1 (MQ=255)
             gAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                          >  1:684331/1‑146 (MQ=255)
              aaaTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAgg                                                                                                                                      <  1:109199/149‑1 (MQ=255)
              aaaTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCTgccc                                                                                                                                      <  1:1064071/149‑4 (MQ=255)
               aaTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAggg                                                                                                                                     <  2:1205863/149‑1 (MQ=255)
                aTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGgc                                                                                                                                    <  2:431653/149‑1 (MQ=255)
                    tATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGt                                                                                                                                <  2:304472/149‑1 (MQ=255)
                     aTGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGtt                                                                                                                               >  3:680325/1‑149 (MQ=255)
                      tGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGttt                                                                                                                              <  1:792190/149‑1 (MQ=255)
                        cATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTa                                                                                                                            >  2:375631/1‑149 (MQ=255)
                           aaGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                         <  1:381876/149‑1 (MQ=255)
                           aaGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                         <  1:1561015/149‑1 (MQ=255)
                           aaGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                         <  3:237903/149‑1 (MQ=255)
                            aGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAgg                                                                                                                        <  2:310074/149‑1 (MQ=255)
                             ggAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAgg                                                                                                                                      >  2:832171/1‑134 (MQ=255)
                             ggAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAgg                                                                                                                                      <  1:832171/134‑1 (MQ=255)
                               aGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                                     <  2:1090520/149‑1 (MQ=255)
                               aGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                                     <  3:634059/149‑1 (MQ=255)
                               aGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                                     <  3:599309/149‑1 (MQ=255)
                               aGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                                     <  3:34669/149‑1 (MQ=255)
                                     agagTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                               <  1:827870/149‑1 (MQ=255)
                                     agagTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                               <  4:246825/149‑1 (MQ=255)
                                         tATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                          <  1:1431741/134‑1 (MQ=255)
                                         tATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                          >  2:1431741/1‑134 (MQ=255)
                                         tATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTcgc                                                                                                           <  4:1370028/149‑1 (MQ=255)
                                             cGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAAACGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGccc                                                                                                       <  4:888759/149‑1 (MQ=255)
                                               ttttCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                               <  4:1053786/139‑1 (MQ=255)
                                                    aTTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCgg                                                                                                <  1:883665/149‑1 (MQ=255)
                                                     tttAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                          <  1:735715/122‑1 (MQ=255)
                                                     tttAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                          >  2:735715/1‑122 (MQ=255)
                                                         aaGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGcc                                                                                           <  1:20775/149‑1 (MQ=255)
                                                          aGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAg                                                                                             <  2:157595/146‑1 (MQ=255)
                                                          aGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAg                                                                                             >  1:157595/1‑146 (MQ=255)
                                                               ttCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAgg                                                                                     >  4:879351/1‑149 (MQ=255)
                                                               ttCTTGCGCTTTAATTACAAACTGCACCGATGATGGTGGCGACAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCCCTCCCCCGACCCGGTAGCCGTAAgg                                                                                     <  4:1281359/149‑1 (MQ=17)
                                                                tCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAggg                                                                                    <  2:738701/149‑1 (MQ=255)
                                                                  ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAt                                                                                  <  1:375631/149‑1 (MQ=255)
                                                                      gcTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGt                                                                              >  4:427941/1‑149 (MQ=255)
                                                                         ttAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGa                                                                           <  2:538942/149‑1 (MQ=255)
                                                                          tAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  3:879351/149‑1 (MQ=255)
                                                                          tAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  3:278849/149‑1 (MQ=255)
                                                                          tAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  2:409108/149‑1 (MQ=255)
                                                                          tAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  4:680325/149‑1 (MQ=255)
                                                                           aaTTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAgg                                                                         <  3:207229/149‑1 (MQ=255)
                                                                             ttACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                       >  2:1095633/1‑149 (MQ=255)
                                                                                cAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGt                                                                              <  4:576518/139‑1 (MQ=255)
                                                                                 aaaCTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACcg                                                                   >  2:845389/1‑149 (MQ=255)
                                                                                         ccGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcgg                                                           <  1:1215265/149‑1 (MQ=40)
                                                                                             tGTTGGCGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAt                                                       <  4:27443/149‑1 (MQ=37)
                                                                                              gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      <  2:118920/149‑1 (MQ=35)
                                                                                              gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      <  2:938766/149‑1 (MQ=35)
                                                                                              gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      <  1:1095633/149‑1 (MQ=35)
                                                                                              gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      <  1:1236316/149‑1 (MQ=35)
                                                                                                tggtggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtttt                                                    <  2:1227370/149‑1 (MQ=35)
                                                                                                   tggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                                 <  2:1589915/149‑1 (MQ=35)
                                                                                                   tggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                                 <  3:858686/149‑1 (MQ=35)
                                                                                                   tggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                                 <  1:177565/149‑1 (MQ=34)
                                                                                                        tCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCa                                            >  2:505353/1‑149 (MQ=34)
                                                                                                        tCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCa                                            <  3:427941/149‑1 (MQ=34)
                                                                                                          aaaaTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCaga                                          >  1:632006/1‑149 (MQ=35)
                                                                                                             aTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGcc                                       >  4:563670/1‑149 (MQ=34)
                                                                                                                gCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgg                                    <  1:1003059/149‑1 (MQ=32)
                                                                                                                  cGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                  <  2:257140/149‑1 (MQ=32)
                                                                                                                  cGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                  <  1:1314898/149‑1 (MQ=32)
                                                                                                                      gCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCg                                               >  1:1518021/1‑132 (MQ=32)
                                                                                                                      gCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCg                                               <  2:1518021/132‑1 (MQ=33)
                                                                                                                        gTTCCCTGATGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                  >  3:1096835/1‑143 (MQ=17)
                                                                                                                        gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                  <  4:1096835/143‑1 (MQ=32)
                                                                                                                         ttCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                                 >  4:1451815/1‑127 (MQ=32)
                                                                                                                         ttCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgg                           >  1:158537/1‑149 (MQ=32)
                                                                                                                              tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  1:777482/97‑1 (MQ=33)
                                                                                                                              tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcggc                                                          >  1:1283886/1‑113 (MQ=32)
                                                                                                                              tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcggc                                                          <  2:1283886/113‑1 (MQ=33)
                                                                                                                               gAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                     >  3:226417/1‑149 (MQ=32)
                                                                                                                                     cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCACGAACCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCAAGAGCCCGGGGGTGCAGGGGGCGGCGGCGa               <  1:96125/149‑1 (MQ=12)
                                                                                                                                       gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt             <  1:367122/149‑1 (MQ=32)
                                                                                                                                        ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          <  4:978161/85‑1 (MQ=11)
                                                                                                                                        ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                          >  3:978161/3‑87 (MQ=11)
                                                                                                                                           aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc         >  2:419887/1‑149 (MQ=32)
                                                                                                                                            gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      <  2:7706/103‑1 (MQ=33)
                                                                                                                                            gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                      >  1:7706/1‑103 (MQ=32)
                                                                                                                                            gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg        <  1:1308237/149‑1 (MQ=32)
                                                                                                                                            gCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                       <  1:1373486/134‑1 (MQ=255)
                                                                                                                                            gCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                       >  2:1373486/1‑134 (MQ=255)
                                                                                                                                              ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc      >  2:619056/1‑149 (MQ=32)
                                                                                                                                                  aTGGATGTGGGCTGAGGGCGCGGTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt  <  4:226417/149‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                  
TTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT  >  NC_000913/914932‑915226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: