Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,080 G→C 7.3% 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 7.3% 173.5 / 8.0 69intergenic (‑223/+272)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base G (22/41);  new base C (2/3);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-02

GTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCC  >  NC_000913/914931‑915219
                                                                                                                                                     |                                                                                                                                           
gTTGCAATTAAAAGTAAATGTTATGCATAAGGAGCAGGAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATg                                                                                                                                             <  1:1157972/150‑1 (MQ=255)
   gCAATTAAAAGGAAATGTTATGCATAAGGAGAAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAAATGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCAACATGGAt                                                                                                                                          >  2:212756/1‑150 (MQ=255)
    cAATTAAAAGGACATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGCTGTTGGGGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAtg                                                                                                                                         <  1:698733/150‑1 (MQ=255)
        tAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTggg                                                                                                                                     <  1:875238/150‑1 (MQ=255)
         aaaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGc                                                                                                                                    <  1:282295/150‑1 (MQ=255)
         aaaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGc                                                                                                                                    >  1:1264645/1‑150 (MQ=255)
          aaatgcaaTTTTATGCATAAGTAGCAGTAGAGTATTCGTTATCATTTAACGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAATATCGCCGAGGCATTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                   <  2:108356/144‑1 (MQ=37)
          aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCt                                                                                                                                   >  1:830276/1‑150 (MQ=255)
          aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCt                                                                                                                                   >  1:355450/1‑150 (MQ=255)
           aaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTg                                                                                                                                  <  1:1016136/150‑1 (MQ=255)
           aaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTg                                                                                                                                  <  2:470199/150‑1 (MQ=255)
           aaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCTg                                                                                                                                  <  1:432697/150‑1 (MQ=255)
            aGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCTgc                                                                                                                                 <  1:1210699/150‑2 (MQ=255)
             ggAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                <  1:723511/150‑1 (MQ=255)
             ggAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATAATCTTGCGCTTTAATAACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                <  2:970586/150‑1 (MQ=255)
              gAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAgg                                                                                                                               <  1:514092/150‑1 (MQ=255)
                aaTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGgc                                                                                                                             <  2:375935/150‑1 (MQ=255)
                aaTGGTATGCATAAGGAGCAGTAGACTATTCGTTGTCATTTAAAGATCTGCTTGCGCATTAATTACAAACTGCACCGAGGTGGGTGGCGGCAAACTCGCCGAGGCGGTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGgc                                                                                                                             <  1:162859/150‑1 (MQ=25)
                 aTGTTGTGCATAAGGCGCAGTAGTGTATTCTTTTTCATTTAATTATATTCTTGCTCGTTCATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGTGGCGGTCCCTGCAGGCCGGGGGAGCCCACATGGATGGGGGCTGAGGgcg                                                                                                                            <  2:661982/150‑1 (MQ=255)
                          aTAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACa                                                                                                                   >  2:26183/1‑150 (MQ=255)
                           tAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                  <  1:333401/150‑1 (MQ=255)
                             agcaGCAGTAGAGTATGCGTGTTCATTTTAAGCTATTCTTGCACTTTAATTACAAACTCCACCGAATTTGCTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGTATGTGGGCTGAGGGCGCGTATTACAggg                                                                                                                <  2:641765/147‑1 (MQ=18)
                              gtatcAGTAGAGTATTCGTGTTCATTTAAATATATTCTTGCGCTTTACTTACAAACTGCACCGATGGTGTGGGCGGCAATAGCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACTTGGATGTGGGCTGAGGGCGCGTTTTACAGGGa                                                                                                               <  2:209094/146‑1 (MQ=21)
                                aGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATg                                                                                                             <  2:557871/150‑1 (MQ=255)
                                  cAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGtt                                                                                                           <  1:995580/150‑1 (MQ=255)
                                  cAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGGTGGTGGCGTCAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTa                                                                                                          >  1:310381/1‑150 (MQ=255)
                                    gTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                                                                                                         >  2:1123673/1‑150 (MQ=255)
                                    gTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCCTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                                                                                                         <  2:98260/150‑1 (MQ=255)
                                    gTAGAGTATTCGTTTTCATTAAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                                                                                                         <  1:1099487/150‑1 (MQ=255)
                                           aTTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCAGGGGCAGCCCACATGGATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCgcgc                                                                                                  >  1:630153/1‑150 (MQ=255)
                                              cGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCg                                                                                               >  1:123102/1‑150 (MQ=255)
                                              cGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACAGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCACTGAAGGCAGGGGCAGACAACATGGATGTGGGCTGAGCGCGCGTTTTACAGGGATGTTACCTCGCGCCCg                                                                                               >  1:61784/1‑150 (MQ=21)
                                                       ttACAGAGATTCTTGCGCGTTAATTACAACCTGCCCCTATGTTGGGGGCGTCAAACTCGCCGAGGCGGTCCCTGCAGGCCGGGGCAGCCCACATGGACGTGGGCTGAGGGCGCGTTTTACAGTGATGTTACCTCGCGCCCGACCCGGGAg                                                                                      <  1:638444/150‑1 (MQ=17)
                                                       ttAAAGATATTCTTGCGCTTTAAATACATACTCCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAg                                                                                      <  1:638442/150‑1 (MQ=255)
                                                         aaaGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGTAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGcc                                                                                    <  2:511873/150‑1 (MQ=255)
                                                               atTCTTGCGCTTTAATTACAAACTGCACCGATGGTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATTGATGTGGGCTGAGGGCGCGTTTTACAGGGAGGTTACCTCGCGCCCGACCCGGTAGCCGTAAgg                                                                              <  1:938625/150‑1 (MQ=255)
                                                                   ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATa                                                                          <  2:937182/150‑1 (MQ=255)
                                                                   ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATa                                                                          >  1:1213868/1‑150 (MQ=255)
                                                                     gcgcTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAAGCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCCGTAGCCGTAAGGGATAAg                                                                        <  1:552955/150‑1 (MQ=255)
                                                                      cgcTTTAATTACCAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGACGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGt                                                                       <  2:763831/150‑1 (MQ=255)
                                                                                   aaCTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCAccaac                                                          >  1:338348/1‑147 (MQ=255)
                                                                                     cTGCACCGATGTTGGTGGCGTCAAAATCGCCGCGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTTGCCGTAAGGGATAAGTCGAGGGCACCGCGCa                                                        >  1:1061190/1‑150 (MQ=255)
                                                                                      tGCACCGCTGTAGGAGGCGTCAACCTCGGCGAGGCGTACCCTGAAGGCCGGGGCAGCCCACATGGAGGTGGGCTGAGGGCGCGTGTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAg                                                       <  1:193797/150‑1 (MQ=14)
                                                                                        cACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcg                                                     >  1:682013/1‑150 (MQ=40)
                                                                                        cACCGATGTTGGTGGCGTCAAAATCGCAGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGTCTCGTTTTACAGGGATGTTACCTCGCGCCAGACACGATAGCCGTAAGTGATAAGTCGAGGGCAACGCGCAgcg                                                     >  2:467194/1‑150 (MQ=14)
                                                                                         aCCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcgg                                                    <  2:302007/150‑1 (MQ=40)
                                                                                          ccGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcggc                                                   >  2:1256358/1‑150 (MQ=39)
                                                                                          ccGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCGGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcggc                                                   <  2:525367/150‑1 (MQ=255)
                                                                                               gTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGCAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAttt                                              <  1:1135109/150‑1 (MQ=25)
                                                                                                    tggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTc                                         <  1:1107622/150‑1 (MQ=34)
                                                                                                     ggcggcATACTCGCCGAGGCGTTCCCTGTAGGCCGGGGCAGCCCACATGGTTGGGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCAGGTAGCCGTAAGGGATAAGTCGAGGGCACCGAGCAGCGGCGATTTTGTTCg                                        <  2:132174/150‑1 (MQ=11)
                                                                                                       cGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGcc                                      >  2:82653/1‑150 (MQ=34)
                                                                                                         gcAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCag                                    <  1:976820/149‑1 (MQ=21)
                                                                                                                  ccGAGGCGTTCCCTGAAGTCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                           <  1:633150/150‑1 (MQ=17)
                                                                                                                  ccGAGGCGTTCCCTGAAGGCGGGGGCAGCCCACATGGATGTGTGCTGTGGGCGCGTTTTACAGGGATGTTACCTCGCGCGCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                           <  1:1083959/150‑1 (MQ=255)
                                                                                                                  ccGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggggg                          >  1:249914/1‑150 (MQ=17)
                                                                                                                         gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgg                    <  2:282313/150‑1 (MQ=32)
                                                                                                                            cccTGAAGGCCGTGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCCGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                 >  2:1224474/1‑150 (MQ=17)
                                                                                                                            cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                 >  1:874847/1‑150 (MQ=32)
                                                                                                                            cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                 <  1:1123673/150‑1 (MQ=32)
                                                                                                                              cTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcg               >  1:1010737/1‑150 (MQ=32)
                                                                                                                               tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg              >  1:448010/1‑150 (MQ=32)
                                                                                                                               tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg              <  2:851669/150‑1 (MQ=32)
                                                                                                                               tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg              <  2:855745/150‑1 (MQ=32)
                                                                                                                                 aaGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcg            <  2:1264645/150‑1 (MQ=32)
                                                                                                                                      cGGGGCAGCCCCCAGGGATGGGGGCGGAGGGCGCGTGTTTCAGGGACGTTACCGCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTGCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt       <  1:203686/150‑1 (MQ=11)
                                                                                                                                        gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg     <  2:1010737/150‑1 (MQ=32)
                                                                                                                                           cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc  <  2:830276/150‑1 (MQ=32)
                                                                                                                                           cAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc  >  2:254735/1‑150 (MQ=17)
                                                                                                                                           cAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGgcc  >  1:1222226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCC  >  NC_000913/914931‑915219

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