Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,526,256 C→T 5.8% G201S (GGC→AGC) ‡ yfeR ← putative LysR‑type DNA‑binding transcriptional regulator YfeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,526,2560CT 5.8% 440.5 / 10.0 138G201S (GGC→AGC) ‡yfeRputative LysR‑type DNA‑binding transcriptional regulator YfeR
Reads supporting (aligned to +/- strand):  ref base C (63/66);  new base T (6/2);  total (70/68)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 146 total aligned reads displayed.

AGGCAATGGTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTG  >  NC_000913/2526107‑2526405
                                                                                                                                                     |                                                                                                                                                     
aGGCAATGGTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGcc                                                                                                                                                       <  1:2289666/150‑1 (MQ=255)
aGGCAATGGTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGcc                                                                                                                                                       <  2:1028422/150‑1 (MQ=255)
 ggCAATGGTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCt                                                                                                                                                      <  2:427394/150‑1 (MQ=255)
       ggTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGAcgcg                                                                                                                                                >  1:1989333/1‑150 (MQ=255)
       ggTAGCGCCAGGGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGAcgcg                                                                                                                                                <  1:1645916/150‑1 (MQ=255)
        gTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGt                                                                                                                                               <  2:1484391/150‑1 (MQ=255)
        gTAGCGCCAGTACGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAGGCAGCCGACCGCTGATTGAcgcag                                                                                                                                               <  1:187413/150‑3 (MQ=255)
             gcCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATcc                                                                                                                                          <  2:1122031/150‑1 (MQ=255)
             gcCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATcc                                                                                                                                          <  1:1019525/150‑1 (MQ=255)
              ccAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCt                                                                                                                                         >  1:221140/1‑150 (MQ=255)
                aGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGc                                                                                                                                       >  2:1689457/1‑150 (MQ=255)
                aGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGc                                                                                                                                       <  1:605598/150‑1 (MQ=255)
                 gTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCa                                                                                                                                      <  2:2195723/150‑1 (MQ=255)
                    cGGGAAGTATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAaca                                                                                                                                   >  1:2083945/1‑150 (MQ=255)
                     gggAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAacac                                                                                                                                  <  1:815893/150‑1 (MQ=255)
                      ggAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACAcc                                                                                                                                 >  1:1088358/1‑150 (MQ=255)
                       gaagaaTACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCa                                                                                                                                <  2:1925076/150‑1 (MQ=255)
                            aTACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTc                                                                                                                           <  1:2119825/150‑1 (MQ=255)
                              aCTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGc                                                                                                                         >  2:639928/1‑150 (MQ=255)
                              aCTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGc                                                                                                                         <  2:987912/150‑1 (MQ=255)
                              aCTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGc                                                                                                                         >  1:1918058/1‑150 (MQ=255)
                               cTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCt                                                                                                                        >  1:593963/1‑150 (MQ=255)
                                tGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTc                                                                                                                       <  2:1809126/150‑1 (MQ=255)
                                 gatgccgatgccTGCGGCTACCATCGGGAACAGCGTCGACGTATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGAGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTcc                                                                                                                      >  2:1653743/1‑150 (MQ=255)
                                  atgccgatgccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTccc                                                                                                                     <  1:1604028/150‑1 (MQ=255)
                                    gccgatgccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTg                                                                                                                   <  2:232821/150‑1 (MQ=255)
                                    gccgatgccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTg                                                                                                                   <  1:2387598/150‑1 (MQ=255)
                                    gccgatgccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTg                                                                                                                   <  1:1425987/150‑1 (MQ=255)
                                      cgatgccTGCGGCTACCATCGGGAACAGCGGCGCCGGATGTCCAATCTCCTGCACAAAATTCGCCTGAATACCATTGCGCGCCATCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACGAATTTCGCTCCCTGTa                                                                                                                 <  2:2336660/150‑1 (MQ=255)
                                          gccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAAtgc                                                                                                             >  1:1211882/1‑150 (MQ=255)
                                           ccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAAtgcc                                                                                                            >  1:2104012/1‑150 (MQ=255)
                                           ccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAAtgcc                                                                                                            >  2:1328405/1‑150 (MQ=255)
                                           ccTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAAtgcc                                                                                                            <  1:897156/150‑1 (MQ=255)
                                                 gCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAg                                                                                                      >  1:127975/1‑150 (MQ=255)
                                                   tACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAggg                                                                                                    <  2:2147178/150‑1 (MQ=255)
                                                    aCCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGt                                                                                                   >  1:2006132/1‑150 (MQ=255)
                                                      cATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTAc                                                                                                 <  1:334285/150‑1 (MQ=255)
                                                       aTCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACa                                                                                                <  2:1726526/150‑1 (MQ=255)
                                                       aTCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACa                                                                                                >  2:406549/1‑150 (MQ=255)
                                                        tCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACAt                                                                                               <  1:1198552/150‑1 (MQ=255)
                                                         cGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATa                                                                                              >  2:2214757/1‑150 (MQ=255)
                                                          gggAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATaa                                                                                             <  1:61943/150‑1 (MQ=255)
                                                          gggAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATaa                                                                                             <  2:2479503/150‑1 (MQ=255)
                                                           ggAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCa                                                                                                       >  2:1604719/1‑139 (MQ=255)
                                                           ggAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAAt                                                                                            >  1:1557449/1‑150 (MQ=255)
                                                               cAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATATTc                                                                                        <  1:732987/150‑1 (MQ=255)
                                                                aGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTcc                                                                                       <  2:1088358/150‑1 (MQ=255)
                                                                     cGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGc                                                                                  <  2:1563788/150‑1 (MQ=255)
                                                                      gCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGcc                                                                                 <  2:238588/150‑1 (MQ=255)
                                                                             gTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAgcgc                                                                          >  2:1031083/1‑150 (MQ=255)
                                                                             gTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAgcgc                                                                          >  1:1257418/1‑150 (MQ=255)
                                                                              tCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCt                                                                         >  2:1775088/1‑150 (MQ=255)
                                                                              tCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCt                                                                         >  2:787674/1‑150 (MQ=255)
                                                                              tCCAATCTCCTGCACAAAATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACCCCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCATGGCCAAAGCGCt                                                                         >  2:2004692/1‑150 (MQ=255)
                                                                               ccAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTa                                                                        >  1:2323658/1‑150 (MQ=255)
                                                                                 aaTCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATc                                                                      >  2:273693/1‑150 (MQ=255)
                                                                                       cTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATg                                                                <  2:2052781/150‑1 (MQ=255)
                                                                                        tGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATgg                                                               <  2:354907/150‑1 (MQ=255)
                                                                                          cacaATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCa                                                             <  2:96413/150‑1 (MQ=255)
                                                                                           acaATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCaa                                                            >  2:1737527/1‑150 (MQ=255)
                                                                                            caATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCaaa                                                           >  1:1461416/1‑150 (MQ=255)
                                                                                            caATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCa                                                             >  2:345268/1‑148 (MQ=255)
                                                                                              atatTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGc                                                         <  1:702348/150‑1 (MQ=255)
                                                                                              atatTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGc                                                         <  2:178589/150‑1 (MQ=255)
                                                                                                atTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCag                                                       >  1:558634/1‑150 (MQ=255)
                                                                                                   cGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCagaaa                                                    >  1:306206/1‑150 (MQ=255)
                                                                                                     ccTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCagaaaga                                                  >  2:691517/1‑150 (MQ=255)
                                                                                                      cTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCagaaagaa                                                 >  1:1103065/1‑150 (MQ=255)
                                                                                                       tGAATACCATTGCACGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAgaaagaaa                                                <  2:786403/150‑1 (MQ=255)
                                                                                                        gAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGaaagaaag                                               >  2:393765/1‑150 (MQ=255)
                                                                                                         aaTACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAgg                                              <  2:1258544/150‑1 (MQ=255)
                                                                                                             ccATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCg                                          >  1:854830/1‑150 (MQ=255)
                                                                                                             ccATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCg                                          >  1:90306/1‑150 (MQ=255)
                                                                                                               aTTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGa                                        <  2:593963/150‑1 (MQ=255)
                                                                                                                  gcgcgcCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAg                                     >  1:2178301/1‑150 (MQ=255)
                                                                                                                     cgcCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAAt                                  >  1:2432845/1‑150 (MQ=255)
                                                                                                                     cgcCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAAt                                  >  2:2365320/1‑150 (MQ=255)
                                                                                                                      gcCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATc                                 <  2:892837/150‑1 (MQ=255)
                                                                                                                           ttacgCGTCAGGCAGCCGACCGCTGATTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTc                            >  2:553062/5‑150 (MQ=255)
                                                                                                                            gcTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTca                           >  2:436532/1‑150 (MQ=255)
                                                                                                                             cTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACAAAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTcac                          <  2:2101399/150‑1 (MQ=255)
                                                                                                                              tGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTcaca                         >  2:472260/1‑150 (MQ=255)
                                                                                                                                cGTCAGGCAGCCGACCGCTGATTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACAtt                       >  1:683988/1‑150 (MQ=255)
                                                                                                                                cGTCAGGCAGCCGACCGCTGATTGACGAGTAATCCTGCAACACGAATTTCGCACCCTGTTATGCCTGGCAGGGTACATATTCTTCCACGGACAAAGCGCTATCGCGATGGCAAAGCAGAAAGCAAGGCTCGGAAAGAATCGCTTCACAtt                       >  2:1142696/1‑150 (MQ=255)
                                                                                                                                  tCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGc                     >  1:1151696/1‑150 (MQ=255)
                                                                                                                                    aaTCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAg                   >  1:1799650/1‑150 (MQ=255)
                                                                                                                                     aTCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGa                  >  1:1336967/1‑150 (MQ=255)
                                                                                                                                       cAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATc                <  1:2216393/150‑1 (MQ=255)
                                                                                                                                         gCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGc              <  1:406549/150‑1 (MQ=255)
                                                                                                                                            gTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCa           >  1:1170353/1‑150 (MQ=255)
                                                                                                                                                gCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTgg       <  1:2115393/150‑1 (MQ=255)
                                                                                                                                                  tGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGc     <  2:90306/150‑1 (MQ=255)
                                                                                                                                                  tGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGc     <  1:366981/150‑1 (MQ=255)
                                                                                                                                                   gCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGcc    >  2:1194347/1‑150 (MQ=255)
                                                                                                                                                   gCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGcc    <  2:2104012/150‑1 (MQ=255)
                                                                                                                                                    ccTGCCGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCt   >  2:1804581/1‑150 (MQ=255)
                                                                                                                                                     cTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTg  <  1:1263043/150‑1 (MQ=255)
                                                                                                                                                     cTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTg  <  2:1557449/150‑1 (MQ=255)
                                                                                                                                                     cTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTg  <  1:1254523/150‑1 (MQ=255)
                                                                                                                                                     cTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTg  <  2:1846587/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCAATGGTAGCGCCAGTGCGGGAAGAATACTGATGCCGATGCCTGCGGCTACCATCGGGAACAGCGTCGCCGGATGTCCAATCTCCTGCACAATATTCGCCTGAATACCATTGCGCGCCAGCGCTGCGTCAATCAGCGGTCGGCTGCCTGACGCGTAATCCTGCAACACCAATTTCGCTCCCTGTAATGCCTGCCAGGGTACATAATCTTCCACGGCCAAAGCGCTATCGCGATGGCAAAGCAGAAAGAAAGGCTCGGAAAGAATCGCTTCACATTGCAGATCGCCCACTGGGCCTG  >  NC_000913/2526107‑2526405

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: