Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,167,904 A→T 5.9% P328P (CCA→CCT) ‡ PP_5407 → TnsD‑like transposition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,167,9040AT 5.9% 202.8 / 2.7 68P328P (CCA→CCT) ‡PP_5407TnsD‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base T (2/2);  total (38/30)
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.89e-01

AAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167762‑6168050
                                                                                                                                               |                                                                                                                                                  
aaaTGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAgg                                                                                                                                              >  6:440145/1‑149 (MQ=255)
   tGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCGCTGGGAGGCTa                                                                                                                                           <  5:428417/149‑1 (MQ=255)
     aCCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCGCTGGGAGGCTag                                                                                                                                          <  6:314123/148‑2 (MQ=255)
             gCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaact                                                                                                                                  >  8:241206/1‑148 (MQ=255)
                ccGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTtaactaac                                                                                                                               >  3:373022/1‑148 (MQ=255)
                  cggggtACTGCCATCCGCTGACGCACCTGCGTCTGCTATGCTGGCTCCTCGGCCCCTTACAATCATCCGTTCACGCCCATCAACCTCCCGCGACCCAGCAGCAACTTCCAGTGACTACCCTCCCAGCGAGGCTAAACTTAACTAACAct                                                                                                                            <  4:358587/144‑1 (MQ=255)
                     gTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActtt                                                                                                                          <  5:401040/148‑1 (MQ=255)
                        aCTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActttctt                                                                                                                       <  1:294039/148‑1 (MQ=255)
                        aCTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActttctt                                                                                                                       <  5:440145/148‑1 (MQ=255)
                          tGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAaca                                                                                                                              <  1:434840/139‑1 (MQ=255)
                          tGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAaca                                                                                                                              >  2:434840/1‑139 (MQ=255)
                             cATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGc                                                                                                                  <  7:219799/148‑1 (MQ=255)
                              aTCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGcc                                                                                                                 >  3:38145/1‑148 (MQ=255)
                                  gTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTAttt                                                                                                            >  6:182842/1‑149 (MQ=255)
                                  gTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTAtt                                                                                                             >  7:347398/1‑148 (MQ=255)
                                   ttGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTAttt                                                                                                            >  8:385105/1‑148 (MQ=255)
                                      aaGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                        <  4:337199/149‑1 (MQ=255)
                                       aGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGcc                                                                                                       <  2:232239/149‑1 (MQ=255)
                                        gCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCTGACATGCCGCGAACCAGCAGCAAATTCCAGTGAATAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGcc                                                                                                       >  1:125829/1‑149 (MQ=255)
                                              gAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaagg                                                                                                  >  8:75376/1‑147 (MQ=255)
                                              gAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaagg                                                                                                  >  8:135612/1‑147 (MQ=255)
                                               aGTCTGATTTGGTTGCTATTCGGTAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCACCAACTTCCAGGGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaaggaag                                                                                               <  6:421925/149‑1 (MQ=255)
                                               aGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaaggaag                                                                                               <  2:307015/149‑1 (MQ=255)
                                               aGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaaggaa                                                                                                <  7:385105/148‑1 (MQ=255)
                                                 tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCCTTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAaggaagga                                                                                             <  4:296623/149‑1 (MQ=255)
                                                      tttGGTGGCTATTCGTGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACt                                                                                        <  4:333174/149‑1 (MQ=255)
                                                      tttGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACt                                                                                        >  3:57449/1‑149 (MQ=255)
                                                      tttGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGAGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACt                                                                                        >  2:330498/1‑149 (MQ=255)
                                                             gCTATTCGGGAGCTTATAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTg                                                                                    <  3:113464/146‑1 (MQ=255)
                                                             gCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                 <  8:225949/149‑1 (MQ=255)
                                                             gCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                 >  2:20354/1‑149 (MQ=255)
                                                             acTATTCTCGAGCTTAAAATCATTCTTTAAGAACTATCAAAATTCCTCGAACAATCACCAACTTCCATTCACTAGACTAACACCGAGGATAAACTTAAATAAAAATTACTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTg                                                                                    >  4:113464/2‑146 (MQ=255)
                                                               tATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCg                                                                               >  8:121533/1‑149 (MQ=255)
                                                                 ttCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaa                                                                             <  4:373022/149‑1 (MQ=255)
                                                                   cGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            <  8:347398/148‑1 (MQ=255)
                                                                    gggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            <  6:382144/147‑1 (MQ=255)
                                                                    gggAGCTTAGAATCATTCGTTAAGGACTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            >  5:382144/1‑147 (MQ=255)
                                                                                tCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAtt                                                              >  1:86117/1‑149 (MQ=255)
                                                                                  aTTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAttt                                                             <  2:125829/148‑1 (MQ=255)
                                                                                   ttCGTTAAGGCCTATCAACATGCCTCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                           <  7:135612/149‑1 (MQ=255)
                                                                                   ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCGAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                           >  4:31621/1‑149 (MQ=255)
                                                                                      gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTta                                                                                                                                     <  8:365618/72‑1 (MQ=255)
                                                                                      gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTta                                                                                                                                     >  7:365618/1‑72 (MQ=255)
                                                                                      gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATg                                                         <  1:20354/148‑1 (MQ=255)
                                                                                               tATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                 >  1:301449/1‑147 (MQ=255)
                                                                                                aTCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                 >  5:167966/1‑146 (MQ=255)
                                                                                                 tCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGACTCCCAGCGATGATAAACTTAACTAACACTATATTTAAGCATATTTCAGCACAAGCAAGGAGAACTCTTGACCCGAAACCAAAAAAACAATTTAATGACTTGAACgag                                              >  1:73463/1‑148 (MQ=255)
                                                                                                         gCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCATATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAGAAACTATTTAATGACTTGAAGGAGAAAAtcct                                      >  8:272409/1‑149 (MQ=255)
                                                                                                                     cagcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgat                          >  1:234265/1‑149 (MQ=255)
                                                                                                                      agcaACTTCCAGTGACTAGCCTCCCATCGAGGCTAAACTTAACTTACACTTTCTTTAAGCCTATTTCATCCCAATGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgat                          <  6:393630/148‑1 (MQ=255)
                                                                                                                           cTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTgg                    >  8:77766/1‑149 (MQ=255)
                                                                                                                            ttCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                   >  6:318473/1‑149 (MQ=255)
                                                                                                                            ttCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                   <  7:75376/149‑1 (MQ=255)
                                                                                                                            ttCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                   >  3:418187/1‑149 (MQ=255)
                                                                                                                               ttaagtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAATCCCAAAAAACTAGTTAATGCCTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgcgc                 >  7:73467/8‑148 (MQ=255)
                                                                                                                               cAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgcgc                 >  7:58051/1‑148 (MQ=255)
                                                                                                                               cAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgcgc                 <  4:57449/148‑1 (MQ=255)
                                                                                                                                   gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACACAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGTCGACAGc             >  4:257669/4‑148 (MQ=255)
                                                                                                                                    aCTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                     >  7:330796/1‑107 (MQ=255)
                                                                                                                                    aCTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                     <  8:330796/107‑1 (MQ=255)
                                                                                                                                    aCTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCa            >  5:248470/1‑148 (MQ=255)
                                                                                                                                        gCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaa         >  2:88406/1‑147 (MQ=255)
                                                                                                                                          cTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGa     <  7:241206/149‑1 (MQ=255)
                                                                                                                                           tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc    >  4:281271/1‑149 (MQ=255)
                                                                                                                                           tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc    >  5:430126/1‑149 (MQ=255)
                                                                                                                                           tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc    >  2:118622/1‑149 (MQ=255)
                                                                                                                                             ccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc    <  3:31621/147‑1 (MQ=255)
                                                                                                                                             ccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGAAAAAAAGACGt  >  4:40598/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167762‑6168050

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