Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,384 A→T 6.1% Q1298L (CAG→CTG)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,3840AT 6.1% 235.6 / 6.4 82Q1298L (CAG→CTG) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base T (3/2);  total (43/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 = 3.15e-01

GCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGG  >  NC_002947/1903235‑1903533
                                                                                                                                                     |                                                                                                                                                     
gCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATgcgca                                                                                                                                                       >  1:71868/1‑150 (MQ=255)
  aaGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCCGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCGGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATgcgcagc                                                                                                                                                     >  1:1595480/1‑150 (MQ=255)
  aaGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATgcgcagc                                                                                                                                                     >  1:775577/1‑150 (MQ=255)
           ggCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGACGTCCGATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCa                                                                                                                                            <  1:115433/150‑1 (MQ=255)
            gCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAg                                                                                                                                           <  1:109462/150‑1 (MQ=255)
             cACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAgg                                                                                                                                          >  1:1654599/1‑150 (MQ=255)
              aCCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAggg                                                                                                                                         >  2:784018/1‑150 (MQ=255)
                cTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGAt                                                                                                                                       >  2:1081319/1‑150 (MQ=255)
                 tGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTCGCATGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGAtg                                                                                                                                      >  1:1277978/1‑150 (MQ=255)
                  gCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCATGGAtgt                                                                                                                                     <  1:1418396/150‑1 (MQ=255)
                  gCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGAtgt                                                                                                                                     >  2:178788/1‑150 (MQ=255)
                    ggCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTgg                                                                                                                                   >  1:920916/1‑150 (MQ=255)
                       tGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATg                                                                                                                                >  1:1777557/1‑150 (MQ=255)
                           cacGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggc                                                                                                                            <  2:1654599/150‑1 (MQ=255)
                            acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggcg                                                                                                                           >  1:1523389/1‑150 (MQ=255)
                            acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggcg                                                                                                                           >  2:1414183/1‑150 (MQ=255)
                            acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCCATGCGCTGCGCATCGCCacacacat                                                                                                                                     >  2:966047/1‑133 (MQ=255)
                            acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCCATGCGCTGCGCATCGCCacacacat                                                                                                                                     <  1:966047/140‑8 (MQ=255)
                            acGGCGTTCTCGCGGGACCAGGTCGAGAAGAACTAAAAACAGCGGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggcg                                                                                                                           >  1:1314875/1‑150 (MQ=255)
                             cGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGt                                                                                                                          <  2:557903/150‑1 (MQ=255)
                              ggCGTTCTCGCGGGACCAGGTCGAGAAGATCGATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGtt                                                                                                                         >  1:348481/1‑150 (MQ=255)
                                   tctcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCggg                                                                                                                    >  1:1010494/1‑150 (MQ=255)
                                     tcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGTGGCGGGGGCGTATTTCTATGTGTGTGGAGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTg                                                                                                                  >  1:128252/1‑150 (MQ=255)
                                     tcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTg                                                                                                                  >  2:759392/1‑150 (MQ=255)
                                      tgcgGGTCCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGAGGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGa                                                                                                                 <  1:282113/149‑1 (MQ=255)
                                       gcgGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGAt                                                                                                                >  2:1648756/1‑150 (MQ=255)
                                           gACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTg                                                                                                              >  2:648028/1‑148 (MQ=255)
                                           gACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTg                                                                                                              <  1:648028/148‑1 (MQ=255)
                                             ccAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGAtt                                                                                                               >  2:1487203/1‑145 (MQ=255)
                                             ccAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGAtt                                                                                                               <  1:1487203/145‑1 (MQ=255)
                                             ccAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGATGCGCCGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGc                                                                                                          >  1:615417/1‑150 (MQ=255)
                                              cAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCACAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCg                                                                                                                      >  1:1276621/1‑137 (MQ=255)
                                              cAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCACAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCg                                                                                                                      <  2:1276621/137‑1 (MQ=255)
                                               aGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCTGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCg                                                                                                        >  1:1075629/1‑150 (MQ=255)
                                                ggTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGa                                                                                                       >  2:594748/1‑150 (MQ=255)
                                                  tCGAGAAGAGCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGc                                                                                                     <  2:615417/150‑1 (MQ=255)
                                                     agaagaTCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATg                                                                                                  <  1:178788/150‑1 (MQ=255)
                                                      gaagaTCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATgg                                                                                                 >  1:1494610/1‑150 (MQ=255)
                                                         gaTCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGcag                                                                                              >  1:177179/1‑148 (MQ=255)
                                                            cTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggat                                                                                           <  1:1862529/150‑1 (MQ=255)
                                                             tATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggatg                                                                                          <  2:1321446/150‑1 (MQ=255)
                                                              aTGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCGAGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCGAGGGAGGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggctgg                                                                                         >  1:1788163/1‑150 (MQ=255)
                                                               tgtgCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggatgga                                                                                        >  1:950765/1‑150 (MQ=255)
                                                                gtgCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggatggat                                                                                       >  2:1327351/1‑150 (MQ=255)
                                                                gtgCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggatggat                                                                                       <  1:220395/150‑1 (MQ=255)
                                                                                tCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCgg                                                                                                                     >  1:1604332/1‑104 (MQ=255)
                                                                                tCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCgg                                                                                                                     <  2:1604332/104‑1 (MQ=255)
                                                                                tCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCGGGATGGATGCGGcggcggcggcgg                                                                       >  2:1574976/1‑150 (MQ=255)
                                                                                 cGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGgcggctgcggc                                                                      <  2:775577/150‑1 (MQ=255)
                                                                                    gCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTa                                                                   <  2:1421133/150‑1 (MQ=255)
                                                                                     cAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTAt                                                                  <  1:1398301/150‑1 (MQ=255)
                                                                                        ggTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTc                                                               <  1:1648756/150‑1 (MQ=255)
                                                                                        ggTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTc                                                               <  1:784018/150‑1 (MQ=255)
                                                                                        ggTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTc                                                               >  2:1290288/1‑150 (MQ=255)
                                                                                           gcgcTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGTCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCCaa                                                            <  2:71868/150‑1 (MQ=255)
                                                                                           gcgcTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGaa                                                            <  2:1777557/150‑1 (MQ=255)
                                                                                                      gCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGATAGTCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCaa                                                 >  1:1367091/1‑150 (MQ=255)
                                                                                                       cAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAg                                                <  2:1494610/150‑1 (MQ=255)
                                                                                                        aaTGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAgg                                               <  1:594748/150‑1 (MQ=255)
                                                                                                           ggCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCATGAGCAAGGCCa                                            <  2:177179/150‑1 (MQ=255)
                                                                                                              tGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGc                                              <  2:965385/145‑1 (MQ=255)
                                                                                                              tGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGc                                              >  1:965385/1‑145 (MQ=255)
                                                                                                              tGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGc                                         <  2:565333/150‑1 (MQ=255)
                                                                                                                     gggggCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCg                                  <  1:1081319/150‑1 (MQ=255)
                                                                                                                        ggCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGATTGAGCAAGGCCCGGCGCTATCGgcg                               >  1:861776/1‑150 (MQ=255)
                                                                                                                            tATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGAt                           <  1:1677949/150‑1 (MQ=255)
                                                                                                                             aTTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGAtg                          >  1:779248/1‑150 (MQ=255)
                                                                                                                              tttCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGAtgt                         >  1:1839926/1‑150 (MQ=255)
                                                                                                                                tCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGTATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATgtgt                       <  2:236209/150‑1 (MQ=255)
                                                                                                                                 cTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTa                      <  2:1277978/150‑1 (MQ=255)
                                                                                                                                  tATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCCGGCGTTATCGGCGCGATGTGTAc                     <  2:950765/150‑1 (MQ=255)
                                                                                                                                  tATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTAc                     <  2:1075629/150‑1 (MQ=255)
                                                                                                                                  tATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCCGCGGCGGCTGCGCCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTAc                     <  2:1052372/150‑1 (MQ=255)
                                                                                                                                  tATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGAGTTGCGGGTGATTGTCGCGGGGCAGGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTAc                     >  2:1655208/1‑150 (MQ=255)
                                                                                                                                      tgtgtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGaa                 >  1:1016313/1‑150 (MQ=255)
                                                                                                                                      tgtgtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGaa                 <  1:644813/150‑1 (MQ=255)
                                                                                                                                      tgtgtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGaa                 <  2:1010494/150‑1 (MQ=255)
                                                                                                                                      tgtgtgGCGATGCGCAGCGCATGGCCAGGGAGGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGaa                 <  1:511283/150‑1 (MQ=255)
                                                                                                                                         gtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTg              >  1:143395/1‑150 (MQ=255)
                                                                                                                                         ctgGCCATGCGCTGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCAGGGGGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTg              >  2:1315768/2‑150 (MQ=255)
                                                                                                                                          tgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGGTGCGGCTGGGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGc             >  2:1402159/1‑150 (MQ=255)
                                                                                                                                           ggCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCt            >  2:362865/1‑150 (MQ=255)
                                                                                                                                            gCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTc           >  1:596356/1‑150 (MQ=255)
                                                                                                                                              caTGCGCTGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCtt                                                       <  1:1448719/103‑1 (MQ=255)
                                                                                                                                              caTGCGCTGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCtt                                                       >  2:1448719/2‑104 (MQ=255)
                                                                                                                                                   gcagcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATg    <  1:759392/150‑1 (MQ=255)
                                                                                                                                                   gcagcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGAACTAAAGTGCTCCAAGATg    >  2:1481401/1‑150 (MQ=255)
                                                                                                                                                     agcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATggg  >  1:1051344/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGG  >  NC_002947/1903235‑1903533

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