Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,602 T→G 8.2% intergenic (+21/+17) 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,6020TG 8.2% 337.2 / 17.7 119intergenic (+21/+17)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (56/53);  new base G (6/4);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 119 total aligned reads displayed.

ACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTT  >  NC_000913/3181457‑3181750
                                                                                                                                                 |                                                                                                                                                    
aCAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTTTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                  >  1:952910/1‑150 (MQ=255)
 cAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCa                                                                                                                                                 >  1:779150/1‑150 (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)
   gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCAGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCACGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATc                                                                                                                                               >  2:1207155/1‑150 (MQ=255)
    cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                              <  2:941872/150‑1 (MQ=255)
    cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                              <  2:1491138/150‑1 (MQ=255)
    cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                              >  2:7821/1‑150 (MQ=255)
      atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCGACATATTTTTGTTGAATAATCCACGATACCGGATGGAACTCTCCATCCgg                                                                                                                                            >  2:870996/1‑150 (MQ=255)
      atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                            <  2:1241009/150‑1 (MQ=255)
      atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                            >  1:1361575/1‑150 (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:1716711/1‑150 (MQ=255)
            tGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTg                                                                                                                                      <  2:1619891/150‑1 (MQ=255)
              aTTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTatcgtgg                                                                                                                                    <  1:543629/150‑7 (MQ=255)
                 ggTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGc                                                                                                                                 <  2:1666777/150‑1 (MQ=255)
                  gTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                                >  1:12171/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)
                          ctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCa                                                                                                                        <  2:1598517/150‑1 (MQ=255)
                          ctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGAACACCAAGGATATGTCTCGGTTTTATCCACATATTTTTGATGAATAAGCCACGATACCGTATGGCACTCGCCATCCGGAAATTTTTAGCCTATCTGCa                                                                                                                        >  1:1291939/1‑150 (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)
                                           ccgcgGACTTGCCATTCGCGACCACCGTGCATAGAACACCCAGGATATGTCTCGCTTTTATCCACAAATTTTTGTTGAATAAGCCACGATACCGGATCGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAgg                                                                                                       <  2:1623528/147‑1 (MQ=255)
                                                    tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGTCTGCCCAtt                                                                                              >  1:100520/1‑150 (MQ=255)
                                                      gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                            <  2:1172779/150‑1 (MQ=255)
                                                       cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                           >  1:1575410/1‑150 (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:1287493/1‑150 (MQ=255)
                                                                   ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                               <  1:976658/150‑1 (MQ=255)
                                                                     gTGCATTGATCACCCAGGATATGTCTCGGTTTTTTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCgg                                                                             <  2:952910/150‑1 (MQ=255)
                                                                      tGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGAAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGa                                                                            >  1:67220/1‑150 (MQ=255)
                                                                        cATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTTCCCATTTCTATACCCTCAACCCGAAt                                                                          >  2:292898/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                                                                  <  2:716860/150‑1 (MQ=255)
                                                                                aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                  >  1:1470227/1‑150 (MQ=255)
                                                                                 cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTAAACCCTCAACCGGAATGGTAATTgg                                                                 <  2:50379/150‑1 (MQ=255)
                                                                                  ccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATTAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACACAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                                >  2:264487/1‑150 (MQ=255)
                                                                                     ggATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATAGCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTcc                                                             <  2:870089/150‑1 (MQ=255)
                                                                                          tGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGccc                                                        >  1:414643/1‑150 (MQ=255)
                                                                                             ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                     >  2:603599/1‑150 (MQ=255)
                                                                                              tcgcaattATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGCCACTCGCCATCCGGTAATTGTTAGCCTAACTGCACGGACAGCATACACACGCTGCCCATTTCAAAACCCTCAACCCGAATGGTAATTGGCTCCTGCCCATcc                                                    <  2:46320/144‑1 (MQ=255)
                                                                                                  ttttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACg                                                <  2:12171/150‑1 (MQ=255)
                                                                                                    ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATAGGCTCCTGCCCATCCCACGCa                                              <  2:1631049/150‑1 (MQ=255)
                                                                                                     tATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAc                                             <  1:351742/150‑1 (MQ=255)
                                                                                                     tATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACACTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAc                                             >  2:1663344/1‑150 (MQ=255)
                                                                                                         cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                         >  2:1549489/1‑150 (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)
                                                                                                                              ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaata                    >  2:1104159/1‑150 (MQ=255)
                                                                                                                              ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaata                    >  2:53677/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:140466/1‑150 (MQ=255)
                                                                                                                                      ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:915524/1‑150 (MQ=255)
                                                                                                                                      ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACCACAACGGCAAATAATGCTCTg            >  2:1074046/1‑150 (MQ=255)
                                                                                                                                      ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:595450/1‑150 (MQ=255)
                                                                                                                                      ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:190298/1‑150 (MQ=255)
                                                                                                                                      ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:184018/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)
                                                                                                                                            ttcACTGGACATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGTATGCCCATTTCTATACGCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAAGACATACAACAACGGCAAATAATGCTCTGGCGTTg      <  2:1291939/148‑1 (MQ=255)
                                                                                                                                             gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg     <  2:184018/150‑1 (MQ=255)
                                                                                                                                               aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGACTGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt   <  1:1047130/150‑1 (MQ=255)
                                                                                                                                                cTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGtt  >  2:358312/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTT  >  NC_000913/3181457‑3181750

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: