Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,017,449 T→A 5.8% E138V (GAA→GTA)  PP_5461 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,017,4490TA 5.8% 210.2 / 2.4 69E138V (GAA→GTA) PP_5461hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/29);  new base A (2/2);  total (38/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 = 4.22e-01

GGTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGT  >  NC_002947/1017305‑1017596
                                                                                                                                                 |                                                                                                                                                   
ggTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCCCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtga                                                                                                                                                 >  4:5909/1‑149 (MQ=255)
 gTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgat                                                                                                                                                <  3:157421/149‑1 (MQ=255)
  tGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgatg                                                                                                                                               <  8:114043/149‑1 (MQ=255)
   gATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgatg                                                                                                                                               >  5:380150/1‑148 (MQ=255)
        ttCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCGTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAacac                                                                                                                                          >  6:311623/1‑148 (MQ=255)
             aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCaaa                                                                                                                                      <  1:430249/147‑1 (MQ=255)
             aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCaa                                                                                                                                       <  4:331559/146‑1 (MQ=255)
                tGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCt                                                                                                                                 >  2:433434/1‑149 (MQ=255)
                 ggTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTc                                                                                                                                >  1:217279/1‑149 (MQ=255)
                 ggTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTc                                                                                                                                >  2:244584/1‑149 (MQ=255)
                   ttCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGt                                                                                                                              >  4:152797/1‑149 (MQ=255)
                      tttATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGa                                                                                                                            >  4:349383/1‑148 (MQ=255)
                         aTGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACAt                                                                                                                         <  1:433434/148‑1 (MQ=255)
                             aTTGAGATTTTGGGGGCCGCTTTGAGTCCCTTTTGCTACAACTCTTACCAACAAAAGCCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTTAGATCAATAAATACAGTTTGGTGTACTGATGAACACCAAAGCCCTCATGACATTttct                                                                                                                    >  5:110561/1‑149 (MQ=255)
                                          gggCCGCTTTGCGTCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAAt                                                                                                       <  3:240389/149‑1 (MQ=255)
                                               gCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAg                                                                                                  <  2:217279/149‑1 (MQ=255)
                                                    gaggCCATTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCGAGTTTGGTGTACTGATGATCACCTAAGGCCTCGTGTGATTTTCTTCGGTATATAAATTGTAGAttt                                                                                              <  6:110561/146‑1 (MQ=255)
                                                        cccTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTggg                                                                                         <  3:5909/149‑1 (MQ=255)
                                                               gCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgtgt                                                                                  >  7:206919/1‑149 (MQ=255)
                                                                 tACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGTGGGAGTGTGc                                                                                >  1:262212/1‑149 (MQ=255)
                                                                 tACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGc                                                                                >  7:62437/1‑149 (MQ=255)
                                                                 tACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTACTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGt                                                                                                    <  2:267803/128‑1 (MQ=255)
                                                                 tACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTACTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGt                                                                                                    >  1:267803/1‑128 (MQ=255)
                                                                   cAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGcc                                                                               >  4:61930/1‑148 (MQ=255)
                                                                    aaCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCaa                                                                             >  4:179446/1‑149 (MQ=255)
                                                                     aCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGGGTGTGCCAAt                                                                            <  1:149000/149‑1 (MQ=255)
                                                                           cGCTACAACAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGc                                                                         <  4:239487/146‑1 (MQ=255)
                                                                           cGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGc                                                                         <  2:328674/146‑1 (MQ=255)
                                                                           cGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTa                                                                      >  8:219345/1‑149 (MQ=255)
                                                                             cTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTact                                                                    >  6:102561/1‑147 (MQ=255)
                                                                                cAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttt                                                                   <  6:380150/147‑1 (MQ=255)
                                                                                cAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttt                                                                   <  3:152797/147‑1 (MQ=255)
                                                                                 taaaGGCTTGTGAAAGGGAGGGCCCGGA‑TGTCGGGCTTGAGTTAAATAAATACAGTTTTGTGTTATGATGAACACCAAGGGCGGCGTGAATTTTTGTTGGGTAGATATATTGTATGTTTTTGTGGGAGTGTGGGAATGGCGTATTTTTc                                                                >  4:359932/2‑149 (MQ=255)
                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGAGTGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttttt                                                                 >  1:27095/1‑149 (MQ=255)
                                                                                 aaaaGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTc                                                                >  5:337488/1‑149 (MQ=255)
                                                                                 aaaaGGCGTGTGCAAGGCAGCGCCCGGA‑TGTCGGCCTTGAGATCAATAAATCCAGTTTGGTGTTCTGAAGAACACCAAAGCCCTCGAGCCATTTTCTTCGGTATATAAATTGTAGATATTTGGGGGAGAGTGCCAATGGCGTATTTTTc                                                                >  6:423139/1‑149 (MQ=255)
                                                                                   aaGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGc                                                              >  5:184266/1‑149 (MQ=255)
                                                                                     ggCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTg                                                            <  8:62437/149‑1 (MQ=255)
                                                                                         tgtgCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCggg                                                        <  7:387995/149‑1 (MQ=255)
                                                                                          gtgCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                       <  6:337488/149‑1 (MQ=255)
                                                                                            gCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCt                                                     >  6:236002/1‑149 (MQ=255)
                                                                                              aaGGCAGGGCCCGTA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGt                                                   <  1:367385/149‑1 (MQ=255)
                                                                                               aGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTa                                                  >  7:1112/1‑149 (MQ=255)
                                                                                                ggCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAg                                                 >  8:404290/1‑149 (MQ=255)
                                                                                                    gggCCCGGA‑TGTCGGGGTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGatat                                             <  1:328811/149‑1 (MQ=255)
                                                                                                    gggCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCTGTTTGGTTTTCTGATGAACACCAAAGGCCTCGTGTCATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTTCCAATGGCGTTTTTTTCGCTGCGGGCCTGTAGatat                                             <  8:111344/149‑1 (MQ=255)
                                                                                                        ccGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGt                                         >  8:145201/1‑149 (MQ=255)
                                                                                                          ggA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGc                                       >  1:200023/1‑149 (MQ=255)
                                                                                                              tGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCa                                    >  1:406863/1‑149 (MQ=255)
                                                                                                               gTCGGGCTTGAGATCAATAAATCAAGTTTGATGTTCTAATGAACACCAAAGACATAATAAAATTTTCTTAAGTATATAAATTGTAAATTTTTTGGGGAGTGTGAAAATGAAGTATTTTTAGCTAAGGGCATGAAAATATGAGTgataa                                    >  8:272806/1‑144 (MQ=255)
                                                                                                                tCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCAt                                   <  5:274605/148‑1 (MQ=255)
                                                                                                                 cGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAg                                 >  8:18850/1‑149 (MQ=255)
                                                                                                                        gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGTGCCTGTAGATATGAGTGCTCATAGTAATGGa                          >  7:101327/1‑149 (MQ=255)
                                                                                                                        gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATgg                           <  3:349383/148‑1 (MQ=255)
                                                                                                                        gagaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGa                          >  8:285021/1‑149 (MQ=255)
                                                                                                                          gaTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACt                        <  3:61930/149‑1 (MQ=255)
                                                                                                                              aaTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATTGACTATc                     <  5:311623/148‑1 (MQ=255)
                                                                                                                              aaTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATc                     >  7:95534/1‑148 (MQ=255)
                                                                                                                               aTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCaa                   >  6:123909/1‑149 (MQ=255)
                                                                                                                               aTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCa                    >  7:201851/1‑148 (MQ=255)
                                                                                                                                  aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAat                 <  7:219345/148‑1 (MQ=255)
                                                                                                                                      cAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatctatc            <  7:145201/149‑1 (MQ=255)
                                                                                                                                       aGTTTGGTGTTCTGATGAACAACAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGATTGCTCATAGTAATGGACTATCAAatctatct           <  2:262212/149‑1 (MQ=255)
                                                                                                                                        gTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCtt          <  5:236002/149‑1 (MQ=255)
                                                                                                                                         tttGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg         >  4:166970/1‑149 (MQ=255)
                                                                                                                                         tttGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg         <  8:201851/149‑1 (MQ=255)
                                                                                                                                         tttGGTGTTCATCAGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg         >  3:267402/1‑149 (MQ=255)
                                                                                                                                          ttGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg         <  7:18850/148‑1 (MQ=255)
                                                                                                                                             gtgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaa      >  4:178633/1‑148 (MQ=255)
                                                                                                                                                ttCTGATCAACACCAAAGGCCTCGAGACATTTTCTTCGGTATATAAATAGTAGATTTTTGGGGGAGAGTGACAATGGCGTATTTATCGCTGCGGGCCTGTAGAAATGAGGGCTCATAGTAAGGCAGTATCAAAACTATCTTGCaaaact  >  8:291582/1‑147 (MQ=255)
                                                                                                                                                 tCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGtatat                                                                                                           <  4:78124/43‑1 (MQ=255)
                                                                                                                                                 tCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGtatat                                                                                                           >  3:78124/1‑43 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGA‑TGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGT  >  NC_002947/1017305‑1017596

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