Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,094,128 G→T 6.7% S13S (TCC→TCA ptsN ← phosphotransferase system subunit IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,094,1280GT 6.7% 240.8 / 6.1 75S13S (TCC→TCAptsNphosphotransferase system subunit IIA
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base T (2/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

ATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTGCAT  >  NC_002947/1093981‑1094276
                                                                                                                                                         |                                                                                                                                                       
aTTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAgg                                                                                                                                                         >  4:380824/1‑148 (MQ=255)
 ttGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAggg                                                                                                                                                        >  7:121750/1‑148 (MQ=255)
  tGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGaa                                                                                                                                                      >  5:60460/1‑149 (MQ=255)
       aaaCCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGcc                                                                                                                                                <  8:313490/149‑1 (MQ=255)
           cGGTATAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCCGTTCCCTCGTGAACGTGCCggga                                                                                                                                            <  8:362220/149‑2 (MQ=255)
                gAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAg                                                                                                                                       <  6:60460/149‑1 (MQ=255)
                   ccGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGAt                                                                                                                                    >  4:233614/1‑149 (MQ=255)
                    cGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATg                                                                                                                                   <  1:331111/149‑1 (MQ=255)
                    cGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATg                                                                                                                                   >  3:330163/1‑149 (MQ=255)
                    cGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCCGTTCCCTCGTGAACGTGCCggg                                                                                                                                             >  4:281506/1‑139 (MQ=255)
                    cGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCCGTTCCCTCGTGAACGTGCCggg                                                                                                                                             <  3:281506/139‑1 (MQ=255)
                           ttCCCCGGCTACCAGCTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCTCTTGTCGCCGTCGGCGATTACCCCTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGAt                                                                                                                                    <  1:125922/147‑1 (MQ=255)
                            tCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCaa                                                                                                                            <  3:275227/148‑1 (MQ=255)
                             cccGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGCAACGGCCGGGGGTCAGGATGGTTTCAAg                                                                                                                           <  1:227657/148‑1 (MQ=255)
                              ccGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTc                                                                                                                         <  3:380824/149‑1 (MQ=255)
                               cGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTc                                                                                                                              >  1:366302/1‑143 (MQ=255)
                               cGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTc                                                                                                                              <  2:366302/143‑1 (MQ=255)
                                    tACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCa                                                                                                                   >  1:240326/1‑149 (MQ=255)
                                    tACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCa                                                                                                                   >  3:169659/1‑149 (MQ=255)
                                          tttttCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATg                                                                                                             <  2:189986/149‑1 (MQ=255)
                                              tctaattcGTCCTGGATCTCCAGTTCTGGCTCCTG‑‑‑‑‑GTCGGAGTTGACAGTGCCTACCTTATGCTGGGCTCGCTTGTGGCTGCATCCCTTCACGTTCACGAGGGAACGGCGGGGGGGCTGTAAGGTTTCAAGTCCGATGATGGATCAgg                                                                                                         <  2:69320/141‑1 (MQ=255)
                                               cGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAggg                                                                                                       <  4:45838/148‑1 (MQ=255)
                                                gAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAgggg                                                                                                      >  8:27152/1‑148 (MQ=255)
                                                gAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAgggg                                                                                                      >  4:40781/1‑148 (MQ=255)
                                                gAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAgggg                                                                                                      >  5:148390/1‑148 (MQ=255)
                                                gAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCCCCCGGCCGTTCCCTCGTGAACGTGCCggga                                                                                                                                            >  4:151440/1‑111 (MQ=37)
                                                gAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCCCCCGGCCGTTCCCTCGTGAACGTGCCggga                                                                                                                                            <  3:151440/112‑2 (MQ=37)
                                                   gACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGAt                                                                                                   >  3:309328/1‑148 (MQ=255)
                                                   gACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGAt                                                                                                   >  2:120995/1‑148 (MQ=255)
                                                      gTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGc                                                                                               >  1:246156/1‑149 (MQ=255)
                                                       tCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCg                                                                                              >  1:6322/1‑149 (MQ=255)
                                                       tCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCg                                                                                              >  1:269458/1‑149 (MQ=255)
                                                           gCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTTACGAGGGAACGGCCGGGGGTCAGGATGGTTTCACGTCGGATCATGGATGAGGGGGATCAGCggg                                                                                            >  6:282641/1‑147 (MQ=255)
                                                             aTCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGc                                                                                        >  5:211379/1‑149 (MQ=255)
                                                                 ccAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCa                                                                                    >  4:240019/1‑149 (MQ=255)
                                                                  cAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCa                                                                                    <  8:186906/148‑1 (MQ=255)
                                                                       cTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGACGGGGGATCAGCGGGCAGCTGCACCtt                                                                                <  1:210876/148‑1 (MQ=255)
                                                                          gCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTgcagc                                                                           <  4:309328/149‑1 (MQ=255)
                                                                          gCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTgcag                                                                            <  2:385/148‑1 (MQ=255)
                                                                            aCTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGcagcag                                                                         >  6:142728/1‑149 (MQ=255)
                                                                              ttG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAg                                                                                         >  4:129771/1‑131 (MQ=255)
                                                                              ttG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAg                                                                                         <  3:129771/131‑1 (MQ=255)
                                                                                       tCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTg                                                                   >  6:381760/1‑149 (MQ=255)
                                                                                       tCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTg                                                                   <  2:269458/149‑1 (MQ=255)
                                                                                               gACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAg                                                                                         >  6:80726/1‑119 (MQ=255)
                                                                                               gACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAg                                                                                         <  5:80726/119‑1 (MQ=255)
                                                                                                aCAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTc                                                          <  3:240019/149‑1 (MQ=255)
                                                                                                    gTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCtt                                                       <  4:330163/149‑1 (MQ=255)
                                                                                                    gTGGCGACCGTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTTACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGTATCATTGTTGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCtt                                                       <  5:282641/149‑1 (MQ=255)
                                                                                                          aCCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGttt                                                 >  4:37815/1‑149 (MQ=255)
                                                                                                           ccTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGtttt                                                <  2:240326/149‑1 (MQ=255)
                                                                                                                tCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCa                                            <  3:40781/148‑1 (MQ=255)
                                                                                                                  cAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGtt                                         <  3:37815/149‑1 (MQ=255)
                                                                                                                   aGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCAGTAGGTAACGGCGGGGGGTAAGTATGGTTTCAAGTCGTATCTTGGATTTGGACGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGtt                                         <  1:120995/149‑1 (MQ=255)
                                                                                                                   aGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGtttt                                                            <  4:46740/129‑1 (MQ=255)
                                                                                                                   aGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGtttt                                                            >  3:46740/1‑129 (MQ=255)
                                                                                                                   aGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGtt                                         <  6:314365/148‑1 (MQ=255)
                                                                                                                   aGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGtt                                         >  5:314365/1‑148 (MQ=255)
                                                                                                                    ggtcgcgcTGCTTGCTGCCTCCCGGCACGTGCCCGAGGGAACGGCCGGGGGGCAGGAGGGGTTCAAGTCGGATCATGGCTGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTgg                                       <  4:210505/146‑1 (MQ=255)
                                                                                                                          gcTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTc                                  >  2:84884/1‑148 (MQ=255)
                                                                                                                              cttGCTGCCTCCCGGCACGTTCACGATGGAACGGCCGTGGGTCAGGATGGTTTCAAGTCGGATCATGTATGAGGGGGATCAGCGGGCATCTGCACCTTGCAGCAGTCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTctatc                              <  4:169659/148‑5 (MQ=255)
                                                                                                                                 gctgcCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTgcagc                                                                           >  2:363671/1‑100 (MQ=255)
                                                                                                                                 gctgcCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTgcagc                                                                           <  1:363671/100‑1 (MQ=255)
                                                                                                                                     ccTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCg                        <  5:157864/147‑1 (MQ=255)
                                                                                                                                     ccTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGc                      <  2:246156/149‑1 (MQ=255)
                                                                                                                                       tCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGcc                     >  6:324488/1‑148 (MQ=255)
                                                                                                                                       tCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGcc                     >  3:357362/1‑148 (MQ=255)
                                                                                                                                         ccGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAg                  >  1:160178/1‑149 (MQ=255)
                                                                                                                                         ccGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAg                  <  6:148390/149‑1 (MQ=255)
                                                                                                                                          cGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAgg                 >  8:197832/1‑149 (MQ=255)
                                                                                                                                            gCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCg               >  1:327648/1‑149 (MQ=255)
                                                                                                                                             cACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGt              <  8:121750/149‑1 (MQ=255)
                                                                                                                                                    aCGAGGGAACGGCCGGGGGACAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGc       >  3:281303/1‑149 (MQ=255)
                                                                                                                                                     cGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCt      <  5:142544/149‑1 (MQ=255)
                                                                                                                                                      gAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTg     >  8:40066/1‑149 (MQ=255)
                                                                                                                                                      gAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTg     >  5:86011/1‑149 (MQ=255)
                                                                                                                                                       aGGGAACGGCCGGGGGTCAGGATGGTCTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTGc    <  3:233614/149‑1 (MQ=255)
                                                                                                                                                         ggAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTGCAt  <  5:324488/149‑1 (MQ=255)
                                                                                                                                                         ggAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGACGGCCAAGGCGTCGATCGCTGCa   >  2:38197/1‑148 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ATTGCCGAAACCGGTAGAACCGAGCTTTTCCCGGGCTACCAGTTTTTCGAAGACGTCCTGCATCTCCAGTTCTGGCACTTG‑‑‑‑‑GTCGGCGATGACAGTGGCGACCTTTTCCAGGGCGCGCTTCTTGCTGCCTCCCGGCACGTTCACGAGGGAACGGCCGGGGGTCAGGATGGTTTCAAGTCGGATCATGGATGAGGGGGATCAGCGGGCAGCTGCACCTTGCAGCAGGCTTTGCTGTTTTTCCTTGTGTTTTTTCAGTTGGCGGTCGAGCTTGTCGGCCAAGGCGTCGATCGCTGCAT  >  NC_002947/1093981‑1094276

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: