Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,074,828 C→T 9.2% A28A (GCC→GCT 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,8280CT 9.2% 210.9 / 9.0 77A28A (GCC→GCTPP_4468Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base T (5/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

ATACCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCACCCCCTACCCAT  >  NC_002947/5074680‑5074976
                                                                                                                                                    |                                                                                                                                                    
ataCCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGcc                                                                                                                                                      <  1:391397/149‑1 (MQ=255)
 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)
                                                                                                                                                    |                                                                                                                                                    
ATACCTCATTTGAGTGTTTATCCGAAAGGCTTCTGTGTGGCTATATCCTGACTTTTGGCCAAACCATGCAGCTGAGAATTGCCCTAGCGGGCACACTTCGAGCCCTACGGCTAAACCGTCGCTTTCGGTACGAGCAGCTTTCTGACGCCAGCGTCAAATCTAAGATCAGCGCCCTTGAACGCGGGGAGGCCGGCATTACGCTGGAAAAATTCGAGAGCCTGGCTGAGGCACTTGAAATAGACCCGTTGGCCTTGTTTGCACTGTGCATGTCCGCACGGCGTGGCACCCCCTACCCAT  >  NC_002947/5074680‑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: