Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,080 G→T 5.6% intergenic (+16/‑182) curA → / → mcbR NADPH‑dependent curcumin/dihydrocurcumin reductase/DNA‑binding transcriptional dual regulator McbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,0800GT 5.6% 508.7 / 11.4 150intergenic (+16/‑182)curA/mcbRNADPH‑dependent curcumin/dihydrocurcumin reductase/DNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base G (67/74);  new base T (3/6);  total (70/80)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 156 total aligned reads displayed.

GTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAG  >  NC_000913/1519933‑1520226
                                                                                                                                                   |                                                                                                                                                  
gTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTa                                                                                                                                                  >  1:1565576/1‑150 (MQ=255)
atGAAAGAGCATAAAATCCACAACAGAGAAGAAATTACTGACAGTTTAGAGAATGATCCAAAAAAAAAAAACGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTa                                                                                                                                                  >  2:1505998/2‑150 (MQ=255)
  gAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                                                                >  1:554288/1‑150 (MQ=255)
  gAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                                                                <  1:1508511/150‑1 (MQ=255)
    aaGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTACgcc                                                                                                                                              >  1:120463/1‑150 (MQ=255)
     agagGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                                                                             >  1:1821920/1‑150 (MQ=255)
      gagGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                                                                            >  1:1437692/1‑150 (MQ=255)
       agGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgcc                                                                                                                                           <  1:1358671/150‑1 (MQ=255)
       agGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                                                                           <  2:1895966/150‑1 (MQ=255)
        ggATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                                                                          >  2:1332447/1‑150 (MQ=255)
        ggATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                                                                          >  1:295988/1‑150 (MQ=255)
         gATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                                                                         <  1:760745/150‑3 (MQ=255)
          aTAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                                                                        <  1:2084244/150‑4 (MQ=255)
          aTAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                                                                        >  1:30736/1‑150 (MQ=255)
           tAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                                                                       <  1:383803/150‑5 (MQ=255)
           tAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                                                                       <  1:659549/150‑5 (MQ=255)
           tAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaa                                                                                                                                       <  1:1664941/150‑1 (MQ=255)
            aaaaTCCACTACCGAGAATAAATTACTGACGGTTTAGAGAATGCGCCACAGAAGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCTCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttga                                                                                                                                      <  2:1747625/150‑6 (MQ=255)
               aTCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAt                                                                                                                                   <  1:2022717/150‑1 (MQ=255)
                tCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATc                                                                                                                                  >  1:199980/1‑150 (MQ=255)
                   aCTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                                                               >  2:1934670/1‑150 (MQ=255)
                   aCTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                                                               <  1:468213/150‑1 (MQ=255)
                               aaaTTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTa                                                                                                                   <  1:2392796/150‑1 (MQ=255)
                                aaTTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAt                                                                                                                  <  1:2329264/150‑1 (MQ=255)
                                     cTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGCTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                                                             >  2:1781823/1‑150 (MQ=255)
                                     cTGACGGTGTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                                                             >  1:2038311/1‑150 (MQ=255)
                                        aCGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCt                                                                                                          <  2:1181018/150‑1 (MQ=255)
                                            tttAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaa                                                                                                      >  1:2296483/1‑150 (MQ=255)
                                              tAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGgcggcg                                                                                                                                                    >  1:1895301/1‑102 (MQ=255)
                                              tAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                                                                    <  1:1777560/150‑1 (MQ=255)
                                                 agaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTTGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                                                                 <  1:1163649/150‑1 (MQ=255)
                                                 agaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCCGCGTTAATtaac                                                                                                 >  2:1777073/1‑150 (MQ=255)
                                                 agaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                                                                 <  2:1159309/150‑1 (MQ=255)
                                                      gcgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                                                                            <  2:1594415/150‑1 (MQ=255)
                                                      gcgcCACAGACGGTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                                                                            <  2:2357199/150‑1 (MQ=255)
                                                         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                                                                         >  1:1957082/1‑150 (MQ=255)
                                                         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                                                                         >  1:170459/1‑150 (MQ=255)
                                                            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                                                                      >  2:2164440/1‑150 (MQ=255)
                                                               aCGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTg                                                                                   <  1:75362/150‑1 (MQ=255)
                                                               aCGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTg                                                                                   <  2:2380317/150‑1 (MQ=255)
                                                                  tttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTc                                                                                <  2:1780594/150‑1 (MQ=255)
                                                                    tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                                                                              <  2:856037/150‑1 (MQ=255)
                                                                    tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                                                                              >  2:1357107/1‑150 (MQ=255)
                                                                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgctg                                                                           <  1:183689/150‑1 (MQ=255)
                                                                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgctg                                                                           <  2:565521/150‑1 (MQ=255)
                                                                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGGGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCACGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgctg                                                                           <  2:546715/150‑1 (MQ=255)
                                                                        ggCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACTGCGGCGTAAGCCTACATAAAAATCCAGGCAAAAAGCATATTTTCAGCTTTAATTAACTAACCACATCATTGATctgctga                                                                          >  2:273196/1‑149 (MQ=255)
                                                                             gctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                                                                     <  2:621138/150‑1 (MQ=255)
                                                                              ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                                                                    >  1:1180070/1‑150 (MQ=255)
                                                                               tgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaac                                                                   >  2:1218570/1‑150 (MQ=255)
                                                                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa                                                             >  2:1197267/1‑150 (MQ=255)
                                                                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa                                                             >  2:141299/1‑150 (MQ=255)
                                                                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                                                            <  1:459154/150‑1 (MQ=255)
                                                                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                                                            <  2:1004960/150‑1 (MQ=255)
                                                                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAATCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg                                                           <  1:1033131/150‑1 (MQ=255)
                                                                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg                                                           <  2:1521519/150‑1 (MQ=255)
                                                                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg                                                           >  1:2188114/1‑150 (MQ=255)
                                                                                        aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt                                                          <  2:1803525/150‑1 (MQ=255)
                                                                                           aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                       <  2:1647060/150‑1 (MQ=255)
                                                                                           aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                       >  1:2053529/1‑150 (MQ=255)
                                                                                           aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                       <  1:696794/150‑1 (MQ=255)
                                                                                           aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTATGCTGCTAATAACAACAAAGGTGAg                                                       >  2:2152755/1‑150 (MQ=255)
                                                                                           aCTTAGTCAAAGTGGTGATCCTCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGAGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTTCTGCTAATAACAACAAAGGTGAg                                                       <  1:1595474/150‑1 (MQ=255)
                                                                                            cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGt                                                      >  1:1000341/1‑150 (MQ=255)
                                                                                            cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTg                                                     <  2:2021482/150‑1 (MQ=255)
                                                                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGGTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt                                                  >  2:764308/1‑150 (MQ=255)
                                                                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGttt                                                 >  2:295252/1‑150 (MQ=255)
                                                                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGttt                                                 >  1:1115356/1‑150 (MQ=255)
                                                                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg                                                           >  2:2292159/1‑138 (MQ=255)
                                                                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt                                               <  1:34535/150‑1 (MQ=255)
                                                                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt                                               <  1:969628/150‑1 (MQ=255)
                                                                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt                                               <  1:1537656/150‑1 (MQ=255)
                                                                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt                                               <  1:1943889/150‑1 (MQ=255)
                                                                                                    aaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATg                                              <  2:735097/150‑1 (MQ=255)
                                                                                                    aaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATg                                              <  2:1565576/150‑1 (MQ=255)
                                                                                                    aaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATg                                              >  2:187992/1‑150 (MQ=255)
                                                                                                     aGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGa                                             <  1:1457864/150‑1 (MQ=255)
                                                                                                      gtggtgACCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGat                                            >  1:2084034/1‑150 (MQ=255)
                                                                                                         gtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattat                                         >  1:2373065/1‑150 (MQ=255)
                                                                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTc                                       >  1:1616083/1‑150 (MQ=255)
                                                                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTc                                       >  1:1074868/1‑150 (MQ=255)
                                                                                                             tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCAc                                     <  2:1181553/150‑1 (MQ=255)
                                                                                                             tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCAc                                     >  1:1003730/1‑150 (MQ=255)
                                                                                                              ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACt                                    <  2:462462/150‑1 (MQ=255)
                                                                                                                    ggCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATa                              <  1:1651223/150‑1 (MQ=255)
                                                                                                                       gggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATacac                           <  2:1432313/150‑1 (MQ=255)
                                                                                                                        ggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACAcc                          <  1:149673/150‑1 (MQ=255)
                                                                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAgg                       <  2:1821920/150‑1 (MQ=255)
                                                                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAgg                       >  2:1035099/1‑150 (MQ=255)
                                                                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAgg                       >  2:51327/1‑150 (MQ=255)
                                                                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGa                    >  1:538492/1‑150 (MQ=255)
                                                                                                                                  aaaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTc                <  2:1115356/150‑1 (MQ=255)
                                                                                                                                   aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTcc               <  2:1853047/150‑1 (MQ=255)
                                                                                                                                   aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTcc               >  1:884234/1‑150 (MQ=255)
                                                                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGAGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAAGTCCtt             >  2:1116228/1‑150 (MQ=255)
                                                                                                                                       ttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCt           >  1:714012/1‑150 (MQ=255)
                                                                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGc         >  2:2393981/1‑150 (MQ=255)
                                                                                                                                          aaCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGcc        <  1:1126545/150‑1 (MQ=255)
                                                                                                                                            cggcggCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCAt      >  2:1447929/1‑150 (MQ=255)
                                                                                                                                             ggcggcGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATg     <  1:295252/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAG  >  NC_000913/1519933‑1520226

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