Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,432,709 T→A 7.8% L15M (TTG→ATG) ‡ PP_3049 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,432,7090TA 7.8% 181.2 / 4.8 64L15M (TTG→ATG) ‡PP_3049hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GTGGTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAG  >  NC_002947/3432566‑3432849
                                                                                                                                               |                                                                                                                                             
gtGGTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGAtc                                                                                                                                          <  5:273903/149‑1 (MQ=255)
   gTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCg                                                                                                                                       <  5:156321/149‑1 (MQ=255)
      gATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTg                                                                                                                                    >  6:133163/1‑149 (MQ=255)
            ccGCCTGATGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTgg                                                                                                                                <  6:158955/147‑1 (MQ=255)
              gCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGa                                                                                                                            >  1:142839/1‑149 (MQ=255)
                 agagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGAGATCAACCCGCTGc                                                                                                                                   <  5:85477/139‑1 (MQ=255)
                 agagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGAGATCAACCCGCTGc                                                                                                                                   >  6:85477/1‑139 (MQ=255)
                  gagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCg                                                                                                                        >  6:323289/1‑149 (MQ=255)
                     gACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCt                                                                                                                                     <  6:1149/133‑1 (MQ=255)
                     gACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCt                                                                                                                                     >  5:1149/1‑133 (MQ=255)
                           ggcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAg                                                                                                                                 >  3:168884/1‑131 (MQ=255)
                           ggcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAg                                                                                                                                 <  4:168884/131‑1 (MQ=255)
                            gcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGccccc                                                                                                               >  6:243875/1‑148 (MQ=255)
                                  ccTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGtt                                                                                                       <  4:49009/149‑1 (MQ=255)
                                    tACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGtttt                                                                                                     >  7:272612/1‑149 (MQ=255)
                                       acCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAg                                                                                                  <  2:142839/149‑1 (MQ=255)
                                         cATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGa                                                                                                <  5:133163/149‑1 (MQ=255)
                                          aTCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGAc                                                                                               >  2:93338/1‑149 (MQ=255)
                                                ggCAGTGAGGCTGGTGGGGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCa                                                                                                                      <  4:3692/121‑1 (MQ=255)
                                                ggCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCa                                                                                                                      >  3:3692/1‑121 (MQ=255)
                                                   aGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTggcc                                                                                                                              >  5:175410/1‑108 (MQ=38)
                                                   aGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTggcc                                                                                                                              <  6:175410/110‑3 (MQ=38)
                                                             gAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTccc                                                                              <  5:243875/147‑1 (MQ=255)
                                                                       ccTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCa                                                                                           <  2:347014/124‑1 (MQ=255)
                                                                       ccTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCa                                                                                           >  1:347014/1‑124 (MQ=255)
                                                                       ccTTGTAACGCTGATCTTGAGGAAGGCGGGATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCa                                                                  <  8:199287/149‑1 (MQ=255)
                                                                         ttGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATg                                                                >  3:74897/1‑149 (MQ=255)
                                                                           gTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAg                                                                                                  >  3:94645/1‑113 (MQ=255)
                                                                           gTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAg                                                                                                  <  4:94645/113‑1 (MQ=255)
                                                                            tAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                              <  8:272612/148‑1 (MQ=255)
                                                                            tAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGACGGCCAGCAGCGCGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTACAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcg                                                             >  5:282988/1‑149 (MQ=255)
                                                                                gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                           >  8:241861/1‑147 (MQ=255)
                                                                                gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                           >  7:140130/1‑147 (MQ=255)
                                                                                gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                           >  6:33403/1‑147 (MQ=255)
                                                                                   gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTTGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGc                                                       >  4:143661/1‑148 (MQ=255)
                                                                                   gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                              >  1:154395/1‑141 (MQ=255)
                                                                                   gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                              <  2:154395/141‑1 (MQ=255)
                                                                                       ttGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACgg                                                  >  7:42948/1‑149 (MQ=255)
                                                                                         gAGGAAGGCGGTATGAGTGTTGGTAGCTCTCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGAc                                                    <  7:241861/145‑1 (MQ=255)
                                                                                         gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                              >  6:6080/1‑135 (MQ=255)
                                                                                         gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACggg                                                 <  5:33403/148‑1 (MQ=255)
                                                                                         gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGGTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                              <  5:6080/135‑1 (MQ=255)
                                                                                                       gAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGAGCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                             <  1:112940/138‑1 (MQ=255)
                                                                                                       gAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                             >  2:112940/1‑138 (MQ=255)
                                                                                                         gtgtTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctg                                >  3:269601/1‑149 (MQ=255)
                                                                                                              ggTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                           >  3:85123/1‑149 (MQ=255)
                                                                                                               gTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTgg                          >  3:191486/1‑149 (MQ=255)
                                                                                                               gTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTgg                          >  2:133572/1‑149 (MQ=255)
                                                                                                                  gCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCtt                       <  4:85123/149‑1 (MQ=255)
                                                                                                                   cTATCACAAGCAGTTGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTc                      <  3:143661/149‑1 (MQ=255)
                                                                                                                    tATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGATTTTCAGGACTTCACGCGGCTGGTCCCCAAAgc                                                                        <  1:46777/99‑1 (MQ=255)
                                                                                                                    tATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAgc                                                                        >  2:46777/1‑98 (MQ=255)
                                                                                                                      tCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                             >  6:160852/1‑123 (MQ=255)
                                                                                                                      tCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                             <  5:160852/123‑1 (MQ=255)
                                                                                                                        acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                           <  2:322374/139‑1 (MQ=255)
                                                                                                                        acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                           >  1:322374/1‑139 (MQ=255)
                                                                                                                        acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAAc                  <  4:269601/148‑1 (MQ=255)
                                                                                                                           aGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCca               >  2:8634/1‑148 (MQ=255)
                                                                                                                                  ggCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGcc        <  7:249394/148‑1 (MQ=255)
                                                                                                                                  ggCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGcc        >  8:249394/1‑148 (MQ=255)
                                                                                                                                      agcagcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGa   <  4:191486/149‑1 (MQ=255)
                                                                                                                                      agcagcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGa   <  1:93338/149‑1 (MQ=255)
                                                                                                                                          gcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCCAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCGTGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAg  >  3:253962/1‑146 (MQ=255)
                                                                                                                                          gcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCGTGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAg  <  4:253962/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GTGGTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGTTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAG  >  NC_002947/3432566‑3432849

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