Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,074,829 A→G 9.1% S29G (AGC→GGC)  PP_4468 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,074,8290AG 9.1% 210.4 / 9.6 77S29G (AGC→GGC) PP_4468Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base G (5/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCACCCCCTACCCAT  >  NC_002947/5074681‑5074976
                                                                                                                                                    |                                                                                                                                                   
taCCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCa                                                                                                                                                     <  8:293469/149‑1 (MQ=255)
 aCCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAg                                                                                                                                                    >  8:186477/1‑149 (MQ=255)
  ccTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTATTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGc                                                                                                                                                   >  1:109018/1‑149 (MQ=255)
      aTTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTc                                                                                                                                                <  6:316155/148‑1 (MQ=255)
      aTTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTc                                                                                                                                                <  1:357871/148‑1 (MQ=255)
          gAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAAt                                                                                                                                            <  7:367108/148‑1 (MQ=255)
             tgtTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTaa                                                                                                                                        <  5:200257/149‑1 (MQ=255)
             tgtTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTaa                                                                                                                                        <  7:195388/149‑1 (MQ=255)
               tttATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAg                                                                                                                                       <  5:258623/148‑1 (MQ=255)
                  aTCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCTGGCGTCagaaag                                                                                                                                          <  4:255668/142‑6 (MQ=255)
                  aTCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCTGGCGTCagaaag                                                                                                                                          >  3:255668/1‑137 (MQ=255)
                    ccGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAAt                                                                                                                                            <  1:256937/138‑1 (MQ=255)
                    ccGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAAt                                                                                                                                            >  2:256937/1‑138 (MQ=255)
                      gAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCTGGCGTCaga                                                                                                                                             <  6:361788/135‑3 (MQ=255)
                      gAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCTGGCGTCaga                                                                                                                                             >  5:361788/1‑133 (MQ=255)
                          ggCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGGCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCtt                                                                                                                           <  8:13057/149‑1 (MQ=255)
                             ttCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAAt                                                                                                                                            <  6:238998/129‑1 (MQ=255)
                             ttCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAAt                                                                                                                                            >  5:238998/1‑129 (MQ=255)
                                  tgtgGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCgg                                                                                                                   <  5:394088/149‑1 (MQ=255)
                                         atatCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCt                                                                                                                                          <  3:95474/119‑1 (MQ=255)
                                         atatCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCt                                                                                                                                          >  4:95474/1‑119 (MQ=255)
                                          tatCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTCGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCg                                                                                                           <  4:54557/149‑1 (MQ=255)
                                          tatCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGcc                                                                                                            <  5:272440/148‑1 (MQ=255)
                                            tCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCgg                                                                                                          >  2:270555/1‑148 (MQ=255)
                                             ccTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGt                                                                                                                                                 <  6:229816/108‑1 (MQ=255)
                                             ccTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGt                                                                                                                                                 >  5:229816/1‑108 (MQ=255)
                                              cTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTGCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGCACGCGGGGAGGCCGGCa                                                                                                        <  7:186477/148‑1 (MQ=255)
                                               tGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCAt                                                                                                       >  2:242345/1‑148 (MQ=255)
                                                   ttttGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGc                                                                                                  >  5:37653/1‑149 (MQ=255)
                                                   ttttGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGc                                                                                                  >  7:86861/1‑149 (MQ=255)
                                                          cAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGaaa                                                                                            >  5:51140/1‑148 (MQ=255)
                                                           aaaCCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGaaaaa                                                                                          >  8:405353/1‑149 (MQ=255)
                                                            aaCCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAAt                                                                                         >  8:257647/1‑149 (MQ=255)
                                                                      cTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGGCGGCATTACGCTGGAAAAATTCGAGAGCCt                                                                               >  1:262991/1‑149 (MQ=255)
                                                                       tGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTg                                                                              >  4:394709/1‑149 (MQ=255)
                                                                       tGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTg                                                                              >  8:139507/1‑149 (MQ=255)
                                                                         agaATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGc                                                                            >  3:228836/1‑149 (MQ=255)
                                                                            aTTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCt                                                                           >  8:130034/1‑147 (MQ=255)
                                                                             ttGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAg                                                                        <  4:228836/149‑1 (MQ=255)
                                                                                cccTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCa                                                                     >  5:407902/1‑149 (MQ=255)
                                                                                 ccTAGCGGGCACAGTTCGAGCCCTAAGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTATGATCAGCGCCCGTGAACGCGGGGAGGCCGGCAGTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGcgg                                                                    >  8:409349/1‑147 (MQ=255)
                                                                                   tAGCGGGAACACTTAGAGCCCTACGGCTAACACGGCGCTTTCTGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACtt                                                                  <  8:390350/149‑1 (MQ=255)
                                                                                     gCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACtt                                                                  <  1:242345/147‑1 (MQ=255)
                                                                                       gggCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCg                                                                                                           <  4:291148/104‑1 (MQ=255)
                                                                                       gggCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCg                                                                                                           >  3:291148/1‑104 (MQ=255)
                                                                                       gggCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGaaa                                                              >  4:365104/1‑149 (MQ=255)
                                                                                       gggCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGaaa                                                              >  7:31869/1‑149 (MQ=255)
                                                                                          cacaCTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAg                                                           <  7:257647/149‑1 (MQ=255)
                                                                                            cacTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCTGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGAc                                                         <  2:109018/149‑1 (MQ=255)
                                                                                            cacTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGAc                                                         <  8:86861/149‑1 (MQ=255)
                                                                                             acTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGAcc                                                        <  2:210937/149‑1 (MQ=255)
                                                                                             acTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGAcc                                                        <  2:217562/149‑1 (MQ=255)
                                                                                             acTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGAcc                                                        >  3:392492/1‑149 (MQ=255)
                                                                                               ttCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCg                                                      >  6:229419/1‑149 (MQ=255)
                                                                                                        tACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTg                                             >  1:404582/1‑149 (MQ=255)
                                                                                                          cGGCTAAACCGTCGCTTTCGGTACGATCAGCTTGCTGACGCGAGCGTCAAATCAAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGt                                            <  3:103938/148‑1 (MQ=255)
                                                                                                          cGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGt                                            <  7:139507/148‑1 (MQ=255)
                                                                                                           ggCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGttt                                          <  7:130034/149‑1 (MQ=255)
                                                                                                           ggCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGttt                                          <  7:405353/149‑1 (MQ=255)
                                                                                                              tAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAg                                                                                                                                  >  1:157469/1‑58 (MQ=255)
                                                                                                              tAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAg                                                                                                                                  <  2:157469/58‑1 (MQ=255)
                                                                                                                aaCCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACt                                     <  2:27547/149‑1 (MQ=255)
                                                                                                                      cGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCAt                               <  1:270555/149‑1 (MQ=255)
                                                                                                                         tttCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTc                            >  6:334920/1‑149 (MQ=255)
                                                                                                                            cGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCg                          <  8:31869/148‑1 (MQ=255)
                                                                                                                              gTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCAc                       >  3:403692/1‑149 (MQ=255)
                                                                                                                              gTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCAc                       <  7:264507/149‑1 (MQ=255)
                                                                                                                               tACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACg                      >  6:126779/1‑149 (MQ=255)
                                                                                                                                   agcagcTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGt                  >  3:173056/1‑149 (MQ=255)
                                                                                                                                    gcagcTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTg                 >  3:429006/1‑149 (MQ=255)
                                                                                                                                    gcagcTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCCGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTg                 <  4:391101/149‑1 (MQ=255)
                                                                                                                                      ttagatttGACGCTGGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGc               >  8:106958/8‑149 (MQ=255)
                                                                                                                                       tagatttGACGCTGGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCa              >  6:349133/7‑149 (MQ=255)
                                                                                                                                         tttCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCAcc            >  6:293699/1‑149 (MQ=255)
                                                                                                                                         tttCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCAcc            >  7:148773/1‑149 (MQ=255)
                                                                                                                                         gatttGACGCTGGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCAcc            >  8:172550/5‑149 (MQ=255)
                                                                                                                                             tGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCAccccc         >  2:390671/1‑148 (MQ=255)
                                                                                                                                                   cAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCACCCCCTACCCAt  >  5:99656/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCACCCCCTACCCAT  >  NC_002947/5074681‑5074976

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: