Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,210,758 C→A 6.1% intergenic (+35/+34) PP_3692 → / ← PP_5621 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,210,7580CA 6.1% 208.5 / 3.1 65intergenic (+35/+34)PP_3692/PP_5621hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (2/2);  total (32/33)
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.49e-01

TTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTC  >  NC_002947/4210610‑4210905
                                                                                                                                                    |                                                                                                                                                     
tttGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAAc                                                                                                                                                       >  8:299957/1‑149 (MQ=255)
      cATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGt                                                                                                                                                 <  7:299957/149‑1 (MQ=255)
          cccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                              <  4:113087/148‑1 (MQ=255)
           ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                            <  4:180473/149‑1 (MQ=255)
            cAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGc                                                                                                                                           >  6:211978/1‑149 (MQ=255)
              gTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTGACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCt                                                                                                                                          >  1:308769/1‑148 (MQ=255)
               tCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTa                                                                                                                                        <  2:193296/149‑1 (MQ=255)
                 ggaggaGCTTCCCAGCAGTGCCTTGCAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGc                                                                                                                                      <  3:119643/149‑1 (MQ=255)
                         ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTgg                                                                                                                              >  5:142033/1‑149 (MQ=255)
                           cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTgagc                                                                                                                                          >  4:135935/1‑132 (MQ=255)
                           cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTgagc                                                                                                                                          <  3:135935/135‑4 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTACCGAATTCTCAGCACGCAAGGAGCAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                            <  5:303799/148‑1 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTga                                                                                                                                            >  3:282117/1‑131 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTga                                                                                                                                            <  4:282117/132‑2 (MQ=255)
                              agcagTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgcagca                                                                                                                         >  7:313638/1‑149 (MQ=255)
                                 agTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATAc                                                                                                                      <  8:313638/149‑1 (MQ=255)
                                         ggAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCaa                                                                                                               >  7:186677/1‑148 (MQ=255)
                                                aaGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGAGTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAgc                                                                                                     <  8:280075/149‑1 (MQ=255)
                                                  ggTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAgcgg                                                                                                   <  2:315467/149‑1 (MQ=255)
                                                         gAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTcc                                                                                            >  4:189887/1‑149 (MQ=255)
                                                                aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTa                                                                                       <  2:400991/147‑1 (MQ=255)
                                                                aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTa                                                                                       <  3:318582/147‑1 (MQ=255)
                                                                  gTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCAt                                                                                   <  3:189887/149‑1 (MQ=255)
                                                                         cactgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAAt                                                                             <  5:211978/148‑1 (MQ=255)
                                                                           ctgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATAcc                                                                          >  1:175250/1‑149 (MQ=255)
                                                                           ctgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATAcc                                                                          >  2:158343/1‑149 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCAc                                                                                                                  <  3:74406/108‑1 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCAc                                                                                                                  >  4:74406/1‑108 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGa                                                                        >  1:189576/1‑148 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGCCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGa                                                                        >  1:358143/1‑148 (MQ=255)
                                                                               aTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                       >  1:274048/1‑148 (MQ=255)
                                                                                tAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTAGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCa                                                                     <  6:150247/149‑1 (MQ=255)
                                                                                   cTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGc                                                                  >  1:123595/1‑149 (MQ=255)
                                                                                        ttAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCg                                                             >  6:233198/1‑149 (MQ=255)
                                                                                         tAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                            >  7:99735/1‑149 (MQ=255)
                                                                                          aaCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                            <  2:358143/148‑1 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                            >  7:56584/1‑69 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                            <  8:56584/69‑1 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGc                                                           >  4:75655/1‑148 (MQ=255)
                                                                                            cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAC‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGa                                                         >  1:143159/1‑149 (MQ=255)
                                                                                             gAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGAt                                                        >  7:42177/1‑149 (MQ=255)
                                                                                               aTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGa                                                      <  8:42177/149‑1 (MQ=255)
                                                                                                        aCGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGa                                             >  4:251981/1‑149 (MQ=255)
                                                                                                                aGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTc                                     >  2:403435/1‑149 (MQ=255)
                                                                                                                   aGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGcc                                  <  1:158343/149‑1 (MQ=255)
                                                                                                                   aGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGc                                   <  3:75655/148‑1 (MQ=255)
                                                                                                                     cATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                                >  1:148993/1‑149 (MQ=255)
                                                                                                                     cATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                                <  2:175250/149‑1 (MQ=255)
                                                                                                                         cGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                             <  4:85992/148‑1 (MQ=255)
                                                                                                                         cGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                             <  2:123595/148‑1 (MQ=255)
                                                                                                                         cGGGGAGCGGGCTCAGAACCGCGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                             <  4:80512/148‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                            <  4:269683/131‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                            >  3:269683/1‑131 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGaa                          <  8:186677/149‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGaa                          <  5:233198/149‑1 (MQ=255)
                                                                                                                              aGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAAAGTGTAGCGGCGATTCCGATTACCATCCCAATACCGATAAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAg                         >  1:325107/1‑149 (MQ=255)
                                                                                                                                cGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGAc                     >  2:258417/1‑149 (MQ=255)
                                                                                                                                  ggCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGAt                                                        <  4:324468/112‑1 (MQ=255)
                                                                                                                                   gCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCa                  <  2:308769/149‑1 (MQ=255)
                                                                                                                                     tCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGc                >  1:355690/1‑149 (MQ=255)
                                                                                                                                     tCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGc                >  8:211449/1‑149 (MQ=255)
                                                                                                                                        gAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGa             >  1:148504/1‑149 (MQ=255)
                                                                                                                                               gTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTg         >  3:396/1‑146 (MQ=255)
                                                                                                                                                  aaCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGt                        >  4:59413/1‑128 (MQ=255)
                                                                                                                                                  aaCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGt                        <  3:59413/128‑1 (MQ=255)
                                                                                                                                                   aCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTc  <  2:148504/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTC  >  NC_002947/4210610‑4210905

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: