Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,966,389 A→T 5.2% intergenic (‑402/+28) PP_5232 ← / ← amtB hypothetical protein/ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,966,3890AT 5.2% 230.1 / 2.9 76intergenic (‑402/+28)PP_5232/amtBhypothetical protein/ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base T (2/2);  total (35/41)
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.11e-01

GCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGG  >  NC_002947/5966246‑5966538
                                                                                                                                               |                                                                                                                                                     
gcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                                 >  1:1669424/1‑150 (MQ=255)
gcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                                 >  2:1346878/1‑150 (MQ=255)
gcgcATTCTAGCGCGCCCGTGGAAAAAACAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                                 >  1:544258/1‑150 (MQ=255)
       cTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTACGGCTTGTAAGTTTTTTACCTACGTCACGCACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTCAAAACGTGGCAAAAAAACGCCCGGCAAGCCGGGCAAGttt                                                                                                                                          <  1:1610266/150‑1 (MQ=255)
        tAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCat                                                                                                                                              <  2:837051/145‑2 (MQ=255)
        tAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCat                                                                                                                                              >  1:837051/1‑144 (MQ=255)
             cgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACcg                                                                                                                                    <  1:431991/150‑1 (MQ=255)
                    ggAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCa                                                                                                                             <  2:1051973/150‑1 (MQ=255)
                         aagcaaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTAc                                                                                                                        <  1:1101324/150‑1 (MQ=255)
                          agcaaagcaCCATAAACGACCGGCTGCTTGTAAGTTTTGTACCCACGGCACGAACGTGTCGAATAAACCGTTTGCCGCGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGTAAGTTTTTACCGCGTCGCAGTTACa                                                                                                                       <  1:844797/150‑1 (MQ=255)
                           gcaaagcaCCATAAACGCCCTGCAACTTTTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAGACCGTTTGCCGTGCGGGAAAGGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAg                                                                                                                      <  1:1754802/150‑1 (MQ=255)
                            caaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAgg                                                                                                                     <  2:1364811/150‑1 (MQ=255)
                             aaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGt                                                                                                                    <  1:1403196/150‑1 (MQ=255)
                             aaagcaCCATAAACGCCCTGCGACTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATCAACCTTTTGCCGTGAGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGt                                                                                                                    <  1:437724/150‑1 (MQ=255)
                               agcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCATtttttt                                                                                                                                        <  1:539307/128‑1 (MQ=255)
                               agcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCATtttttt                                                                                                                                        >  2:539307/1‑128 (MQ=255)
                                    cATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                                             >  2:1150563/1‑150 (MQ=255)
                                     aTAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCcg                                                                                                            >  2:1217373/1‑150 (MQ=255)
                                      tAAACGCCCTGCGGCTTGTAAGTTTTTTTTTTTTTTTACGAACGTGTCGAATAAACCGTTTGTCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACTGCGTCGTATTTACTGTTTGTAGCcgc                                                                                                           >  2:829905/1‑150 (MQ=255)
                                      tAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCcgc                                                                                                           >  2:1332030/1‑150 (MQ=255)
                                        aaCGCCCTGCGGCTTGTAAGTTTTTAGCCGACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCgcgc                                                                                                         >  1:1717290/1‑150 (MQ=255)
                                          cGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTc                                                                                                       >  2:1205661/1‑150 (MQ=255)
                                             ccTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAACCCGCGAGCGGTAAATTCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                                    >  1:1304721/1‑150 (MQ=255)
                                             ccTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACTGGTTGTAGCCGCGCTCGtt                                                                                                    >  2:1279483/1‑150 (MQ=255)
                                                    cTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                                             <  1:834909/150‑1 (MQ=255)
                                                       gTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGcc                                                                                          <  2:1147294/150‑1 (MQ=255)
                                                         aaGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGAGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGATGTCCTAAACCAg                                                                                        >  1:91352/1‑150 (MQ=255)
                                                           gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                                      <  1:1332030/150‑1 (MQ=255)
                                                           gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                                      >  1:1420593/1‑150 (MQ=255)
                                                            ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTc                                                                                     <  2:1641989/150‑1 (MQ=255)
                                                              ttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGa                                                                                   <  1:11028/150‑1 (MQ=255)
                                                                ttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAgg                                                                                 <  2:544258/150‑1 (MQ=255)
                                                                  aCCTACGTCACGAACGTGTCGGATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGcc                                                                               >  2:1640436/1‑150 (MQ=255)
                                                                   ccTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgtg                                                                                                 >  2:743772/1‑131 (MQ=255)
                                                                   ccTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgtg                                                                                                 <  1:743772/131‑1 (MQ=255)
                                                                    cTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAGATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccga                                                                             >  1:1253561/1‑150 (MQ=255)
                                                                      aCGTCACGAACGTGTCAAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccgacc                                                                           <  2:165427/150‑1 (MQ=255)
                                                                          cACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                                       >  1:1623745/1‑150 (MQ=255)
                                                                           accaaCGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGTCGACCGActc                                                                      <  2:291377/147‑1 (MQ=255)
                                                                           aCGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGActc                                                                      <  2:384778/150‑1 (MQ=255)
                                                                             gAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtctt                                                                    >  2:1409805/1‑150 (MQ=255)
                                                                                    tCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttcttct                                                               >  1:929570/1‑148 (MQ=255)
                                                                                    tCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttcttct                                                               <  2:929570/148‑1 (MQ=255)
                                                                                      gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGt                                                           >  1:1082216/1‑150 (MQ=255)
                                                                                      gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGt                                                           <  1:1346878/150‑1 (MQ=255)
                                                                                       aaTAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGtt                                                          <  1:361513/150‑1 (MQ=255)
                                                                                        aTAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTg                                                         <  1:1409805/150‑1 (MQ=255)
                                                                                           ctcggTTTGCGGTGGGTGAAATGCTGAAAAACTGTCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCTCAGTTACTGGTTGTAGCCGCGCTCGTTGTGCTGAGCTAGGTCGAGGCCGACCGACTCGTCGTCTTCGTTGAcc                                                      <  2:569007/146‑1 (MQ=255)
                                                                                            aCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAccc                                                     <  2:1183535/150‑1 (MQ=255)
                                                                                              cGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAcccag                                                   >  1:569007/1‑148 (MQ=255)
                                                                                                tttGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                                             >  2:519468/1‑90 (MQ=255)
                                                                                                tttGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                                             <  1:519468/90‑1 (MQ=255)
                                                                                                       tGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATc                                          <  2:1551950/150‑1 (MQ=255)
                                                                                                       tGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATc                                          <  2:1623745/150‑1 (MQ=255)
                                                                                                        gCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccg                                                                              <  1:1695006/113‑1 (MQ=255)
                                                                                                        gCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccg                                                                              >  2:1695006/1‑113 (MQ=255)
                                                                                                           ggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcacca                                      >  2:1656422/1‑150 (MQ=255)
                                                                                                               aTGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATc                                  <  2:1304721/150‑1 (MQ=255)
                                                                                                                tGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAAACTCTCTGGTTGTGCTGAGCCTGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATcc                                 <  1:1205661/150‑1 (MQ=255)
                                                                                                                  cTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                               <  2:1669424/150‑1 (MQ=255)
                                                                                                                   tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                                             <  1:1863480/87‑1 (MQ=255)
                                                                                                                   tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                                             >  2:1863480/1‑87 (MQ=255)
                                                                                                                   tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttc                                                                   >  1:1701279/1‑113 (MQ=255)
                                                                                                                   tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttc                                                                   <  2:1701279/113‑1 (MQ=255)
                                                                                                                   tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGc                              >  1:1342392/1‑150 (MQ=255)
                                                                                                                     aaaaaCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCAc                            <  2:629471/150‑1 (MQ=255)
                                                                                                                         aCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTg                        <  2:1420593/150‑1 (MQ=255)
                                                                                                                             gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCAc                            >  1:1340414/1‑142 (MQ=255)
                                                                                                                             gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCAc                            <  2:1340414/142‑1 (MQ=255)
                                                                                                                              cAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGtt                                                          <  2:1151744/111‑1 (MQ=255)
                                                                                                                              cAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGtt                                                          >  1:1151744/1‑111 (MQ=255)
                                                                                                                               aaaaaacTGCCGGGCAAGCCGGGCAAGTGTTTACGGCGGCGCATTTACAGGGTCTAGCCTCGCTCGTTGTGCTGAGCAAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTCAGAATc                  >  2:321041/1‑150 (MQ=255)
                                                                                                                                         ccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGa        <  1:1150563/150‑1 (MQ=255)
                                                                                                                                          cGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGAc       >  1:1715278/1‑150 (MQ=255)
                                                                                                                                           ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTCTGCTGAGCCCGGGCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGCCAATCACGTAGGTGACg      >  1:1746974/1‑150 (MQ=255)
                                                                                                                                            gCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGCCGGGCTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGa     >  1:1789821/1‑150 (MQ=255)
                                                                                                                                               aGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATgg  <  1:829905/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGG  >  NC_002947/5966246‑5966538

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