Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 846,736 G→C 6.1% T155S (ACT→AGT)  PP_0729 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947846,7360GC 6.1% 249.0 / 3.7 78T155S (ACT→AGT) PP_0729hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base C (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAGCTAG  >  NC_002947/846589‑846881
                                                                                                                                                   |                                                                                                                                                 
gTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGt                                                                                                                                                  <  2:15834/149‑1 (MQ=255)
 tGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCAAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGt                                                                                                                                                  >  5:114618/1‑148 (MQ=255)
  gATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGa                                                                                                                                                >  2:880/1‑149 (MQ=255)
   aTCGCGTACTCTATGTCCCAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTATTTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGa                                                                                                                                                <  2:274570/148‑1 (MQ=255)
      gcgAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCAGTACTGacggt                                                                                                                                            <  2:288130/149‑4 (MQ=255)
      gcgAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCAGTACTGacggt                                                                                                                                            >  5:102267/1‑146 (MQ=255)
       cgAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCg                                                                                                                                           >  4:293688/1‑149 (MQ=255)
        gAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGt                                                                                                                                          >  5:194005/1‑149 (MQ=255)
            tctATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAg                                                                                                                                       >  3:10051/1‑148 (MQ=255)
              tATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTAc                                                                                                                                    >  4:57280/1‑149 (MQ=255)
              tATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTAc                                                                                                                                    >  6:257238/1‑149 (MQ=255)
              tATGTCCAAGGCATCGCGCCCGGTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTACTGACACCGTCAGTACt                                                                                                                                   >  7:9244/1‑149 (MQ=255)
                  tCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTgg                                                                                                                                 >  6:205892/1‑148 (MQ=255)
                  tCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGc                                                                                                                                >  5:245836/1‑149 (MQ=255)
                   ccAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGt                                                                                                                                          >  3:48109/1‑138 (MQ=255)
                   ccAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGt                                                                                                                                          <  4:48109/138‑1 (MQ=255)
                    cAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCa                                                                                                                               <  5:257238/148‑1 (MQ=255)
                    cAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCAc                                                                                                                              >  2:158189/1‑149 (MQ=255)
                    cAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCAc                                                                                                                              <  1:880/149‑1 (MQ=255)
                    cAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCAGTACTGACGGTGTCACTACCGGCaa                                                                                                                              <  6:102267/149‑2 (MQ=255)
                    cAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCAGTACTGACGGTGTCACTACCGGCaa                                                                                                                              <  5:278562/149‑2 (MQ=255)
                     aaGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCAcc                                                                                                                             >  7:279803/1‑149 (MQ=255)
                       ggCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTg                                                                                                                           <  2:49281/149‑1 (MQ=255)
                       ggCATCGCGCCCGGTTAGAACTTGTACCGCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTg                                                                                                                           >  1:292850/1‑149 (MQ=255)
                        gCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTg                                                                                                                           >  8:141522/1‑148 (MQ=255)
                        gCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTg                                                                                                                           <  7:141522/148‑1 (MQ=255)
                          aTCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTa                                                                                                                        <  4:10051/149‑1 (MQ=255)
                          aTCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTa                                                                                                                        >  6:245721/1‑149 (MQ=255)
                              cgcCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATa                                                                                                                    <  3:276134/149‑1 (MQ=255)
                              cgcCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATa                                                                                                                    <  3:293688/149‑1 (MQ=255)
                              cgcCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATa                                                                                                                    <  8:279803/149‑1 (MQ=255)
                               gcCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAg                                                                                                                   >  4:276874/1‑149 (MQ=255)
                                  cGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTc                                                                                                                 >  1:157134/1‑148 (MQ=255)
                                         aaCTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGt                                                                                                         <  6:114618/149‑1 (MQ=255)
                                          aCTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTa                                                                                                        <  7:29510/149‑1 (MQ=255)
                                           cTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTaa                                                                                                       <  6:194005/149‑1 (MQ=255)
                                           cTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTaa                                                                                                       >  3:256138/1‑149 (MQ=255)
                                              gTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAata                                                                                                     <  2:157134/148‑1 (MQ=255)
                                              gTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAata                                                                                                     <  1:287258/148‑1 (MQ=255)
                                               tACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGTGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATTGTCAACTACGTAatat                                                                                                    <  7:79299/148‑1 (MQ=255)
                                               tACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAatat                                                                                                    <  5:205892/148‑1 (MQ=255)
                                                 ccTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACg                                                                                                          <  7:280152/140‑1 (MQ=255)
                                                 ccTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACg                                                                                                          >  8:280152/1‑140 (MQ=255)
                                                  ctctTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGc                                                                                                 >  6:287112/1‑148 (MQ=255)
                                                      tCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAAcc                                                                                            >  8:208752/1‑149 (MQ=255)
                                                        gAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAAccc                                                                                           <  7:208752/148‑1 (MQ=255)
                                                                 aaaCAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGAt                                                                                 >  8:82336/1‑149 (MQ=255)
                                                                 aaaCAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGAt                                                                                 <  4:273889/149‑1 (MQ=255)
                                                                   acaacCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGAGACCGTCAGTACTGGCACGTGGTAGATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCa                                                                               >  8:123539/1‑149 (MQ=255)
                                                                    caacCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGTAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCAc                                                                              <  1:158189/149‑1 (MQ=255)
                                                                      acCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACgt                                                                            <  8:101050/149‑2 (MQ=255)
                                                                      acCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACgg                                                                            <  3:276874/149‑1 (MQ=255)
                                                                      acCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACgg                                                                            <  6:245836/149‑1 (MQ=255)
                                                                          tGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCa                                                                        >  4:78061/1‑149 (MQ=255)
                                                                                ttCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCa                                                                  >  3:195328/1‑149 (MQ=255)
                                                                                  ccGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACt                                                                >  4:146007/1‑149 (MQ=255)
                                                                                    gCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTcc                                                              >  6:255645/1‑149 (MQ=255)
                                                                                        gCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTg                                                          >  6:1346/1‑149 (MQ=255)
                                                                                                  cGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATActc                                                <  3:146007/149‑1 (MQ=255)
                                                                                                    cTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATActcct                                              <  7:82336/149‑1 (MQ=255)
                                                                                                      tCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCAGTCTCTGGATCACGGGTCATCGCCACTACGCTGGCGAATActcctc                                             >  8:50584/1‑148 (MQ=255)
                                                                                                         cTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATGTGGCAAACCCGTCTCTGGATCACggg                                                                           >  4:268873/1‑115 (MQ=255)
                                                                                                         cTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATGTGGCAAACCCGTCTCTGGATCACggg                                                                           <  3:268873/115‑1 (MQ=255)
                                                                                                            ccAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCAGGCCAATACTCCTCTCGCCAt                                      >  8:38837/1‑149 (MQ=255)
                                                                                                                gggAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATctgc                                  >  7:12654/1‑149 (MQ=255)
                                                                                                                                 acaGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACatcatc                 >  5:56119/1‑149 (MQ=255)
                                                                                                                                 acaGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACatcatc                 <  4:176688/149‑1 (MQ=255)
                                                                                                                                 acaGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACatcatc                 <  4:195328/149‑1 (MQ=255)
                                                                                                                                  caGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGACACTCCGCTGGCCAATACTCCTCTCGCCATctgcagaac                             >  4:269805/1‑131 (MQ=255)
                                                                                                                                  caGTATCGTTGCCGGTAGTGACACCGTCAGTACTGGCTCCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCAc                             <  3:269805/136‑1 (MQ=255)
                                                                                                                                   aGTCTCGTTGCCGGTAGTTACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATcgc               <  7:50584/149‑1 (MQ=255)
                                                                                                                                        cGTTGCCGGTAGTGACACCGTCGGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTcg          >  1:220355/1‑149 (MQ=255)
                                                                                                                                             ccGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAg      <  8:12654/148‑1 (MQ=255)
                                                                                                                                             ccGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAGc     >  3:222976/1‑149 (MQ=255)
                                                                                                                                                gTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCAc                                                                              >  3:37961/1‑73 (MQ=255)
                                                                                                                                                gTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCAc                                                                              <  4:37961/73‑1 (MQ=255)
                                                                                                                                                gTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAGCTAg  <  3:57280/149‑1 (MQ=255)
                                                                                                                                                gTAGTGACACCGTCAGTACTGGCACCTGGCACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAGCTAg  <  5:287112/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCGTTGCCGGTAGTGACACCGTCAGTACTGGCACCTGGTACATAGTCAACTACGTAATATGGCAAACCCGTCTCTGGATCACGGGTCATCGCCACTCCGCTGGCCAATACTCCTCTCGCCATCTGCTGCACCAGCACATCATCGCGCTCGCGAGCTAG  >  NC_002947/846589‑846881

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: