Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,601 C→A 8.5% intergenic (+20/+18) 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,6010CA 8.5% 349.3 / 17.0 118intergenic (+20/+18)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (55/53);  new base A (6/4);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 120 total aligned reads displayed.

GGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTT  >  NC_000913/3181453‑3181750
                                                                                                                                                    |                                                                                                                                                     
ggCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACCGGATGGCAca                                                                                                                                                      >  2:64982/1‑149 (MQ=255)
    aCAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                  <  2:399157/150‑1 (MQ=255)
     cAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCa                                                                                                                                                 >  1:779150/1‑150 (MQ=255)
      aGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAt                                                                                                                                                <  2:1289667/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)
       gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCAGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCACGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATc                                                                                                                                               >  2:1207155/1‑150 (MQ=255)
        cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                              <  2:941872/150‑1 (MQ=255)
        cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                              >  2:7821/1‑150 (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)
             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)
                      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)
                              ctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCa                                                                                                                        <  2:1598517/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                                                                                                         >  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)
                                                                                      ccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATTAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACACAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                                >  2:264487/1‑150 (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)
                                                                                                         tATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAc                                             <  1:351742/150‑1 (MQ=255)
                                                                                                         tATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACACTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAc                                             >  2:1663344/1‑150 (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)
                                                                                                                          tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                            <  1:1716711/150‑1 (MQ=255)
                                                                                                                          tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                            >  2:1312325/1‑150 (MQ=255)
                                                                                                                                aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                      >  1:82771/1‑150 (MQ=255)
                                                                                                                                  ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaata                    >  2:1104159/1‑150 (MQ=255)
                                                                                                                                  ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCGCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaata                    >  1:1700993/1‑150 (MQ=255)
                                                                                                                                     cGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATg                 <  1:1160100/150‑1 (MQ=255)
                                                                                                                                       aTACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGct               <  2:1329099/150‑1 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCATGGTTCCAATTTCTATACCCTCAACCGGAATGGTATTTGGCTCCTGCCCATCCCGCGCACCGAACCCCTACAACAACGGCAAATAATGCTCTg            >  2:666539/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:915524/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  2:140466/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  2:1372105/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACCACAACGGCAAATAATGCTCTg            >  2:1074046/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg            >  1:595450/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)
                                                                                                                                                    |                                                                                                                                                     
GGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTT  >  NC_000913/3181453‑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: