Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,846,627 G→A 6.0% T38T (ACC→ACT PP_3395 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,846,6270GA 6.0% 192.9 / 2.8 66T38T (ACC→ACTPP_3395LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CGCAGCAGTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTG  >  NC_002947/3846481‑3846761
                                                                                                                                                  |                                                                                                                                      
cgcAGCAGTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTg                                                                                                                                     >  1:1105555/1‑150 (MQ=255)
cgcAGCAGTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTg                                                                                                                                     >  1:1105667/1‑150 (MQ=255)
 gcagcaGTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGc                                                                                                                                    >  2:4916/1‑150 (MQ=255)
  cagcagTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGcc                                                                                                                                   <  1:1160167/150‑1 (MQ=255)
      agTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGAGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCcttct                                                                                                                               >  2:358513/1‑150 (MQ=255)
          cAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGc                                                                                                                           >  2:345194/1‑150 (MQ=255)
           aaGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCa                                                                                                                          >  1:98470/1‑150 (MQ=255)
            aGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAg                                                                                                                         <  1:81869/150‑1 (MQ=255)
             gCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAga                                                                                                                        >  2:134024/1‑150 (MQ=255)
             gCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAga                                                                                                                        <  2:1144450/150‑1 (MQ=255)
              ccccGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCAATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAgac                                                                                                                       <  2:453609/150‑1 (MQ=255)
                         tGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACCGCACCATCGACGCGCGCTCGCATGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACa                                                                                                            <  2:437599/150‑1 (MQ=255)
                           ggCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAgg                                                                                                          >  2:894941/1‑150 (MQ=255)
                             cTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGtt                                                                                                        >  2:962650/1‑150 (MQ=255)
                                cACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAATGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCa                                                                                                                          <  1:412924/129‑1 (MQ=255)
                                cACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCa                                                                                                                          >  2:412924/1‑129 (MQ=255)
                                      gcgcTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAg                                                                                               <  2:367015/150‑1 (MQ=255)
                                       cgcTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAgg                                                                                              <  1:894941/150‑1 (MQ=255)
                                        gcTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAggt                                                                                             >  2:578353/1‑150 (MQ=255)
                                               aCGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAg                                                                                      <  2:308319/150‑1 (MQ=255)
                                                cGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAgg                                                                                     <  1:573308/150‑1 (MQ=255)
                                                 gTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGt                                                                                                                                      >  1:155787/1‑100 (MQ=255)
                                                 gTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGt                                                                                                                                      <  2:155787/100‑1 (MQ=255)
                                                 gTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAgg                                                                                     <  2:768939/149‑1 (MQ=255)
                                                 gTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAgg                                                                                     >  1:768939/1‑149 (MQ=255)
                                                  tttGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGGAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAgggg                                                                                   >  2:1281428/1‑150 (MQ=255)
                                                  tttGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAgggg                                                                                   >  2:474476/1‑150 (MQ=255)
                                                      aGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTc                                                                                         <  2:457499/140‑1 (MQ=255)
                                                      aGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTc                                                                                         >  1:457499/1‑140 (MQ=255)
                                                       gCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACTACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCa                                                                              >  2:1058062/1‑150 (MQ=255)
                                                       gCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCa                                                                              >  1:599499/1‑150 (MQ=255)
                                                         aaGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATc                                                                            <  2:98470/150‑1 (MQ=255)
                                                          aGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCg                                                                           <  1:962650/150‑1 (MQ=255)
                                                               gTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACTATGAACAGGTTCAGTGAAAGGTAGTCTATGGGTTGCATCGCTTgg                                                                      >  2:1030858/1‑150 (MQ=255)
                                                               gTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTgg                                                                      <  1:70240/150‑1 (MQ=255)
                                                               gTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGACGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTgg                                                                      <  1:851903/150‑1 (MQ=255)
                                                                tGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGCGGACGAGGGGATGCATCGGTTGGc                                                                     >  2:657442/1‑150 (MQ=255)
                                                                  ccAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAAc                                                                                                             <  1:509898/108‑1 (MQ=255)
                                                                  ccAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAAc                                                                                                             >  2:509898/1‑108 (MQ=255)
                                                                  ccAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTc                                                                   <  1:578353/150‑1 (MQ=255)
                                                                     aaGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACTACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGc                                                                <  1:1058062/150‑1 (MQ=255)
                                                                       gcgcgcCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCtt                                                              >  2:1068348/1‑150 (MQ=255)
                                                                         gcgcCTCAATGTCTGGAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTgg                                                            <  1:1281428/150‑1 (MQ=255)
                                                                               cAATGTCGGAAATACGTGTACTGACGGAGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAgg                                                      <  2:530405/150‑1 (MQ=255)
                                                                                    tCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGAcgacga                                                                                                                  >  2:63082/1‑85 (MQ=255)
                                                                                    tCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGAcgacga                                                                                                                  <  1:63082/85‑1 (MQ=255)
                                                                                         aaTACGTTTACTGACGGCGGACGGCGCCACGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTa                                            <  1:129410/150‑1 (MQ=255)
                                                                                                 tACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTa                                    <  1:134024/150‑1 (MQ=255)
                                                                                                 tACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTa                                    <  1:345194/150‑1 (MQ=255)
                                                                                                    tGACGGTGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCa                                 >  2:1026860/1‑150 (MQ=255)
                                                                                                     gacggcggacggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAg                                                       <  1:489600/127‑1 (MQ=255)
                                                                                                     gacggcggacggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAg                                                       >  2:489600/1‑127 (MQ=255)
                                                                                                           ggacggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCg                          <  1:4916/150‑1 (MQ=255)
                                                                                                           ggacggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGCGCCTTCTTCGGCGACGACGATCAACAGGTTCAGGGAAAGGTGGTCGAGGGGCTGCATCGCTTGGCTCGGCTTGGGGAAGGGAACGCCGTAGCACTGTAGACAAAACCg                          >  1:1124077/1‑150 (MQ=255)
                                                                                                             acggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGt                                     >  1:955278/1‑137 (MQ=255)
                                                                                                             acggcgCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGt                                     <  2:955278/137‑1 (MQ=255)
                                                                                                                         aCGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCa                        <  2:133740/138‑1 (MQ=255)
                                                                                                                         aCGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCa                        >  1:133740/1‑138 (MQ=255)
                                                                                                                              cgcTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCgg       <  1:909453/150‑1 (MQ=255)
                                                                                                                                 tCGCCTGCGCGGGCGATAGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGtt    <  1:81814/150‑1 (MQ=255)
                                                                                                                                 tCGCCTGCGCGGGCGATAGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGtt    <  1:81863/150‑1 (MQ=255)
                                                                                                                                 tCGCCTGCGCGGGCGATAGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGtt    >  2:81814/1‑150 (MQ=255)
                                                                                                                                 tCGCCTGCGCGGGCGATAGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGtt    >  2:81863/1‑150 (MQ=255)
                                                                                                                                  cGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCACAAACCGCATTGGCTTCTAAAACCGGGttt   <  1:358513/150‑1 (MQ=255)
                                                                                                                                  cGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGttt   <  1:1026860/150‑1 (MQ=255)
                                                                                                                                   gagcgcgcGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTg  <  2:504497/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
CGCAGCAGTGCAAGCCCCGCCGGGGTGGGCTCCACGCCGCGCTTGCTACGTTTGAGCAAGGCGGTGCCAAAGCGCGCCTCAATGTCTGAAATACGTTTACTGACGGCGGACGGCGCCATGAACGCGCGCTCGCCTGCGCGGGCGATGGTGCCTTCTTCGCAGACGACGATGAACAGGTTCAGGGAAAGGTGGTCGAGGGGTTGCATCGCTTGGCTCGGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTG  >  NC_002947/3846481‑3846761

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