Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,787,340 Δ1 bp 14.1% coding (1295/1392 nt) PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,787,3400T.14.1% 186.4 / 26.8 64coding (1295/1392 nt)PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base . (4/5);  total (32/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 = 6.61e-01

GCATGAATGATGCATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAA  >  NC_002947/3787194‑3787485
                                                                                                                                                    |                                                                                                                                                 
gcATGAATGATGCATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTc                                                                                                                                                  >  8:40699/1‑148 (MQ=255)
            cATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGGGGGCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTg                                                                                                                                     <  2:151218/149‑1 (MQ=255)
             aTGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGAGGGAGATCTTCGCCGCCGGCCACCATGATTGCCGAGCCCACGTCGTAGTACTTAAGGTCAAGCTGCAGGGTCAGCtt                                                                                                                                      >  1:293853/1‑149 (MQ=255)
                tgtgGCCCTCTGTTGTTGTTATCGGTCGGTTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGa                                                                                                                                  <  6:304236/148‑1 (MQ=255)
                  tgGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTgg                                                                                                                              <  5:137905/149‑1 (MQ=255)
                      ccTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCaa                                                                                                                           <  7:40699/149‑1 (MQ=255)
                      ccTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCGTCGCCGCCAGCCACCGTGATGGCCGTGCCCACGGAGTAGTACGTACGGTCACTGTCCAGGGGCAGGTTGGGAGTGGc                                                                                                                             >  4:291631/1‑147 (MQ=255)
                                 gttATCGGTCGGGTGGTCAAAGGCGGAACTTGGCGCGTACTGAGACAAACGACGAATCTTCGACGCCAGCCTCCTTGATGGCCGTGACCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCATGTTGGGAGTGGCAAGCCAGcacgc                                                                                                                 <  3:230095/148‑2 (MQ=255)
                                  ttATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGc                                                                                                               >  4:121327/1‑149 (MQ=255)
                                    aTCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGt                                                                                                                                <  4:10288/130‑1 (MQ=255)
                                    aTCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGt                                                                                                                                >  3:10288/1‑130 (MQ=255)
                                    aTCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTg                                                                                                             <  5:233665/148‑1 (MQ=255)
                                        gtcggGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgcgc                                                                                                         >  2:101690/1‑149 (MQ=255)
                                        gtcggGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgcgc                                                                                                         >  1:1957/1‑149 (MQ=255)
                                             gtggtCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTAATTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCgg                                                                                                       >  8:266515/1‑146 (MQ=255)
                                                   aaaGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgc                                                                                                           <  6:68525/136‑1 (MQ=255)
                                                   aaaGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgc                                                                                                           >  5:68525/1‑136 (MQ=255)
                                                   aaaGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTa                                                                                              <  3:121327/149‑1 (MQ=255)
                                                    aaGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgcgc                                                                                                         >  6:325117/1‑137 (MQ=255)
                                                    aaGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgcgc                                                                                                         <  5:325117/137‑1 (MQ=255)
                                                       gCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGa                                                                                          >  4:186136/1‑149 (MQ=255)
                                                           aaCGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATg                                                                                      >  4:48291/1‑149 (MQ=255)
                                                             cttGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATgg                                                                                     <  3:246104/146‑1 (MQ=255)
                                                                ggCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATgg                                                                                     >  4:246104/1‑145 (MQ=255)
                                                                 gcgcGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGcgcg                                                                                                                  >  2:249172/1‑115 (MQ=255)
                                                                 gcgcGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGcgcg                                                                                                                  <  1:249172/115‑1 (MQ=255)
                                                                          gagaCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTAGCGACATGGCGACCCGTGCcgcg                                                                       >  5:324448/1‑149 (MQ=255)
                                                                           agaCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCcgc                                                                        <  3:186136/147‑1 (MQ=255)
                                                                            gaCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGGGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCggg                                                                     <  4:308557/149‑1 (MQ=255)
                                                                             aCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgt                                                                    <  1:101690/149‑1 (MQ=255)
                                                                               aaaCGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTg                                                                                                                                     <  4:101078/82‑1 (MQ=255)
                                                                               aaaCGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTg                                                                                                                                     >  3:101078/1‑82 (MQ=255)
                                                                               aaaCGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCACGCCAGCGCGTGCTg                                                                                                             <  7:310424/106‑1 (MQ=255)
                                                                               aaaCGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGtt                                                                 >  2:67/1‑149 (MQ=255)
                                                                                 acgacgAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGTGTCAGGTTGGGAGTGGCAACCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGCGGCACGGgtcgcc                                                                <  6:340007/149‑5 (MQ=255)
                                                                                   gacgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGc                                                              <  7:65746/149‑1 (MQ=255)
                                                                                    acgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCCGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCa                                                            <  3:302742/149‑1 (MQ=255)
                                                                                    acgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGcc                                                             <  3:48291/149‑1 (MQ=255)
                                                                                    acgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGc                                                              >  4:190123/1‑148 (MQ=255)
                                                                                       gaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAg                                                           >  5:75491/1‑149 (MQ=255)
                                                                                       gaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAg                                                           >  4:40565/1‑148 (MQ=255)
                                                                                                   gccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACAt                                                                                       <  4:4018/110‑1 (MQ=255)
                                                                                                   gccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACAt                                                                                       >  3:4018/1‑110 (MQ=255)
                                                                                                     cAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAg                                                           <  8:257372/135‑1 (MQ=255)
                                                                                                     cAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAg                                                           >  7:257372/1‑135 (MQ=255)
                                                                                                         ctcctTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCa                                          >  3:258641/1‑149 (MQ=255)
                                                                                                            ctTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                       <  3:40565/148‑1 (MQ=255)
                                                                                                               gATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTaca                                    >  1:77807/1‑149 (MQ=255)
                                                                                                                    ccGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGa                                                                                                   >  7:63069/1‑81 (MQ=255)
                                                                                                                    ccGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGa                                                                                                   <  8:63069/81‑1 (MQ=255)
                                                                                                                    ccGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACTCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                       >  1:3142/1‑141 (MQ=255)
                                                                                                                    ccGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACTCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                       <  2:3142/141‑1 (MQ=255)
                                                                                                                         cccACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTc                          <  2:1957/149‑1 (MQ=255)
                                                                                                                         cccACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTc                          <  7:266515/149‑1 (MQ=255)
                                                                                                                              gtcgtAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCAc                     <  2:77807/149‑1 (MQ=255)
                                                                                                                                 gtagtaCTTCAGGTCAAGGCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTg                 >  6:338030/1‑149 (MQ=255)
                                                                                                                                  tagtaCTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGACGTCTTCCTCACGCTg                 >  5:366230/1‑149 (MQ=255)
                                                                                                                                        ttCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGc                                           <  2:222035/117‑1 (MQ=255)
                                                                                                                                        ttCAGGTCAAGGTCCAGGGTCAGGTTGGGAGAGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGc                                           >  1:222035/1‑117 (MQ=255)
                                                                                                                                             gTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACtt      >  2:92249/1‑149 (MQ=255)
                                                                                                                                              tCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTc     >  3:139222/1‑149 (MQ=255)
                                                                                                                                               cAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGtt               >  6:62361/1‑138 (MQ=255)
                                                                                                                                               cAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGtt               <  5:62361/138‑1 (MQ=255)
                                                                                                                                                 aGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGaa  <  1:92249/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCATGAATGATGCATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAA  >  NC_002947/3787194‑3787485

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: