Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,782 C→A 5.3% A312A (GCG→GCT nhaB ← Na(+):H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,7820CA 5.3% 254.3 / 2.6 75A312A (GCG→GCTnhaBNa(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base A (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAACGC  >  NC_000913/1233633‑1233929
                                                                                                                                                      |                                                                                                                                                   
gagaAAAACACCGTCAACAGTGCGGTGAATTGCAGAGCTACGGTGAACGCTTTACCGATAGCATGCAAAACGCTGACACCCGTCAATGATGTTCGCAGAATAATGACTCACAAACCAATCAACCCAACATCCGCCAGATGCAACCGCAgc                                                                                                                                                     >  2:465851/1‑150 (MQ=255)
gagaAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAgc                                                                                                                                                     >  2:53821/1‑150 (MQ=255)
     aaaCACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGt                                                                                                                                                >  1:616483/1‑150 (MQ=255)
           cGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccag                                                                                                                                          <  1:401348/150‑1 (MQ=255)
              cAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                                                                                                                                       >  2:444836/1‑150 (MQ=255)
              cAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAAAGATGTTGCCAGAATAATGACTGACAAACCAATAAAGCCAGCTTCCGCGAGATGAAACGCCAGCGCAGTCAccagcca                                                                                                                                       >  2:493667/1‑150 (MQ=255)
                   gTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGc                                                                                                                                  <  2:673464/150‑1 (MQ=255)
                   gTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGc                                                                                                                                  >  1:937527/1‑150 (MQ=255)
                     gCGGTTAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGGTCATCGGTGACACCGGTCAATGATTTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCTGTCACCATCCAGACGCCa                                                                                                                                <  2:1095553/150‑1 (MQ=255)
                       ggTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaat                                                                                                                              <  1:162515/150‑1 (MQ=255)
                           aaTGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATc                                                                                                                          >  1:890241/1‑150 (MQ=255)
                              ggCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGcc                                                                                                                       >  2:113817/1‑150 (MQ=255)
                              ggCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGcc                                                                                                                       >  2:874192/1‑150 (MQ=255)
                              ggCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGcc                                                                                                                       >  2:976312/1‑150 (MQ=255)
                                 agagATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTgg                                                                                                                    <  2:1582755/150‑1 (MQ=255)
                                   agaTTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGAc                                                                                                                  <  1:679643/150‑1 (MQ=255)
                                    gaTTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACa                                                                                                                 >  1:214794/1‑150 (MQ=255)
                                      ttCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAAt                                                                                                               <  1:1107438/150‑1 (MQ=255)
                                        cGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCa                                                                                                             <  1:976312/150‑1 (MQ=255)
                                               ggCTTTACCTATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGAt                                                                                                      >  1:1096978/1‑150 (MQ=255)
                                                  gtcaCCGATGAGCATGCTATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGAtttt                                                                                                   >  1:165884/4‑150 (MQ=255)
                                                      ccGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATcc                                                                                               <  1:409069/150‑1 (MQ=255)
                                                              cATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACggg                                                                                        <  1:16083/150‑1 (MQ=255)
                                                              cATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCAAATAATCTCCTGGACAATCACACGGAATTTAGCCTGACggg                                                                                        >  2:651526/1‑150 (MQ=255)
                                                                tGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTa                                                                                      <  1:1071271/150‑1 (MQ=255)
                                                                    cATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                                                                                                                 <  1:1137117/103‑1 (MQ=255)
                                                                    cATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                                                                                                                 >  2:1137117/1‑103 (MQ=255)
                                                                      tCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTg                                                                                                                     <  1:562973/113‑1 (MQ=255)
                                                                      tCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTg                                                                                                                     >  2:562973/1‑113 (MQ=255)
                                                                        ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCtggtg                                                                              <  2:741288/150‑1 (MQ=255)
                                                                            acacCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTggcgg                                                                          <  2:230556/150‑1 (MQ=255)
                                                                             cacCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTggcggc                                                                         >  1:1541185/1‑150 (MQ=255)
                                                                                        gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCa                                                              <  2:691365/150‑1 (MQ=255)
                                                                                          tgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCaaa                                                            >  1:794886/1‑150 (MQ=255)
                                                                                          tgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCCCAGTCACCAGCCAGCCGCCAATAATCGCATGGACAATCAGACGGATATTATCCTGACCGGTACGCTCGTGGCGGCATTGATCGTCaaa                                                            >  2:574244/1‑150 (MQ=255)
                                                                                                cAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTg                                                      <  1:913223/150‑1 (MQ=255)
                                                                                                 aGAATAATGACTGACAAACCATTCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGc                                                     <  1:444836/150‑1 (MQ=255)
                                                                                                 aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGc                                                     >  2:682150/1‑150 (MQ=255)
                                                                                                       atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGAGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACt                                               >  2:474664/1‑150 (MQ=255)
                                                                                                       atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGAGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACt                                               >  2:469091/1‑150 (MQ=255)
                                                                                                        tgactgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACtt                                              >  2:283202/1‑150 (MQ=255)
                                                                                                            tgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTcgcg                                          >  1:78326/1‑150 (MQ=255)
                                                                                                             gacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTTGACAATCAGAAGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCgg                                         <  2:1545312/150‑1 (MQ=255)
                                                                                                             gacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCgg                                         >  1:511193/1‑150 (MQ=255)
                                                                                                             gacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCgg                                         >  2:1304152/1‑150 (MQ=255)
                                                                                                               cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGAc                                       >  1:937437/1‑150 (MQ=255)
                                                                                                               cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGAc                                       >  1:444334/1‑150 (MQ=255)
                                                                                                                aaaCCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACt                                      >  2:308396/1‑150 (MQ=255)
                                                                                                                 aaCCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACtt                                     >  2:809672/1‑150 (MQ=255)
                                                                                                                 aaCCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACtt                                     >  2:67404/1‑150 (MQ=255)
                                                                                                                  aCCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACttt                                    >  1:192364/1‑150 (MQ=255)
                                                                                                                    cAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTtct                                  <  1:488183/150‑1 (MQ=255)
                                                                                                                    cAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTtct                                  <  1:865761/150‑1 (MQ=255)
                                                                                                                     aaTCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTtctc                                 <  1:327695/150‑1 (MQ=255)
                                                                                                                       tCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCg                               <  2:78326/150‑1 (MQ=255)
                                                                                                                       tCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCg                               <  1:95991/150‑1 (MQ=255)
                                                                                                                         aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAGCTGTTGCAACACTTCGCGGACTTTCTCCGGc                             <  2:214794/150‑1 (MQ=255)
                                                                                                                         aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGc                             >  1:222709/1‑150 (MQ=255)
                                                                                                                            ccAACTTCAGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGCCGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGc                          >  1:1223462/1‑150 (MQ=255)
                                                                                                                             cAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCg                         <  1:1401414/150‑1 (MQ=255)
                                                                                                                             cAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCg                         >  1:1371389/1‑150 (MQ=255)
                                                                                                                                cTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGttt                      >  1:686607/1‑150 (MQ=255)
                                                                                                                                  tCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCa                    >  2:916539/1‑150 (MQ=255)
                                                                                                                                  tCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGCAGCCCTGCTGGCGCCTATGATCCTCAACCTGTTGCAACACTTAGCCGACTTTCTCCGGCAGCgaatca                    >  1:579231/1‑145 (MQ=255)
                                                                                                                                    cGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCAcc                  <  2:1541185/150‑1 (MQ=255)
                                                                                                                                    cGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCAcc                  >  2:861466/1‑150 (MQ=255)
                                                                                                                                    cGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCAcc                  <  2:616483/150‑1 (MQ=255)
                                                                                                                                          aTGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCACCACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGcc            >  1:1454232/1‑150 (MQ=255)
                                                                                                                                           tGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTTACCGTAGCCa           <  1:1304700/150‑1 (MQ=255)
                                                                                                                                            gCAACGCCAGAGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCaa          <  1:469091/150‑1 (MQ=255)
                                                                                                                                            gCAACGCCAGAGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCaa          <  1:474664/150‑1 (MQ=255)
                                                                                                                                              aaCGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAAc        <  2:937527/150‑1 (MQ=255)
                                                                                                                                                  ccAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAAc    >  1:1357816/1‑150 (MQ=255)
                                                                                                                                                  ccAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAAc    >  1:393076/1‑150 (MQ=255)
                                                                                                                                                   cAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAACg   <  1:223865/150‑1 (MQ=255)
                                                                                                                                                   cAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAACg   >  2:1402211/1‑150 (MQ=255)
                                                                                                                                                    aGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAACGc  <  2:937437/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GAGAAAAACACCGTCAACAGTGCGGTGAATGGCAGAGATTCGGTGAAGGCTTTACCGATAGCATGCTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCTGACGGGTACGCTGGTGGCGGCTTTGATCGTCAAACTGTTGCAACACTTCGCGGACTTTCTCCGGCAGCGTTTCACCGTAGCCAAACCAACGC  >  NC_000913/1233633‑1233929

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