Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,996,897 G→T 5.3% P89H (CCT→CAT)  PP_3525 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,996,8970GT 5.3% 245.7 / 2.9 75P89H (CCT→CAT) PP_3525hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base T (2/2);  total (40/35)
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.85e-01

ACCAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAG  >  NC_002947/3996750‑3997043
                                                                                                                                                    |                                                                                                                                                   
aCCAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAgg                                                                                                                                                    >  6:250883/1‑149 (MQ=255)
 ccAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAgg                                                                                                                                                    >  8:274497/1‑148 (MQ=255)
   aGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCAACACGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGAATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGCCa                                                                                                                                                  >  1:116662/1‑149 (MQ=255)
   aGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCAt                                                                                                                                                 <  8:263617/149‑1 (MQ=255)
     cACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgc                                                                                                                                               >  8:41542/1‑149 (MQ=255)
     cACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgc                                                                                                                                               <  4:334822/149‑1 (MQ=255)
         cccTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACTTGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTg                                                                                                                                           <  1:315692/149‑1 (MQ=255)
               tgcAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTc                                                                                                                                     >  7:110742/1‑149 (MQ=255)
                      gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                               <  8:11990/148‑1 (MQ=255)
                      gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                               <  5:248505/148‑1 (MQ=255)
                       cAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGa                                                                                                                             <  6:111415/149‑1 (MQ=255)
                        aaCCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCGGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGaa                                                                                                                            <  4:79821/149‑1 (MQ=255)
                        aaCCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGaa                                                                                                                            <  8:66631/149‑1 (MQ=255)
                        aaCCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGATTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGaa                                                                                                                            >  7:309230/1‑149 (MQ=255)
                            aCTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATa                                                                                                                          >  2:68958/1‑147 (MQ=255)
                            aCTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCACCCATCACCGAAGAGCAACTCCCGGTAGCTGAAATAGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATa                                                                                                                          >  5:271501/1‑147 (MQ=255)
                              tGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgc                                                                                                                      >  4:31853/1‑149 (MQ=255)
                                cACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCATGAATACAgcg                                                                                                                     <  3:362696/148‑1 (MQ=255)
                                 aCCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                   >  7:251687/1‑149 (MQ=255)
                                  ccTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggcg                                                                                                                  >  5:354094/1‑149 (MQ=255)
                                  ccTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggcg                                                                                                                  <  5:129834/149‑1 (MQ=255)
                                   cTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGgcgc                                                                                                                 >  3:373600/1‑149 (MQ=255)
                                   cTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGgcgc                                                                                                                 >  4:391813/1‑149 (MQ=255)
                                      gCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGcgcg                                                                                                                <  8:204573/147‑1 (MQ=255)
                                          tcctGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTc                                                                                                          >  3:12737/1‑149 (MQ=255)
                                             tGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAg                                                                                                       >  8:76197/1‑149 (MQ=255)
                                                ggtgggTAGCGTCTGGCG‑ATTCGCCCATCACGTTCGAGCAGCACCCGGGAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTa                                                                                                    <  6:271501/146‑1 (MQ=255)
                                                ccGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                      >  2:177623/1‑147 (MQ=255)
                                                 cGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGgcgc                                                                                                                 <  7:345544/135‑1 (MQ=255)
                                                 cGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGgcgc                                                                                                                 >  8:345544/1‑135 (MQ=255)
                                                      tACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCtt                                                                                              <  2:116662/149‑1 (MQ=255)
                                                         cGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                           <  5:250883/149‑1 (MQ=255)
                                                           gCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGgct                                                                                         >  6:86136/1‑149 (MQ=255)
                                                            cTGACC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGgctg                                                                                        >  1:114225/1‑149 (MQ=255)
                                                               gCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑c                                                                                    <  7:274497/149‑1 (MQ=255)
                                                                cc‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTg                                                                                             >  4:338612/1‑140 (MQ=255)
                                                                cc‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTg                                                                                             <  3:338612/140‑1 (MQ=255)
                                                                cc‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑Cg                                                                                   <  8:110742/149‑1 (MQ=255)
                                                                   aCTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGc                                                                                 >  2:312906/1‑149 (MQ=255)
                                                                   aCTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGc                                                                                 >  2:91834/1‑149 (MQ=255)
                                                                       gCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGt                                                                             <  2:114225/149‑1 (MQ=255)
                                                                                 gTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGa                                                                   <  4:373600/149‑1 (MQ=255)
                                                                                 gTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGa                                                                   <  4:12737/149‑1 (MQ=255)
                                                                                    gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCa                                                                <  5:314691/149‑1 (MQ=255)
                                                                                      gcagcaCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGc                                                              >  3:288831/1‑149 (MQ=255)
                                                                                             ccGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCAc                                                       >  6:394952/1‑149 (MQ=255)
                                                                                                   gCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTCCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTcc                                                  <  1:68958/148‑1 (MQ=255)
                                                                                                      gTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACGAGCTTGGCACCGTCCatat                                              >  8:200708/1‑149 (MQ=255)
                                                                                                         aTCGCGCATGGCCGGGCCCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgc                                           >  5:292439/1‑149 (MQ=255)
                                                                                                         aTCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTa                                                                                                    <  4:196581/93‑1 (MQ=255)
                                                                                                         aTCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTa                                                                                                    >  3:196581/1‑93 (MQ=255)
                                                                                                            gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgcgcg                                        >  6:180601/1‑149 (MQ=255)
                                                                                                            gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgcgcg                                        >  1:106344/1‑149 (MQ=255)
                                                                                                                 tGGCCGGGACCGCGAACAGTTTGCTTATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                    <  1:91834/149‑1 (MQ=255)
                                                                                                                  ggCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCgcgcgc                                       >  8:183413/1‑144 (MQ=255)
                                                                                                                  ggCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCgcgcgc                                       <  7:183413/144‑1 (MQ=255)
                                                                                                                   gCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTgccgc                                                                        <  6:226108/110‑1 (MQ=255)
                                                                                                                       ggACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTaaa                             <  7:41542/149‑1 (MQ=255)
                                                                                                                        gACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTaaaa                            >  1:171259/1‑149 (MQ=255)
                                                                                                                            gcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTg                        >  8:95745/1‑149 (MQ=255)
                                                                                                                            gcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCCCCACCAGTAAAATCTg                        >  8:254264/1‑149 (MQ=255)
                                                                                                                             cgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGTGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGc                       <  5:152067/149‑1 (MQ=255)
                                                                                                                             cgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGc                       >  3:29002/1‑149 (MQ=255)
                                                                                                                             ccaaCCGTTTGCTGATCAGCGCCGGCAGGCGCTGAATGTCGGCCACCAATACAGCGGCGCGTGCTTCCCGGTAGCCCTTGGGCTGCT‑CGGCCCGTGCCGCCTTCCCCAGCTTGCCCCCGCCCATATCCCCCCCCACCACTAAAATcccc                       <  4:29002/147‑4 (MQ=255)
                                                                                                                               aaCAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCgt                     >  8:359506/1‑149 (MQ=255)
                                                                                                                                    tttGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCgg                >  4:256683/1‑149 (MQ=255)
                                                                                                                                     ccgacattCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCgg                >  6:305399/8‑148 (MQ=255)
                                                                                                                                          gATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAcc            >  2:308701/1‑147 (MQ=255)
                                                                                                                                          gATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAcc            >  6:188781/1‑147 (MQ=255)
                                                                                                                                            tCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAgg        >  5:299827/1‑149 (MQ=255)
                                                                                                                                            tCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAgg        <  2:106344/149‑1 (MQ=255)
                                                                                                                                             cAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAgg        <  6:299827/148‑1 (MQ=255)
                                                                                                                                             cAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGc       >  5:308434/1‑149 (MQ=255)
                                                                                                                                              aGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCt      <  8:251687/149‑1 (MQ=255)
                                                                                                                                                cgcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTg     >  4:115999/1‑148 (MQ=255)
                                                                                                                                                 gcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTa   <  3:31853/149‑1 (MQ=255)
                                                                                                                                                  cAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAg  <  1:177623/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCC‑ACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCT‑CGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAG  >  NC_002947/3996750‑3997043

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