Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,985,345 G→A 6.9% S20F (TCC→TTC)  PP_2611 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,985,3450GA 6.9% 171.0 / 2.4 58S20F (TCC→TTC) PP_2611membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base A (2/2);  total (26/32)
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.77e-01

TCTATATTTTCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGTAATTAATCGTTAGCGTGGC  >  NC_002947/2985198‑2985491
                                                                                                                                                   |                                                                                                                                                  
tCTATATTTTCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTAGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAAATGGTTGTTTGAATAAGTTGTCCgg                                                                                                                                                    >  5:92534/1‑148 (MQ=255)
tCTATATTTTCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATAAAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCgg                                                                                                                                                    >  4:227114/1‑148 (MQ=255)
      ttttCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAAtttt                                                                                                                                             <  3:214391/149‑1 (MQ=255)
       tttCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAAttttt                                                                                                                                            >  2:347415/1‑149 (MQ=255)
        ttCATTAATGAGATTGTAACATTAAAATACCCAGGTAAATATAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGTGTTGCGCAGGCATAAGTGTTAACTGGTTGTTTGAATAAGTTGTCCGGAAAttttt                                                                                                                                            <  7:346866/148‑1 (MQ=255)
        ttCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAAttttt                                                                                                                                            <  2:197179/148‑1 (MQ=255)
         tCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGCCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTg                                                                                                                                           <  6:96544/148‑1 (MQ=255)
            ttAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCtt                                                                                                                                       >  5:307761/1‑149 (MQ=255)
                           cATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGtgtg                                                                                                                        <  3:227114/149‑1 (MQ=255)
                               ataataCCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGAAAATTTTTGGCTTGAATCCTGAAGTGTGGGCt                                                                                                                    >  2:229553/1‑149 (MQ=255)
                               ataataCCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGAAAATTTTTGGCTTGAATCCTGAAGTGTGGGCt                                                                                                                    >  2:235784/1‑149 (MQ=255)
                                  ataCCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTg                                                                                                                  <  6:230079/148‑1 (MQ=255)
                                     cccAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGttt                                                                                                               <  1:347415/148‑1 (MQ=255)
                                      ccAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTa                                                                                                              >  3:270989/1‑148 (MQ=255)
                                            aaaTAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCtt                                                                                                                   >  6:123834/1‑137 (MQ=255)
                                            aaaTAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCtt                                                                                                                   <  5:123834/137‑1 (MQ=255)
                                             aaTAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTTAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGcc                                                                                                       <  6:53542/148‑1 (MQ=255)
                                               tAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGc                                                                                                        >  3:10838/1‑145 (MQ=255)
                                               tAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGc                                                                                                        >  1:176329/1‑145 (MQ=255)
                                                 gcaccATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTTTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATa                                                                                                  <  8:272493/145‑1 (MQ=255)
                                                 gAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGAt                                                                                                   <  8:334930/148‑1 (MQ=255)
                                                 gAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATa                                                                                                  <  5:20284/149‑1 (MQ=255)
                                                             tgtAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGc                                                                                             >  2:288525/1‑142 (MQ=255)
                                                             tgtAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGCCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGc                                                                                             <  1:288525/142‑1 (MQ=255)
                                                                 gATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGAAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAg                                                                                  <  1:229553/149‑1 (MQ=255)
                                                                 gATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGAAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAg                                                                                  <  1:235784/149‑1 (MQ=255)
                                                                  aTTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGttttct                                                                                 >  1:145052/1‑146 (MQ=255)
                                                                     tttCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGtt                                                                              >  5:256242/1‑149 (MQ=255)
                                                                       tCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTc                                                                             <  2:285691/148‑1 (MQ=255)
                                                                         ccATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTc                                                                             >  3:117926/1‑146 (MQ=255)
                                                                          cATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTc                                                                             >  5:181952/1‑145 (MQ=255)
                                                                           aTGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGAt                                                                        <  7:160619/149‑1 (MQ=255)
                                                                           aTGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGAt                                                                        >  5:210164/1‑149 (MQ=255)
                                                                                 aCAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAg                                                                  >  3:350591/1‑149 (MQ=255)
                                                                                 aCAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAg                                                                  >  2:52339/1‑149 (MQ=255)
                                                                                   agagTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTg                                                                >  2:246865/1‑149 (MQ=255)
                                                                                    gagTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTgg                                                                                                                       >  5:223748/1‑93 (MQ=255)
                                                                                    gagTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTgg                                                                                                                       <  6:223748/93‑1 (MQ=255)
                                                                                             acCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGt                                                                                                >  8:229548/1‑107 (MQ=255)
                                                                                             acCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAg                                                      <  2:176329/149‑1 (MQ=255)
                                                                                               ccAAGGGTTGCGGAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCTGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGa                                                    <  4:10838/149‑1 (MQ=255)
                                                                                               ccAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAg                                                      >  2:317653/1‑147 (MQ=255)
                                                                                                  aGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAgtg                                                 <  4:117926/149‑1 (MQ=255)
                                                                                                     gTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGaa                                                                     >  3:177792/1‑126 (MQ=255)
                                                                                                     gTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGaa                                                                     <  4:177792/126‑1 (MQ=255)
                                                                                                     gTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTgtcg                                              <  3:146340/149‑1 (MQ=255)
                                                                                                            aGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCa                                       <  1:52339/149‑1 (MQ=255)
                                                                                                             ggCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAg                                      <  4:270989/149‑1 (MQ=255)
                                                                                                                 tAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAATAGAAGGAGTGTCGTCATCCAGTCCt                                  <  6:210164/149‑1 (MQ=255)
                                                                                                                   aGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCttg                                <  1:317653/149‑1 (MQ=255)
                                                                                                                    gtgtGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCttgg                               <  6:307761/149‑1 (MQ=255)
                                                                                                                      gtgAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGTTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCttggtt                             >  1:23022/1‑149 (MQ=255)
                                                                                                                        gAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATc                                         >  2:340688/1‑135 (MQ=255)
                                                                                                                        gAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATc                                         <  1:340688/135‑1 (MQ=255)
                                                                                                                         aaCTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCttggttgg                           >  8:185900/1‑148 (MQ=255)
                                                                                                                           cTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCa                         >  3:121360/1‑148 (MQ=255)
                                                                                                                                  tttGAATAAGTTGTCCGGAAATTTTTGGCTTGATTCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGtaat                  >  5:93061/1‑148 (MQ=255)
                                                                                                                                  tttGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGtaatt                 <  6:92534/149‑1 (MQ=255)
                                                                                                                                                 cGGAAATTTTTGGTTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGTAATTAATCGTTAGCGTgg   <  2:23022/148‑1 (MQ=255)
                                                                                                                                                 cGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGTAATTAATCGTTAGCGTGGc  <  6:181952/149‑1 (MQ=255)
                                                                                                                                                  ggAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGGATAAGGAGTGTCGTCATACAGTCCTTGGTTGGCATAGtaa                   >  6:105117/1‑131 (MQ=255)
                                                                                                                                                   gAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGTAATTAATCGTTAGCGTgg   >  8:104962/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTATATTTTCATTAATGAGATTGTAACATGATAATACCCAGGTAAATAGAAACATGGCAGTGTAGATTTTTCCCATGAATACAGAGTCGGACACCCAAGGGTTGCGCAGGCATAAGTGTGAACTGGTTGTTTGAATAAGTTGTCCGGAAATTTTTGGCTTGAATCCTGAAGTGTGGGCTTGTTTAAAAAGCCAGATAGTCGCCATAGGTTTAGCGTTCCAGATGAATAGTGGTAGTAGAAGGAGTGTCGTCATCCAGTCCTTGGTTGGCATAGTAATTAATCGTTAGCGTGGC  >  NC_002947/2985198‑2985491

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