Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,600 A→G 8.5% intergenic (+19/+19) 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,6000AG 8.5% 320.1 / 17.5 118intergenic (+19/+19)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (56/52);  new base G (6/4);  total (62/56)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 120 total aligned reads displayed.

TCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181451‑3181749
                                                                                                                                                     |                                                                                                                                                     
tCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCa                                                                                                                                                       >  1:1709367/1‑150 (MQ=255)
  ggCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACCGGATGGCAca                                                                                                                                                     >  2:64982/1‑149 (MQ=255)
      aCAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTTTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                 >  1:952910/1‑150 (MQ=255)
      aCAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                 <  2:399157/150‑1 (MQ=255)
        aGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAt                                                                                                                                               <  2:1289667/150‑1 (MQ=255)
        aGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAt                                                                                                                                               <  2:1296361/150‑1 (MQ=255)
         gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCACCCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCGCAGATTTTTGTTGAAGAAGCCACGATACCGGATGGCACTCGCCagc                                                                                                                                              >  2:1047130/1‑148 (MQ=255)
         gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATCGCACTCGCCATc                                                                                                                                              >  1:1637642/1‑150 (MQ=255)
         gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGACTCGGTTTTATCCACATTTTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATc                                                                                                                                              >  1:1631049/1‑150 (MQ=255)
          cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             <  2:1491138/150‑1 (MQ=255)
          cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             >  2:7821/1‑150 (MQ=255)
          cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             <  2:941872/150‑1 (MQ=255)
            atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCGACATATTTTTGTTGAATAATCCACGATACCGGATGGAACTCTCCATCCgg                                                                                                                                           >  2:870996/1‑150 (MQ=255)
            atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                           <  1:258729/150‑1 (MQ=255)
            atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                           <  2:1241009/150‑1 (MQ=255)
             tatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGt                                                                                                                                          <  1:1280101/150‑1 (MQ=255)
               tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTat                                                                                                                                        <  2:267095/150‑2 (MQ=255)
               tGCTGATTGGAAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCGCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTat                                                                                                                                        <  1:1101902/150‑2 (MQ=255)
                  tGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTg                                                                                                                                     <  2:1619891/150‑1 (MQ=255)
                  tGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTg                                                                                                                                     >  2:1716711/1‑150 (MQ=255)
                    aTTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTatcgtgg                                                                                                                                   <  1:543629/150‑7 (MQ=255)
                       ggTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGc                                                                                                                                <  2:1666777/150‑1 (MQ=255)
                        gTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                               >  1:50379/1‑150 (MQ=255)
                        gTAGCTGGCTGGTGAACGATCAACCCGCCCGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                               >  2:585948/1‑150 (MQ=255)
                            ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATc                                                                                                                           >  2:684133/1‑150 (MQ=255)
                              ggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTg                                                                                                                         <  1:1463459/150‑1 (MQ=255)
                               gctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGc                                                                                                                        >  1:1523144/1‑150 (MQ=255)
                               gctggGGAACGATCAACCCGCCGGATTTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGc                                                                                                                        <  2:527377/150‑1 (MQ=255)
                                ctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGAACACCAAGGATATGTCTCGGTTTTATCCACATATTTTTGATGAATAAGCCACGATACCGTATGGCACTCGCCATCCGGAAATTTTTAGCCTATCTGCa                                                                                                                       >  1:1291939/1‑150 (MQ=255)
                                 tggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAc                                                                                                                      <  2:1279443/150‑1 (MQ=255)
                                  ggTGAACGATCAACCCGGCGGAATTGTGATTCCTGAAGACAGTGCATTGATCATGCAGGTTATGTCTCGGTTTTATCAACATATTTTTGTTGAATAAGCCACGATACCGGATGGGACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                     <  1:758161/150‑1 (MQ=255)
                                  ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                     <  1:758155/150‑1 (MQ=255)
                                    tGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGCCCGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATATTTGTTGAATAAGCCGCGATACCGCATGGCACTCGCCATCCGGTTATTGTTAGCCTATCTGCACCGa                                                                                                                   >  1:275436/1‑150 (MQ=255)
                                    tGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGa                                                                                                                   >  2:844613/1‑150 (MQ=255)
                                        cGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGc                                                                                                               >  2:1377063/1‑150 (MQ=255)
                                         gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGAATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGAACCGACAGCa                                                                                                              >  1:1070155/1‑150 (MQ=255)
                                          aTCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCAt                                                                                                             <  1:237722/150‑1 (MQ=255)
                                            cAACCCCCCGTCATTCGCATTCGTGAAGACCGGGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATAc                                                                                                           <  2:243770/150‑1 (MQ=255)
                                               cccGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCa                                                                                                        <  1:7821/150‑1 (MQ=255)
                                               cccGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCa                                                                                                        >  2:503278/1‑150 (MQ=255)
                                                cggccgGAAGGTGCATTCGTGACGACCGTGCAGTGAGCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAg                                                                                                       <  1:278814/148‑1 (MQ=255)
                                                 cgccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAgg                                                                                                      >  1:1438187/1‑150 (MQ=255)
                                                 ccgcgGACTTGCCATTCGCGACCACCGTGCATAGAACACCCAGGATATGTCTCGCTTTTATCCACAAATTTTTGTTGAATAAGCCACGATACCGGATCGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAgg                                                                                                      <  2:1623528/147‑1 (MQ=255)
                                                  gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCATGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  2:1053021/150‑1 (MQ=255)
                                                          tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGTCTGCCCAtt                                                                                             >  1:100520/1‑150 (MQ=255)
                                                            gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  2:1172779/150‑1 (MQ=255)
                                                              aTTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCta                                                                                         <  2:1116441/150‑1 (MQ=255)
                                                                  gTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCTGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCTCATTTCTATAcc                                                                                     <  1:568527/150‑1 (MQ=255)
                                                                   tGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAccc                                                                                    <  1:1319937/150‑1 (MQ=255)
                                                                        aCCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAc                                                                               <  2:53401/150‑1 (MQ=255)
                                                                         ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                              <  1:976658/150‑1 (MQ=255)
                                                                         ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                              >  1:1287493/1‑150 (MQ=255)
                                                                           gTGCATTGATCACCCAGGATATGTCTCGGTTTTTTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCgg                                                                            <  2:952910/150‑1 (MQ=255)
                                                                            tGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGAAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGa                                                                           >  1:67220/1‑150 (MQ=255)
                                                                              cATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTTCCCATTTCTATACCCTCAACCCGAAt                                                                         >  2:292898/1‑150 (MQ=255)
                                                                              cATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAAt                                                                         >  2:247419/1‑150 (MQ=255)
                                                                                  gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                     >  2:792750/1‑150 (MQ=255)
                                                                                   aTCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTaa                                                                    >  2:542514/1‑150 (MQ=255)
                                                                                    tCAGCCATGAAATGTCTCGGTTTAATCCAGAATTTTTTGTTGAATAAGCCAAGATAACGGAGCGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTAAACCCTCAACCGGAATGGTAAt                                                                   <  2:325492/150‑1 (MQ=255)
                                                                                      aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 >  1:1470227/1‑150 (MQ=255)
                                                                                      aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 <  2:716860/150‑1 (MQ=255)
                                                                                       cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTAAACCCTCAACCGGAATGGTAATTgg                                                                <  2:50379/150‑1 (MQ=255)
                                                                                        ccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                               <  2:534484/150‑1 (MQ=255)
                                                                                           ggATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATAGCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTcc                                                            <  2:870089/150‑1 (MQ=255)
                                                                                                tGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGccc                                                       >  1:414643/1‑150 (MQ=255)
                                                                                                   ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                    >  2:603599/1‑150 (MQ=255)
                                                                                                    tcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATcc                                                   >  1:124033/1‑150 (MQ=255)
                                                                                                        ttttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACg                                               <  2:12171/150‑1 (MQ=255)
                                                                                                          ttATCCACATATTTTTGTTGAATAAGCCTCGATAGCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATAGGCTCCTGCCCATCCCACGCa                                             <  2:1637642/150‑1 (MQ=255)
                                                                                                          ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATAGGCTCCTGCCCATCCCACGCa                                             <  2:1631049/150‑1 (MQ=255)
                                                                                                              ccACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTa                                         <  1:684133/150‑1 (MQ=255)
                                                                                                               cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                        >  2:1549493/1‑150 (MQ=255)
                                                                                                               cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                        >  2:1549489/1‑150 (MQ=255)
                                                                                                                acaTATTTTTGTTGAATAATCCACGATACCGTATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTACTTGGCTCCTGCCCATCCCACGCACCTAac                                       <  2:1143334/150‑1 (MQ=255)
                                                                                                                acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  1:1514115/150‑1 (MQ=255)
                                                                                                                 caTTTTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAaca                                      <  2:1112926/150‑1 (MQ=255)
                                                                                                                 caTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAaca                                      <  2:1523144/150‑1 (MQ=255)
                                                                                                                 caTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAaca                                      >  2:433305/1‑150 (MQ=255)
                                                                                                                   tatTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGTCCATTTCTATACCCTCACCCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAAcaca                                    >  1:819160/1‑150 (MQ=255)
                                                                                                                    atTTTTTTATAATATGCCACGATAGCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACAt                                   <  1:870996/150‑1 (MQ=255)
                                                                                                                     tttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCAGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATa                                  <  1:666144/150‑1 (MQ=255)
                                                                                                                            tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           <  1:1716711/150‑1 (MQ=255)
                                                                                                                            tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           >  2:1312325/1‑150 (MQ=255)
                                                                                                                                tAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCa                       >  2:1215847/1‑150 (MQ=255)
                                                                                                                                  aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                     >  2:1104453/1‑150 (MQ=255)
                                                                                                                                  aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                     >  1:82771/1‑150 (MQ=255)
                                                                                                                                    ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCGCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaata                   >  1:1700993/1‑150 (MQ=255)
                                                                                                                                       cGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATg                <  1:1160100/150‑1 (MQ=255)
                                                                                                                                        attACCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAATTAATgg               >  1:411200/3‑149 (MQ=255)
                                                                                                                                         aTACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGct              <  2:1329099/150‑1 (MQ=255)
                                                                                                                                           aCCGGATGGCGAGTGCCATCAGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGTCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCGTGGCCATCCCACGCAACTAACACATACAACATCGGCTATTAATGctct            >  2:1608838/1‑150 (MQ=255)
                                                                                                                                            ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCATGGTTCCAATTTCTATACCCTCAACCGGAATGGTATTTGGCTCCTGCCCATCCCGCGCACCGAACCCCTACAACAACGGCAAATAATGCTCTg           >  2:666539/1‑150 (MQ=255)
                                                                                                                                            ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  2:1372105/1‑150 (MQ=255)
                                                                                                                                            ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACCACAACGGCAAATAATGCTCTg           >  2:1074046/1‑150 (MQ=255)
                                                                                                                                            ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  1:190298/1‑150 (MQ=255)
                                                                                                                                             cGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTgg          <  1:844613/150‑1 (MQ=255)
                                                                                                                                              ggATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGc         >  2:621783/1‑150 (MQ=255)
                                                                                                                                                 tGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCCCCTACCACATACAACAACGGCAAATAATGCTCTGGCGtt      >  2:1508134/1‑150 (MQ=255)
                                                                                                                                                   gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg    <  2:184018/150‑1 (MQ=255)
                                                                                                                                                     aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGACTGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt  <  1:1047130/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181451‑3181749

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: