Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,603 C→T 5.3% intergenic (+22/+16) ygiC → / ← ygiD ATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,6030CT 5.3% 397.1 / 5.5 129intergenic (+22/+16)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (64/58);  new base T (2/5);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 131 total aligned reads displayed.

CGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCG  >  NC_000913/3181455‑3181752
                                                                                                                                                    |                                                                                                                                                     
cGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                      <  2:585441/150‑1 (MQ=255)
cGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                      >  1:939677/1‑150 (MQ=255)
       tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCTCCCTGGATAGGTCTCGGTTTTATCCACATATATTTCTTGAATAAGCCACGATACCGCATGGGACTCGCCATCCg                                                                                                                                               >  2:351823/1‑150 (MQ=255)
       tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCg                                                                                                                                               >  1:1024876/1‑150 (MQ=255)
       tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGCTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCg                                                                                                                                               <  1:247362/150‑1 (MQ=255)
        atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                              <  2:625408/150‑1 (MQ=255)
           tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTAGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                           >  2:248570/1‑150 (MQ=255)
            gCTGATTGTTAGCTGCCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACAGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACAAAATATTCTTGAATAAGACACGATACCGGATGGCACTCGCCATCCGGTAAt                                                                                                                                          >  2:1438632/1‑150 (MQ=255)
            gCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTatc                                                                                                                                          <  2:1615131/150‑3 (MQ=255)
            gCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCAACATATTTTTGTTGAATAAGCAACGATACCGTATGGCACTCGCCATCCGGTAAt                                                                                                                                          >  2:286853/1‑150 (MQ=255)
               gATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGt                                                                                                                                       <  2:1281888/150‑1 (MQ=255)
                aTTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGtt                                                                                                                                      <  2:850456/150‑1 (MQ=255)
                   ggTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGCTCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGc                                                                                                                                   <  1:1636833/150‑1 (MQ=255)
                    gTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                                  <  2:68205/150‑1 (MQ=255)
                        ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTAta                                                                                                                              >  2:108348/1‑149 (MQ=255)
                        ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATc                                                                                                                              <  1:1665869/150‑1 (MQ=255)
                         tggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCt                                                                                                                             <  2:172997/150‑1 (MQ=255)
                          ggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTg                                                                                                                            >  1:80712/1‑150 (MQ=255)
                           gctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGc                                                                                                                           <  2:957773/150‑1 (MQ=255)
                             tggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAc                                                                                                                         >  2:416992/1‑150 (MQ=255)
                             tggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAc                                                                                                                         >  2:1523026/1‑150 (MQ=255)
                              ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATAACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                        <  2:71543/150‑1 (MQ=255)
                                 gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGAc                                                                                                                     >  2:207349/1‑150 (MQ=255)
                                  aaCGATCAACCCGCCTGAATTTGCATTCGTGAAGTCCGTTCATTGATCTCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACa                                                                                                                    <  1:1201875/150‑1 (MQ=255)
                                   aCGATCAACCCGCCGGAATTGACATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAg                                                                                                                   <  1:169815/150‑1 (MQ=255)
                                     gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCa                                                                                                                 <  1:207349/150‑1 (MQ=255)
                                     gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGACACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCa                                                                                                                 >  2:757742/1‑150 (MQ=255)
                                      aTCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCACCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCAt                                                                                                                <  2:1656539/150‑1 (MQ=255)
                                         aaCCCGCCGGAATTGGCATTCGTGAAGACCCTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACAAATTTTTGTAGAATAAGCCACGATACAGAATGGCACTCTACATCCGGTAATACTAAGACTATCTGCACAGACAGAATaat                                                                                                             >  1:127134/1‑148 (MQ=255)
                                           cccGCCGGAATTGGCATTCGTGAAGCCCGTGCATTGATCACCCAGGATATGACTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGCCACTCGCCATCCTGTAATTGTTAGCCTATCTGCACCTACAGCATACTCa                                                                                                           >  1:926507/1‑150 (MQ=255)
                                           cccGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCa                                                                                                           <  2:1000896/150‑1 (MQ=255)
                                            ccgccgGAATTGGCATTCGTGATGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAg                                                                                                          >  1:308853/1‑150 (MQ=255)
                                              gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                        <  1:486442/150‑1 (MQ=255)
                                              gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                        <  2:348503/150‑1 (MQ=255)
                                                 ggAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAATGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGc                                                                                                     >  2:1471821/1‑150 (MQ=255)
                                                       ggCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAttt                                                                                               >  1:924519/1‑150 (MQ=255)
                                                        gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                              <  1:249669/150‑1 (MQ=255)
                                                          aTTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCta                                                                                            <  2:905026/150‑1 (MQ=255)
                                                             cGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAc                                                                                         >  2:546020/1‑150 (MQ=255)
                                                                gAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCt                                                                                      <  2:939677/150‑1 (MQ=255)
                                                                  aGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                    >  2:501694/1‑150 (MQ=255)
                                                                     ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGCTGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                                 >  1:739075/1‑150 (MQ=255)
                                                                          cATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAAt                                                                            >  1:1180126/1‑150 (MQ=255)
                                                                          aaTTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAAt                                                                            <  2:1063576/149‑1 (MQ=255)
                                                                            ttGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATgg                                                                          <  1:1043452/150‑1 (MQ=255)
                                                                            ttGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATgg                                                                          <  2:1506648/150‑1 (MQ=255)
                                                                            ttGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATgg                                                                          <  2:522950/150‑1 (MQ=255)
                                                                             tGATCACCCAGTATATGTCTGGGTTTTATCCACATATTTTTGGTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGt                                                                         <  2:80712/150‑1 (MQ=255)
                                                                             tGATCACCCAGGATATGTCTCGGTTTTATCCATATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGt                                                                         >  2:792122/1‑150 (MQ=255)
                                                                             tGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGt                                                                         >  1:503512/1‑150 (MQ=255)
                                                                             tGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGt                                                                         <  1:1741266/150‑1 (MQ=255)
                                                                                 cACCCAGGATATGTCTCGGTTTTATCCACATTTTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                     <  1:810991/150‑1 (MQ=255)
                                                                                  aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                    >  2:1756579/1‑150 (MQ=255)
                                                                                  aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                    <  1:1704552/150‑1 (MQ=255)
                                                                                  aCCCAGGATATGTCTCGGGTTTATCCAAATATTTTTGGTTAATCAGCCTCGATAGCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCATCCGGAATGGTAATTg                                                                    <  2:1127618/150‑1 (MQ=255)
                                                                                     cAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCt                                                                 <  2:16089/150‑1 (MQ=255)
                                                                                       ggATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTcc                                                               <  1:1816210/150‑1 (MQ=255)
                                                                                           ttGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATTCTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGcc                                                           <  2:1036439/149‑1 (MQ=255)
                                                                                           atGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTCTCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGcc                                                           <  2:1036435/150‑1 (MQ=255)
                                                                                              tctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCAt                                                        >  1:1237379/1‑150 (MQ=255)
                                                                                              tctcGGATTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATAGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGATAATTGGCTCCTGCCCAt                                                        >  1:1552069/1‑150 (MQ=255)
                                                                                                tcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATcc                                                      <  1:416992/150‑1 (MQ=255)
                                                                                                acGGTTTTATCCACATATTTTTGTTGAATACGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATcc                                                      <  2:926434/149‑1 (MQ=255)
                                                                                                 cGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCACCCGGAATGGTAATTGGCTCCTGCCCATccc                                                     <  1:721494/150‑1 (MQ=255)
                                                                                                 cGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATccc                                                     >  1:1731646/1‑150 (MQ=255)
                                                                                                 cGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCGGCCCATccc                                                     >  2:1707420/1‑150 (MQ=255)
                                                                                                  ggTTTTATCCACATATTTATGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCa                                                    <  1:217840/150‑1 (MQ=255)
                                                                                                      ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCa                                                <  1:608569/150‑1 (MQ=255)
                                                                                                         tCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAGCCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAcca                                             >  1:620405/1‑149 (MQ=255)
                                                                                                         accACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATACGGTAATTGTTAGCCTATCTGAACCTACAGCATACTCAGTCTGCCCATTTCTATACCCTCCTCCGGAATGGTAAATGGATCCTGCCCATCCAACGCAcca                                             >  1:620408/2‑149 (MQ=255)
                                                                                                           cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCACACGCACCTaa                                           >  1:795567/1‑150 (MQ=255)
                                                                                                           cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACACTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                           >  1:1360303/1‑150 (MQ=255)
                                                                                                           cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACACTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                           >  1:725057/1‑150 (MQ=255)
                                                                                                            acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGACTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACACTCAACCGGAATGGTAATTCGCTCCTGCCAATCCCACGCACCTacc                                          >  1:50944/1‑148 (MQ=255)
                                                                                                            acaTAATTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCCTTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                          <  2:1192405/150‑1 (MQ=255)
                                                                                                             aatttttttGTTGAATAAGCCACGATACCGGATGTCACTTGCCATGCGGTAATTGTTAGCCTATCTGCACCGACATCATACTCAGGCTGCCCATTTCTATACCCTCTACCGGATTGGTAATTGGCTCCTGCCCATCCCACGGACCTAaca                                         <  1:1619663/146‑1 (MQ=255)
                                                                                                                atTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACAt                                      >  2:1181777/1‑150 (MQ=255)
                                                                                                                   tttGTTGAATAAGCCACGATAACGGATTGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGATGCCCATTTCTATACACTCAACCGTAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATaca                                   <  1:286853/150‑1 (MQ=255)
                                                                                                                    ttgttgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATacaa                                  <  1:1169303/150‑1 (MQ=255)
                                                                                                                        tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCCCCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCCCCTAACACATACaacaac                              >  1:1396616/1‑150 (MQ=255)
                                                                                                                           aTAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGc                           >  2:636778/1‑150 (MQ=255)
                                                                                                                            tAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACTCATACAACAACGGCa                          >  1:1671092/1‑150 (MQ=255)
                                                                                                                            tAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCa                          <  2:523583/150‑1 (MQ=255)
                                                                                                                              aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                        >  1:177310/1‑150 (MQ=255)
                                                                                                                                ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATTCTCAGGCTGCCCAGTTCTATACCCTCACCCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAGCACATACCACAACGGCAaata                      >  1:84954/1‑150 (MQ=255)
                                                                                                                                 cACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaataa                     >  1:215985/1‑150 (MQ=255)
                                                                                                                                 cACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaataa                     >  2:627758/1‑150 (MQ=255)
                                                                                                                                  aCGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCCCCTAACACATACAACAACGCCAaataat                    >  2:467265/1‑150 (MQ=255)
                                                                                                                                   cGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATg                   <  1:1449602/150‑1 (MQ=255)
                                                                                                                                    taTACTGGTTGGCAATCGGCATCCGGTAATTTTTTTACTATCTGCACCGTATGAATACTCAGGCTGGCCATTTCTATAACATCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAAAGGCAATTAATGc                  <  2:127134/149‑1 (MQ=255)
                                                                                                                                      tACCGGATGGGACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctc                <  2:209102/150‑1 (MQ=255)
                                                                                                                                        ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACATCAACGGCAAATAATGCTCTg              >  1:1671997/1‑150 (MQ=255)
                                                                                                                                        ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGATTGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAAAGCCTAATAATGCTCTg              >  2:350013/1‑150 (MQ=255)
                                                                                                                                         cGGATGGCACTCGCCATCCTGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATATCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTgg             <  1:359684/150‑1 (MQ=255)
                                                                                                                                          ggATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCAATACCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGc            >  1:315689/1‑150 (MQ=255)
                                                                                                                                           gATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCg           <  2:1537806/150‑1 (MQ=255)
                                                                                                                                                 aatCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACATCATACTCAGGCTGACCATTTCTATACCCTCAACCGGATTGGTAATTTGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGGTCTGGCGTTgggg     <  2:926507/148‑2 (MQ=255)
                                                                                                                                                 aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt     >  2:1461216/1‑150 (MQ=255)
                                                                                                                                                    cGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCg  <  2:974592/150‑1 (MQ=255)
                                                                                                                                                    cGCCATCCGGTAATTGTTAGCCTATCTGAACCGACAGCATACTCAGTCTGACCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGTCGTTGGGTTCg  <  1:39081/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCG  >  NC_000913/3181455‑3181752

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