Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,010,532 G→T 5.3% V14V (GTG→GTT)  PP_4415 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,010,5320GT 5.3% 241.7 / 3.5 75V14V (GTG→GTT) PP_4415hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base T (2/2);  total (43/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 = 8.34e-01

GGTATACTCGCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGATCGACTA  >  NC_002947/5010384‑5010672
                                                                                                                                                     |                                                                                                                                            
ggTATACTCGCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTg                                                                                                                                              <  7:2114/149‑1 (MQ=255)
  tataCTCGCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCGATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGgc                                                                                                                                            >  8:146782/1‑149 (MQ=255)
       tCGCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTc                                                                                                                                        >  1:207166/1‑148 (MQ=255)
       tCGCAACATACTGATATCATGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTc                                                                                                                                        >  8:62173/1‑148 (MQ=255)
         gCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAAc                                                                                                                                     >  7:60858/1‑149 (MQ=255)
             cATACTGATATCTTGTCGGATTCCTCCCCGATCGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCg                                                                                                                                 <  2:302608/149‑1 (MQ=255)
                    atatCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACaaa                                                                                                                          >  1:239368/1‑149 (MQ=255)
                      atCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGAACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAAc                                                                                                                         >  2:30410/1‑148 (MQ=255)
                         ttGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTc                                                                                                                               >  4:106571/1‑139 (MQ=255)
                         ttGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTc                                                                                                                               <  3:106571/139‑1 (MQ=255)
                         ttGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACg                                                                                                                     >  8:133050/1‑149 (MQ=255)
                          tGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACgt                                                                                                                    >  6:363065/1‑149 (MQ=255)
                              ggATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACgtggt                                                                                                                 <  8:60858/148‑1 (MQ=255)
                               gATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTTGAGCACCACGGCGCca                                                                                                                              >  3:69290/1‑134 (MQ=255)
                               gATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTTGAGCACCACGGCGCca                                                                                                                              <  4:69290/134‑1 (MQ=255)
                                    ctccCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAAt                                                                                                          <  7:62173/149‑1 (MQ=255)
                                       cccGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGc                                                                                                       <  5:363065/149‑1 (MQ=255)
                                               gCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGgc                                                                                                                                            >  2:285540/1‑104 (MQ=255)
                                               gCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGgc                                                                                                                                            <  1:285540/104‑1 (MQ=255)
                                               gCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATg                                                                                               >  1:243850/1‑149 (MQ=255)
                                                 ttGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGAGGGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCGCAAACTACTTGGTGCTGAATTGCCTGCAATg                                                                                               >  1:270061/1‑148 (MQ=255)
                                                 ttGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTg                                                                                                                                              >  7:148230/1‑100 (MQ=255)
                                                 ttGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTg                                                                                                                                              <  8:148230/100‑1 (MQ=255)
                                                      cATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCg                                                                                        >  2:105490/1‑149 (MQ=255)
                                                      cATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCg                                                                                        >  1:335135/1‑149 (MQ=255)
                                                      cATTCTACATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTACCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCg                                                                                        >  2:85382/1‑149 (MQ=255)
                                                       aTTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTACAATGGAAGTCGc                                                                                       >  1:123229/1‑149 (MQ=255)
                                                        ttCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTTCCTGCAAt                                                                                                >  4:144332/1‑139 (MQ=255)
                                                        ttCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAAt                                                                                                <  3:144332/139‑1 (MQ=255)
                                                             cATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCcac                                                                                    <  2:239368/146‑1 (MQ=255)
                                                              aTCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGc                                                                                       >  1:126111/1‑142 (MQ=255)
                                                              aTCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGc                                                                                       <  2:126111/142‑1 (MQ=255)
                                                              aTCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCt                                                                                >  7:284878/1‑149 (MQ=255)
                                                                  gtcgtgTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTc                                                                               >  2:180370/1‑146 (MQ=255)
                                                                      tcTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGc                                                                        <  8:218124/149‑1 (MQ=255)
                                                                          gATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTg                                                                      <  7:146782/147‑1 (MQ=255)
                                                                          gATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGc                                                                    >  4:252073/1‑149 (MQ=255)
                                                                              ttCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTcg                                                                >  5:55751/1‑149 (MQ=255)
                                                                                  ggTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTcggcgg                                                            <  2:243850/149‑1 (MQ=255)
                                                                                    tACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCggcggc                                                           >  3:352515/1‑148 (MQ=255)
                                                                                     aCCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGcc                                                          >  5:348943/1‑148 (MQ=255)
                                                                                         aTTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCa                                                     <  2:207166/149‑1 (MQ=255)
                                                                                          tttGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCag                                                    <  1:180370/149‑1 (MQ=255)
                                                                                          tttGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCag                                                    <  1:105490/149‑1 (MQ=255)
                                                                                             ggggA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACGTCGAGATGCTGGCTTCGGCGGCCAAGCagag                                                  <  7:303781/148‑1 (MQ=255)
                                                                                                gA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTc                                              >  3:148386/1‑149 (MQ=255)
                                                                                                      aCCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGACGCCGTTGAGCGCCACGGCgcc                                                                                                                               >  2:254049/1‑61 (MQ=21)
                                                                                                      aCCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGACGCCGTTGAGCGCCACGGCgcc                                                                                                                               <  1:254049/63‑3 (MQ=21)
                                                                                                       ccTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGaa                                                                                            <  4:329198/97‑1 (MQ=255)
                                                                                                       ccTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGaa                                                                                            >  3:329198/1‑97 (MQ=255)
                                                                                                       ccTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGt                                         >  2:248566/1‑148 (MQ=255)
                                                                                                        cTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGCCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGa                                       >  4:18061/1‑149 (MQ=255)
                                                                                                             tGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTAc                                  <  4:352515/149‑1 (MQ=255)
                                                                                                               cccACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACg                                 <  1:248566/148‑1 (MQ=255)
                                                                                                                  aCTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGaa                             >  7:211640/1‑149 (MQ=255)
                                                                                                                  aCTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGaa                             >  5:317133/1‑149 (MQ=255)
                                                                                                                    ttCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAAc                            >  3:268172/1‑148 (MQ=255)
                                                                                                                      cAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATgg                                                                                              <  3:142884/80‑1 (MQ=255)
                                                                                                                      cAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATgg                                                                                              >  4:142884/1‑80 (MQ=255)
                                                                                                                      cAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACg                           <  4:268172/147‑1 (MQ=255)
                                                                                                                       aaTTTTGTACCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTTCGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCTCACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGc                        <  1:85382/149‑1 (MQ=255)
                                                                                                                             gTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTcg                                             <  2:162520/122‑1 (MQ=255)
                                                                                                                             gTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTcg                                             >  1:162520/1‑122 (MQ=255)
                                                                                                                               tCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCa                 <  3:18061/148‑1 (MQ=255)
                                                                                                                                ccAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAg                >  3:337606/1‑148 (MQ=255)
                                                                                                                                 cAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGc                        >  6:6071/1‑139 (MQ=255)
                                                                                                                                 cAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGc                        <  5:6071/139‑1 (MQ=255)
                                                                                                                                  aGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAg             <  2:335135/149‑1 (MQ=255)
                                                                                                                                  aGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAg             >  4:14195/1‑149 (MQ=255)
                                                                                                                                  aGCTTGAGCATGGCGCCGTGGCGCTCAACGCCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCATCCCGGTGATTTACGAGAACGAGCATGACCAGGAg             >  6:41247/1‑149 (MQ=255)
                                                                                                                                    cTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCt           <  7:133050/149‑1 (MQ=255)
                                                                                                                                     ttGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCgag                                                                            <  1:203029/83‑1 (MQ=255)
                                                                                                                                     ttGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCgag                                                                            >  2:203029/1‑83 (MQ=255)
                                                                                                                                     ttGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCt           >  3:200241/1‑148 (MQ=255)
                                                                                                                                        aGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGAt        <  4:148386/148‑1 (MQ=255)
                                                                                                                                          cATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGATCGa     <  1:266998/149‑1 (MQ=255)
                                                                                                                                            tGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGATCGAc    >  6:137433/1‑148 (MQ=255)
                                                                                                                                             ggCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGATCGACTa  >  7:339846/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GGTATACTCGCAACATACTGATATCTTGTCGGATTCCTCCCCGATAGGCTTGACCATTCTCCATCCGTCGTCTTGATTTTCAGGTACCAATTTGGGGA‑GGTACCTGCATGCCCACTTCAATTTTGTTCCAGCTTGAGCATGGCGCCGTGGCGCTCAACGGCGTCACAAACTACGTGGTGCTGAATTGCCTGCAATGGAAGTCGCCACACCTCGAGATGCTGGCTTCGGCGGCCAAGCAGAGCTTCGCGGTGATTTACGAGAACGAGCTTGACCAGGAGCTGATCGACTA  >  NC_002947/5010384‑5010672

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