Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,690,931 C→T 6.6% A22V (GCA→GTA)  PP_4151 → excisionase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,690,9310CT 6.6% 176.0 / 2.2 60A22V (GCA→GTA) PP_4151excisionase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base T (2/2);  total (34/26)
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.18e-01

GCGATTTGCGAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATTCTCA  >  NC_002947/4690784‑4691071
                                                                                                                                                   |                                                                                                                                            
gCGATTTGCGAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGc                                                                                                                                              <  4:80297/148‑1 (MQ=255)
         gAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGt                                                                                                                                      <  6:333635/147‑1 (MQ=255)
         gAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGt                                                                                                                                      >  5:333635/1‑147 (MQ=255)
         gAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTc                                                                                                                                     >  8:329223/1‑148 (MQ=255)
          aaTATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCaa                                                                                                                                   <  5:144491/149‑1 (MQ=255)
                          cTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCAAGTCAAGCCGCACGCTGGCAGc                                                                                                                   >  6:140533/1‑149 (MQ=255)
                                     gCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGa                                                                                                        >  6:255244/1‑149 (MQ=255)
                                     gCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGa                                                                                                        <  5:255244/149‑1 (MQ=255)
                                      cTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGa                                                                                                        >  4:239883/1‑148 (MQ=255)
                                             aTCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCa                                                                                                <  1:180083/149‑1 (MQ=255)
                                                     cacaGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         >  6:239370/1‑148 (MQ=255)
                                                     cacaGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         >  5:75902/1‑148 (MQ=255)
                                                         gCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACgg                                                                                    <  1:136938/149‑1 (MQ=255)
                                                             cGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCCACGGATc                                                                                 >  8:304941/1‑148 (MQ=255)
                                                             cGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATc                                                                                 >  1:44042/1‑148 (MQ=255)
                                                             cGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATc                                                                                 >  4:147357/1‑148 (MQ=255)
                                                               cAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGTAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         >  8:152335/1‑138 (MQ=255)
                                                               cAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGTAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         <  7:152335/138‑1 (MQ=255)
                                                               cAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGTAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         <  7:157060/138‑1 (MQ=255)
                                                               cAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGTAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATAccc                                                                                         >  8:157060/1‑138 (MQ=255)
                                                               cAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGaa                                                                              <  7:304941/149‑1 (MQ=255)
                                                                    gATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGc                                                                          >  5:277578/1‑148 (MQ=255)
                                                                         cGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAAc                                                                     >  1:23842/1‑148 (MQ=255)
                                                                          gCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGa                                                                   <  6:98708/149‑1 (MQ=255)
                                                                               ggAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGaa                                                               >  6:256022/1‑148 (MQ=255)
                                                                                   aTGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAAc                                                                                      <  7:240529/121‑1 (MQ=255)
                                                                                   aTGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCAGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCATAAACCAAAc                                                                                      >  8:240529/1‑121 (MQ=255)
                                                                                    tGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCa                                                          <  6:277578/148‑1 (MQ=255)
                                                                                        cATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAgcgc                                                      >  5:121486/1‑148 (MQ=255)
                                                                                         aTGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAgcgc                                                      >  7:131440/1‑147 (MQ=255)
                                                                                          tGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAg                                                   <  7:329223/149‑1 (MQ=255)
                                                                                             ccaccGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTAAACGAGGTAAATCAGCGCCAGgtt                                                >  8:148783/1‑149 (MQ=255)
                                                                                              caccGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGgtt                                                <  8:282753/148‑1 (MQ=255)
                                                                                              caccGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGgtt                                                >  7:282753/1‑148 (MQ=255)
                                                                                               accGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGa                                                                                                        <  6:92381/91‑1 (MQ=255)
                                                                                               accGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGgttgt                                              <  5:46489/149‑1 (MQ=255)
                                                                                                  gCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAATTCAGCGCAAGGttgttg                                            >  8:128598/1‑148 (MQ=255)
                                                                                                  gCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAgtattttg                                            >  8:123586/1‑141 (MQ=255)
                                                                                                   ccAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGaa                                          <  2:307034/149‑1 (MQ=255)
                                                                                                   ccAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGaa                                          <  2:44042/149‑1 (MQ=255)
                                                                                                    cAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGa                                                                                              >  2:300605/1‑96 (MQ=255)
                                                                                                       ccTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGc                                                                          >  3:166183/1‑113 (MQ=255)
                                                                                                       ccTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGc                                                                          <  4:166183/113‑1 (MQ=255)
                                                                                                       ccTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTAc                                      <  3:239883/149‑1 (MQ=255)
                                                                                                       ccTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCAAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTAc                                      <  5:140533/149‑1 (MQ=255)
                                                                                                           tctAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTAc                                  >  3:129051/1‑149 (MQ=255)
                                                                                                                            tGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTGCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGTAAAATCAGCTCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTa                  >  7:337486/1‑148 (MQ=255)
                                                                                                                            tGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTAt                 >  5:329844/1‑149 (MQ=255)
                                                                                                                            tGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTAt                 >  5:182659/1‑149 (MQ=255)
                                                                                                                             gACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTAtt                >  6:305779/1‑149 (MQ=255)
                                                                                                                              aCGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTAtt                <  6:121486/148‑1 (MQ=255)
                                                                                                                                gagaAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCgg             >  7:213473/1‑149 (MQ=255)
                                                                                                                                   aaagaaaTTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCg           >  5:59597/1‑148 (MQ=255)
                                                                                                                                   aaagaaaTTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTACTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCAAGGTTGATGAAGTACCTACTTTCGACCTACGCCTATTCGGCg           >  4:101984/1‑148 (MQ=255)
                                                                                                                                     agaaaTTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGaaa        >  3:59258/1‑149 (MQ=255)
                                                                                                                                      gaaaTTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAAt       >  5:60209/1‑149 (MQ=255)
                                                                                                                                         aTTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATtct    <  5:256022/149‑1 (MQ=255)
                                                                                                                                          ttGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATtctc   <  6:59597/149‑1 (MQ=255)
                                                                                                                                            gCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATTCTCa  <  5:239370/148‑1 (MQ=255)
                                                                                                                                            gCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATTCTCa  <  3:147357/148‑1 (MQ=255)
                                                                                                                                            gCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATTCTCa  <  4:59258/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCGATTTGCGAATATTTCGAATATAACTCGAACTCAGGCTGTTGCATCGTCGCCACAGCCCCGCAATCGATACCGCTCTGGAGATGCCCATGACCACCGCCAACCTCTCTAGCCATGTGCTTCCTGACGAGAAAGAAATTGCCGTGGCAGTCGAGTCAAGCCGCACGCTGGCAGCGTTCCTCTCGACAAAATCAGATACCCAACGGATCGAACTGCTCAACGAGGAAAATCAGCGCCAGGTTGTTGAAGTACCTACTTTCGCCCTACGCCTATTCGGCGAAATTCTCA  >  NC_002947/4690784‑4691071

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: ‑