Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,640,095 T→G 6.2% L396W (TTG→TGG)  PP_4950 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,640,0950TG 6.2% 203.5 / 3.7 64L396W (TTG→TGG) PP_4950hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base G (2/2);  total (33/31)
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.49e-01

ATCAGCAACCTGTGGTACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGC  >  NC_002947/5639950‑5640242
                                                                                                                                                 |                                                                                                                                                   
aTCAGCAACCTGTGGTACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCt                                                                                                                                                  >  8:55918/1‑149 (MQ=255)
             ggTACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTg                                                                                                                                     >  7:40436/1‑149 (MQ=255)
              gTACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTg                                                                                                                                     <  3:217700/148‑1 (MQ=255)
                 ccTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTggag                                                                                                                                  <  3:308622/148‑1 (MQ=255)
                   tGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCAATTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                               >  1:310815/1‑149 (MQ=255)
                     aaCTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                               <  2:319246/147‑1 (MQ=255)
                     aaCTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCt                                                                                                                             >  4:252894/1‑149 (MQ=255)
                       ctactCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGa                                                                                                                           <  4:196916/149‑1 (MQ=255)
                          ctCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGAccc                                                                                                                        <  6:266204/149‑1 (MQ=255)
                          ctCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGAcc                                                                                                                         >  6:324537/1‑148 (MQ=255)
                          ctCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGAcc                                                                                                                         >  6:180066/1‑148 (MQ=255)
                           tCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTgg                                                                                                                                                >  2:140451/1‑124 (MQ=255)
                           tCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTgg                                                                                                                                                <  1:140451/124‑1 (MQ=255)
                             aGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAg                                                                                                                      >  1:233311/1‑148 (MQ=255)
                                ccTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAgc                                                                                                                  >  8:146322/1‑149 (MQ=255)
                                    ccAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAgcgc                                                                                                                >  3:334479/1‑147 (MQ=255)
                                           cgcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                       <  5:147307/149‑1 (MQ=255)
                                               cATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGt                                                                                                   <  2:233311/149‑1 (MQ=255)
                                                  gACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                       >  8:98418/1‑142 (MQ=255)
                                                  gACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                       <  7:98418/142‑1 (MQ=255)
                                                   aCCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCgaga                                                                                                            >  2:82744/1‑136 (MQ=255)
                                                   aCCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGAc                                                                                               <  3:252894/149‑1 (MQ=255)
                                                   aCCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGAc                                                                                               <  7:146322/149‑1 (MQ=255)
                                                   aCCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACACCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCgaga                                                                                                            <  1:82744/136‑1 (MQ=255)
                                                       ggCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAg                                                                                            >  2:1204/1‑148 (MQ=255)
                                                              ccGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCTTGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGCGCTCGAGATGATCCTGTCGACCAGTACTTATg                                                                                    >  4:267048/1‑149 (MQ=255)
                                                                gCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                               >  2:367032/1‑104 (MQ=255)
                                                                gCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                               <  1:367032/104‑1 (MQ=255)
                                                                  ttCGCCACTTTCTGGGTAGAAGACGGGAAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACa                                                                                <  3:29736/149‑1 (MQ=255)
                                                                     tccACTTTCAGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGAGCTTCGATGACAGCTAGTACAGCTTGATTCGCAGCGAGCTGGAGGACCTGACCCAGGAGCGCGAGAAGATCCTGTCGAGCAGTACATATGGACcgcg                                                                             >  8:154311/2‑149 (MQ=255)
                                                                     gCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                       >  8:241710/1‑123 (MQ=255)
                                                                     gCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                       <  7:241710/123‑1 (MQ=255)
                                                                     gCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAgca                                                                             >  8:187652/1‑148 (MQ=255)
                                                                      ccACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAgcg                                                                             >  3:357042/1‑148 (MQ=255)
                                                                            tCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTAcc                                                                      >  8:42451/1‑149 (MQ=255)
                                                                             cTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTAcc                                                                      <  2:310815/148‑1 (MQ=255)
                                                                                ggTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGt                                                                  >  7:182790/1‑149 (MQ=255)
                                                                                    gaagaCGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGa                                                               <  1:1204/148‑1 (MQ=255)
                                                                                          gggCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCAttt                                                        >  1:63493/1‑149 (MQ=255)
                                                                                                  ccAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCgggg                                                >  3:15546/1‑149 (MQ=255)
                                                                                                  ccAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCggg                                                 <  2:63493/148‑1 (MQ=255)
                                                                                                   cAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCgggg                                                <  5:180066/148‑1 (MQ=255)
                                                                                                    aGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCa                                              >  5:76425/1‑149 (MQ=255)
                                                                                                               aGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAgg                                   >  2:328125/1‑149 (MQ=255)
                                                                                                               aGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAg                                    <  8:40436/148‑1 (MQ=255)
                                                                                                               aGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAg                                    <  4:357042/148‑1 (MQ=255)
                                                                                                                gTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGc                                  <  4:334479/149‑1 (MQ=255)
                                                                                                                   ccATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcct                                >  2:57963/1‑148 (MQ=255)
                                                                                                                      tGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgacc                            <  7:42451/149‑1 (MQ=255)
                                                                                                                            tCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCgg                                                                    <  2:153001/103‑1 (MQ=255)
                                                                                                                            tCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCgg                                                                    >  1:153001/1‑103 (MQ=255)
                                                                                                                               aTGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATtg                   <  1:328125/149‑1 (MQ=255)
                                                                                                                                tGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATtgt                  <  6:76425/149‑1 (MQ=255)
                                                                                                                                  aCAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCTTTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGa                >  7:72911/1‑149 (MQ=255)
                                                                                                                                          tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgaccct                          >  7:317353/3‑131 (MQ=255)
                                                                                                                                          tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgaccct                          <  8:137182/129‑1 (MQ=255)
                                                                                                                                          tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgaccct                          >  7:137182/3‑131 (MQ=255)
                                                                                                                                          tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgaccct                          <  8:317353/129‑1 (MQ=255)
                                                                                                                                            cAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGc      >  7:16427/1‑149 (MQ=255)
                                                                                                                                             aGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCt                                                                                                     <  6:79518/53‑1 (MQ=255)
                                                                                                                                             aGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCt                                                                                                     >  5:79518/1‑53 (MQ=255)
                                                                                                                                             aGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGcc     >  7:362547/1‑149 (MQ=255)
                                                                                                                                               cTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCgg   <  8:72911/149‑1 (MQ=255)
                                                                                                                                               cTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGAAATTGTGACTGGACTTGCCg    <  8:362547/148‑1 (MQ=255)
                                                                                                                                                ttGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAg                                    >  5:97915/1‑115 (MQ=255)
                                                                                                                                                ttGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAg                                    <  6:97915/115‑1 (MQ=255)
                                                                                                                                                ttGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGc  >  5:218853/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCAGCAACCTGTGGTACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGC  >  NC_002947/5639950‑5640242

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: ‑