Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,080 G→C 5.6% 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 5.6% 244.8 / 5.9 71intergenic (‑223/+272)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (2/2);  total (34/37)
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.43e-01

TTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCC  >  NC_000913/914932‑915224
                                                                                                                                                    |                                                                                                                                                
ttGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATg                                                                                                                                                  >  4:673815/1‑149 (MQ=255)
  gCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGa                                                                                                                                                <  3:656691/149‑1 (MQ=255)
      ttAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtgtg                                                                                                                                            >  2:60090/1‑149 (MQ=255)
      ttAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtgtg                                                                                                                                            >  1:580048/1‑149 (MQ=255)
       tAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTgg                                                                                                                                           <  2:694609/149‑1 (MQ=255)
        aaaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGAGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTggg                                                                                                                                          <  2:453840/149‑1 (MQ=255)
         aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGc                                                                                                                                         <  4:384647/149‑1 (MQ=255)
         aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTCAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGc                                                                                                                                         >  1:107928/1‑149 (MQ=255)
          aaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCt                                                                                                                                        >  2:695012/1‑149 (MQ=255)
          aaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                        <  4:350331/149‑1 (MQ=255)
           aGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTg                                                                                                                                       >  2:97040/1‑149 (MQ=255)
            ggAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGa                                                                                                                                      <  3:337434/149‑1 (MQ=255)
            ggAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGa                                                                                                                                      <  1:88191/149‑1 (MQ=255)
              aaaTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAgg                                                                                                                                    <  1:587016/149‑1 (MQ=255)
                aTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGgc                                                                                                                                  <  2:354567/149‑1 (MQ=255)
                     aTGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGtt                                                                                                                             <  3:673815/149‑1 (MQ=255)
                       gCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGtttt                                                                                                                           <  3:680334/149‑1 (MQ=255)
                       gCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGtttt                                                                                                                           >  1:471045/1‑149 (MQ=255)
                          tAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                        >  2:202097/1‑149 (MQ=255)
                           aaGGAGCAGTGGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                       >  2:73338/1‑149 (MQ=255)
                           aaGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                       >  1:422068/1‑149 (MQ=255)
                           aaGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                       <  2:366165/149‑1 (MQ=255)
                              gAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGGGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGa                                                                                                                    <  1:86394/149‑1 (MQ=255)
                                    tggagTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                                                                                                              <  1:73338/147‑1 (MQ=255)
                                     agagTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTTGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                             >  4:127459/1‑149 (MQ=255)
                                     agagTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                             <  3:294301/149‑1 (MQ=255)
                                             cGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGccc                                                                                                     >  1:272763/1‑149 (MQ=255)
                                               ttttCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGa                                                                                                   >  3:346117/1‑149 (MQ=255)
                                                       tAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGccc                                                                                                     <  4:175592/139‑1 (MQ=255)
                                                       tAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGccc                                                                                                     >  3:175592/1‑139 (MQ=255)
                                                            atatTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTa                                                                                      >  3:102172/1‑149 (MQ=255)
                                                             tatTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                     <  1:695012/149‑1 (MQ=255)
                                                              atTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAg                                                                                    <  1:411162/149‑1 (MQ=255)
                                                                  ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAt                                                                                >  4:611289/1‑149 (MQ=255)
                                                                  ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAt                                                                                <  2:422068/149‑1 (MQ=255)
                                                                    gcgcTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATaa                                                                              <  1:21813/149‑1 (MQ=255)
                                                                       cTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTc                                                                           >  3:196679/1‑149 (MQ=255)
                                                                         ttAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGa                                                                         <  4:102172/149‑1 (MQ=255)
                                                                           aaTTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAgg                                                                       >  2:699964/1‑149 (MQ=255)
                                                                             ttACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                     >  3:162327/1‑149 (MQ=255)
                                                                             ttACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                     >  1:5233/1‑149 (MQ=255)
                                                                             ttACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                     >  1:308848/1‑149 (MQ=255)
                                                                             ttACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                     >  1:390382/1‑149 (MQ=255)
                                                                                       cACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgc                                                           <  3:127459/149‑1 (MQ=40)
                                                                                        aCCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                     <  4:301273/122‑1 (MQ=39)
                                                                                        aCCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                     >  3:301273/1‑122 (MQ=40)
                                                                                            aTGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGa                                                      <  3:405537/149‑1 (MQ=37)
                                                                                              gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtt                                                    <  2:5233/149‑1 (MQ=35)
                                                                                                tggtggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAtttt                                                  <  1:123388/149‑1 (MQ=35)
                                                                                                  gtggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAt                                                                                >  2:557407/1‑117 (MQ=33)
                                                                                                  gtggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAt                                                                                <  1:557407/117‑1 (MQ=33)
                                                                                                    ggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTc                                              <  2:341696/149‑1 (MQ=34)
                                                                                                       gTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGAccc                                                                                                >  1:545978/1‑96 (MQ=34)
                                                                                                       gTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGAccc                                                                                                <  2:545978/96‑1 (MQ=34)
                                                                                                        tCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCa                                          >  1:340471/1‑149 (MQ=35)
                                                                                                          aaaaTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCaga                                        <  2:497378/149‑1 (MQ=35)
                                                                                                                      gCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGc                            >  2:396814/1‑149 (MQ=32)
                                                                                                                            ccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                      >  4:300208/1‑149 (MQ=32)
                                                                                                                             cTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                     >  1:574050/1‑149 (MQ=32)
                                                                                                                                   gCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcggc               >  4:513435/1‑149 (MQ=32)
                                                                                                                                    ccGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGgcggcggcg              <  2:526111/149‑1 (MQ=32)
                                                                                                                                       gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt           <  3:339035/149‑1 (MQ=32)
                                                                                                                                       gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt           <  4:162327/149‑1 (MQ=32)
                                                                                                                                        ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg          <  1:64752/149‑1 (MQ=32)
                                                                                                                                        ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg          >  3:334626/3‑149 (MQ=12)
                                                                                                                                          cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc        <  2:580048/149‑1 (MQ=32)
                                                                                                                                            gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                     <  3:681603/134‑1 (MQ=32)
                                                                                                                                            gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                     >  4:681603/1‑134 (MQ=32)
                                                                                                                                              ccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc    >  4:259273/1‑149 (MQ=17)
                                                                                                                                                acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc  >  1:660731/1‑149 (MQ=32)
                                                                                                                                                  aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                        <  4:580163/77‑1 (MQ=35)
                                                                                                                                                  aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGAGGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAg                                                                        >  3:580163/1‑77 (MQ=21)
                                                                                                                                                    |                                                                                                                                                
TTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCC  >  NC_000913/914932‑915224

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