Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,966,333 A→T 10.7% 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,3330AT10.7% 186.5 / 9.0 70intergenic (‑402/+28)PP_5232/amtBhypothetical protein/ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base T (3/5);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCAC  >  NC_002947/5966187‑5966468
                                                                                                                                                  |                                                                                                                                       
ggCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGc                                                                                                                                       >  6:281074/1‑149 (MQ=255)
 gcgcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGcc                                                                                                                                      <  2:221382/149‑1 (MQ=255)
   gcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCgg                                                                                                                                    <  8:284073/149‑1 (MQ=255)
   gcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCgg                                                                                                                                    >  7:328881/1‑149 (MQ=255)
    cgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGGCCGGCAAGCCggg                                                                                                                                   <  8:194953/149‑1 (MQ=255)
       aTTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAATAACTGGCAAAAAAATGCCCGGCTTGCCGGGCat                                                                                                                                <  1:230125/149‑2 (MQ=255)
             gcgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGt                                                                                                                              <  3:19188/145‑1 (MQ=255)
             gcgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGt                                                                                                                              >  4:19188/1‑145 (MQ=255)
              cgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTa                                                                                                                         <  6:108655/149‑1 (MQ=255)
                cgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTAcc                                                                                                                       <  8:328881/149‑1 (MQ=255)
                 ccccGTCCAAAAACCAAACCACCAAAAACCCCCAGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAACTTGCCCGGCTTGCCGGGCATTTTTTTGcca                                                                                                                      <  2:143160/148‑2 (MQ=37)
                     gTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTc                                                                                                                  <  5:134599/149‑1 (MQ=255)
                        gAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCTGGCAAGTTTTTACCGCGTCGCa                                                                                                               <  4:215026/149‑1 (MQ=255)
                        gAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCa                                                                                                               <  2:74951/149‑1 (MQ=255)
                        gAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCa                                                                                                               <  2:346116/149‑1 (MQ=255)
                          aaaaGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCAAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCgg                                                                                                                                    <  8:68968/126‑1 (MQ=255)
                          aaaaGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCAAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCgg                                                                                                                                    >  7:68968/1‑126 (MQ=255)
                           aaagcaatgcACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGTGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGGAAGCCGGGCAAGTTTTTACCGCGTCGCAGtt                                                                                                            <  6:112684/149‑1 (MQ=255)
                                        aTAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                               >  8:204956/1‑149 (MQ=255)
                                             cGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCt                                                                                          >  3:265269/1‑149 (MQ=255)
                                              gCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAACTTGCCCGGCTTGCCGGGCATttttttgc                                                                                                                        >  6:284134/1‑116 (MQ=38)
                                              gCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAACTTGCCCGGCTTGCCGGGCATttttttgc                                                                                                                        <  5:284134/118‑3 (MQ=38)
                                                ccTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGt                                                                                       >  5:72017/1‑149 (MQ=255)
                                                 cTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGTCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                      <  1:272420/149‑1 (MQ=255)
                                                    cGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtg                                                                                     <  5:281074/147‑1 (MQ=255)
                                                     ggCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCTTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCTATCGtttt                                                                                    <  2:137813/147‑3 (MQ=255)
                                                     ggCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGc                                                                                  <  7:204956/149‑1 (MQ=255)
                                                     ggCTTGTAAGTTTTTTACCTACGTCACGAAAGTGTCGAATAAACCGTATGCCGTGCGGGAAATGATGAAAAACTGGAAATAAAATGCCATGCAAGCAGGGCAAGTTTTTACCGAGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgt                                                                                    >  1:137813/1‑147 (MQ=255)
                                                          gTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                             >  3:231580/1‑149 (MQ=255)
                                                               ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                        >  6:213051/1‑149 (MQ=255)
                                                                tttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAgg                                                                                                       <  7:27019/117‑1 (MQ=255)
                                                                tttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAgg                                                                                                       >  8:27019/1‑117 (MQ=255)
                                                                   ttACCTACGTCACGAACGTGTCGAATAACCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGa                                                                     <  4:305268/148‑1 (MQ=255)
                                                                     aCCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAgg                                                                   <  4:265269/148‑1 (MQ=255)
                                                                       cTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccg                                                                <  1:254/149‑1 (MQ=255)
                                                                              aCGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAACTGCTGAAAAACTGGCAAAAACATGCCCGGCAAGCCGGTCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                         <  3:194779/149‑1 (MQ=255)
                                                                                         gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAg                                                                    >  6:149330/1‑127 (MQ=255)
                                                                                         gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAg                                                                    <  5:149330/127‑1 (MQ=255)
                                                                                         gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCttcttcttc                                               <  6:72017/148‑1 (MQ=255)
                                                                                         gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCg                                              <  4:231580/149‑1 (MQ=255)
                                                                                         gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCg                                              >  7:186994/1‑149 (MQ=255)
                                                                                            tAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                               <  7:24647/113‑1 (MQ=255)
                                                                                            tAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                               >  8:24647/1‑113 (MQ=255)
                                                                                                 cGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCCGGTCGAGGCCGCCCGACTCTTCTTCTTCGTTGACCCg                                      >  2:289801/1‑149 (MQ=255)
                                                                                                     tGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAg                                   <  5:213051/148‑1 (MQ=255)
                                                                                                        cGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                         >  1:43312/1‑123 (MQ=255)
                                                                                                        cGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                         <  2:43312/123‑1 (MQ=255)
                                                                                                              ggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttcttct                                                 <  8:330545/125‑1 (MQ=255)
                                                                                                              ggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttcttct                                                 >  7:330545/1‑125 (MQ=255)
                                                                                                                   tGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                         <  6:315003/112‑1 (MQ=255)
                                                                                                                   tGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                         >  5:315003/1‑112 (MQ=255)
                                                                                                                   tGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACAGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATc                    >  7:102178/1‑149 (MQ=255)
                                                                                                                      tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                 >  5:141357/1‑149 (MQ=255)
                                                                                                                       gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAgg                                  >  2:33583/1‑131 (MQ=255)
                                                                                                                       gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCCACCGACTCTTCTTCTTCGTTGACCCGCAgg                                  <  1:33583/131‑1 (MQ=255)
                                                                                                                        aaaaaCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATc                            <  2:30994/136‑1 (MQ=255)
                                                                                                                        aaaaaCTGGCAAAAAAATGCCCGGCAAGCCCGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATc                            >  1:30994/1‑136 (MQ=255)
                                                                                                                           aaCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCt            >  2:212110/1‑149 (MQ=255)
                                                                                                                              tGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCtt           >  7:126699/1‑147 (MQ=255)
                                                                                                                               ggtaaaaactTGCCCGGCTTGCCGGGCATTTTTTTGCCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcac                          <  7:169737/128‑1 (MQ=37)
                                                                                                                               ggtaaaaactTGCCCGGCTTGCCGGGCATTTTTTTGCCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcac                          >  8:169737/4‑131 (MQ=37)
                                                                                                                                gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGt                                             >  8:204157/1‑111 (MQ=255)
                                                                                                                                gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGt                                             <  7:204157/111‑1 (MQ=255)
                                                                                                                                   aaaacttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGc                                 <  4:279558/114‑1 (MQ=255)
                                                                                                                                   aaaacttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGc                                 >  3:279558/7‑120 (MQ=255)
                                                                                                                                     aaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCAc  <  8:102178/149‑1 (MQ=255)
                                                                                                                                          gCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTc                                                                       >  1:329656/1‑75 (MQ=255)
                                                                                                                                          gCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTc                                                                       <  2:329656/75‑1 (MQ=255)
                                                                                                                                           cccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCAt                             >  1:287290/1‑116 (MQ=255)
                                                                                                                                           cccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCAt                             <  2:287290/116‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
GGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCAC  >  NC_002947/5966187‑5966468

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