Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,498,826 C→T 8.4% R393K (AGA→AAA)  PP_3988 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,498,8260CT 8.4% 165.7 / 3.5 60R393K (AGA→AAA) PP_3988hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base T (3/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

TATCTCGCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGC  >  NC_002947/4498678‑4498972
                                                                                                                                                      |                                                                                                                                                  
tATCTCGCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGActctc                                                                                                                                                    <  7:42619/149‑1 (MQ=255)
    tcGCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTg                                                                                                                                                  >  7:221118/1‑147 (MQ=255)
      gCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCGAAAGAGTGGGACTCTCTGGACt                                                                                                                                              >  4:148718/1‑149 (MQ=255)
        tCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCAA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACtt                                                                                                                                             >  1:309507/1‑148 (MQ=255)
          aaaTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAg                                                                                                                                          >  2:102382/1‑149 (MQ=255)
           aaTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGc                                                                                                                                         <  3:62564/149‑1 (MQ=255)
            aTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCt                                                                                                                                        <  6:23048/149‑1 (MQ=255)
                 gTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTtata                                                                                                                                    <  1:202822/148‑1 (MQ=255)
                 gTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTtata                                                                                                                                    <  8:221118/148‑1 (MQ=255)
                      gCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGTGACTCTCTGGACTTTAGCTTATAGCGCt                                                                                                                               <  2:366050/148‑1 (MQ=255)
                          gACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCATCCATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCtgtg                                                                                                                            <  1:102382/149‑1 (MQ=255)
                             ccACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATc                                                                                                                       >  3:324918/1‑149 (MQ=255)
                               acaAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAgcg                                                                                                                                 <  8:321637/137‑1 (MQ=255)
                               acaAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAgcg                                                                                                                                 >  7:321637/1‑137 (MQ=255)
                                caAAAACCAACGGCTTACTCACGGCATGTTCTAGGGATTTGAACTTCGCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CTTACACATTAAATTCCAAAGAGTGGGACTCTGTGGACTTTAGCTTATAGCGCTGTGCGATCGTc                                                                                                                    <  5:237222/149‑1 (MQ=255)
                                        aaCGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGcc                                                                                                            >  7:343755/1‑149 (MQ=255)
                                                 cTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGaa                                                                                                    <  6:218492/148‑1 (MQ=255)
                                                   cACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGca                                                                                                 >  1:224416/1‑149 (MQ=255)
                                                               tAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTg                                                                                     >  5:136247/1‑149 (MQ=255)
                                                               tAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGAGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTg                                                                                     >  6:299031/1‑149 (MQ=255)
                                                                aGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTa                                                                                                                                           >  2:220527/1‑94 (MQ=255)
                                                                aGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTa                                                                                                                                           <  1:220527/94‑1 (MQ=255)
                                                                 gggATTTGAACTTCCCGTAGCTCAACGGTCGATAGAGACGAACACCGTCCA‑‑CATACAAATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCTTCTGATCGCCCTCGAGAAGCACATCAACCGCTGcc                                                                                   >  1:161167/1‑149 (MQ=255)
                                                                 gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                >  2:275147/1‑136 (MQ=255)
                                                                 gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                >  2:280092/1‑136 (MQ=255)
                                                                 gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                <  1:275147/136‑1 (MQ=255)
                                                                 gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                <  1:280092/136‑1 (MQ=255)
                                                                              cccGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGttt                                                                        <  4:324918/147‑1 (MQ=255)
                                                                                        aaCGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAAttt                                                             <  2:309507/148‑1 (MQ=255)
                                                                                              cGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAAt                                                      >  7:219475/1‑149 (MQ=255)
                                                                                                aTAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATc                                                     >  4:297191/1‑148 (MQ=255)
                                                                                                  agagACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGATGTGCGATCGTATGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGa                                                   >  1:328483/1‑148 (MQ=255)
                                                                                                     gaCGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGATTGTTATCAATATGAAAAATCGAAGcc                                               >  6:14878/1‑149 (MQ=255)
                                                                                                          acacCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGt                                           >  2:47773/1‑148 (MQ=255)
                                                                                                           cacCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCAAATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCa                                         >  3:93179/1‑149 (MQ=255)
                                                                                                              cGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCTTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTg                                      <  2:224416/149‑1 (MQ=255)
                                                                                                              cGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTg                                      >  1:324719/1‑149 (MQ=255)
                                                                                                               gTCCA‑‑CATACACATTAAATTCCAAAGAGTCGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTgg                                     <  3:40782/149‑1 (MQ=255)
                                                                                                                tCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTgg                                     >  2:90749/1‑148 (MQ=255)
                                                                                                                tCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATAgg                                     >  4:83319/1‑148 (MQ=255)
                                                                                                                 ccA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCa                                   >  7:112426/1‑149 (MQ=255)
                                                                                                                       aTACACATTAAATTCCAAAGAGTGGGACTATCTGGACTTTAGCTTATAGCGCTGTGCGATAGTCTGATCGCCATCGAGAAGCACATCAACCGCTGACCCTCGCCCTTTGTTATCAATTTGAAAAATCGATGCCCAGTCATTGGCATCa                                >  8:127220/1‑148 (MQ=255)
                                                                                                                        tACACATTAAATTCCAAAGAGTGGGACTCTCAGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTCTTATCAATTTGAAAAATCGATGCCCCGTCATTGCcaacaa                               >  7:65766/1‑148 (MQ=255)
                                                                                                                           acaTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATc                           <  8:112426/149‑1 (MQ=255)
                                                                                                                           acaTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATc                           >  1:23045/1‑149 (MQ=255)
                                                                                                                             gctAAAGTCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  8:287651/3‑148 (MQ=38)
                                                                                                                             gctAAAGTCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  1:124702/3‑148 (MQ=38)
                                                                                                                             aTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  7:115408/1‑148 (MQ=255)
                                                                                                                             aTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  1:242082/1‑148 (MQ=255)
                                                                                                                                 aaTTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAg                      <  6:136247/148‑1 (MQ=255)
                                                                                                                                    tCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTg                  >  1:192606/1‑149 (MQ=255)
                                                                                                                                      cAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGct                >  3:248689/1‑149 (MQ=255)
                                                                                                                                       aaaGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGctt               >  2:118443/1‑149 (MQ=255)
                                                                                                                                          gagTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCATATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGcttctt            <  4:93179/149‑1 (MQ=255)
                                                                                                                                           agTTGGACTCTCTGGACTTTATCATATAGCGCTGTGCGATCGTCTGATCTCCCTCGAGAATCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCTTCATTGGCATCAATATCGAAAGGGTGCTTCTTa           <  7:127220/149‑1 (MQ=255)
                                                                                                                                            gTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTaa          <  4:248689/149‑1 (MQ=255)
                                                                                                                                            gTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTaa          <  8:343755/149‑1 (MQ=255)
                                                                                                                                             tcccaCTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAg         <  7:287651/145‑1 (MQ=255)
                                                                                                                                               ccaCTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGAt       <  2:192606/147‑1 (MQ=255)
                                                                                                                                                 aCTCTCTGGACTTTTGGTTATAGCGCTGTGCGATCGTCTGATCGCCGTCGAGAAGCACATGTACCGCTGCCGCTCGGCGTTTGTTATCATTTTGAATAATCGATGCCCCGTCATTGGGATCAATATCGAAAGGGTGCTTCTTAAGATAg     <  8:65766/149‑1 (MQ=255)
                                                                                                                                                 aCTCTCTGGACTTTAGCTTATAGCGCTGTGCTATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAg     <  8:219475/149‑1 (MQ=255)
                                                                                                                                                  ctctctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCACTCGACGTTTTTTATAAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAgt    <  2:161167/149‑2 (MQ=255)
                                                                                                                                                     tctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGc  >  2:276311/1‑148 (MQ=255)
                                                                                                                                                     tctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGc  >  7:157433/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TATCTCGCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGC  >  NC_002947/4498678‑4498972

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