Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,081 G→C 7.2% intergenic (‑224/+271) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,0810GC 7.2% 175.2 / 6.6 69intergenic (‑224/+271)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 = 4.54e-01

GCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGC  >  NC_000913/914934‑915230
                                                                                                                                                   |                                                                                                                                                     
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                                                                                     >  1:1213868/1‑150 (MQ=255)
                                                                ttGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATa                                                                                     <  2:937182/150‑1 (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:1123673/150‑1 (MQ=32)
                                                                                                                         cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                            >  1:874847/1‑150 (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)
                                                                                                                                                   gATGGGGGCGGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc  <  1:1062146/150‑1 (MQ=17)
                                                                                                                                                   |                                                                                                                                                     
GCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGC  >  NC_000913/914934‑915230

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: