Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,601 C→A 6.6% 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 6.6% 393.0 / 13.6 137intergenic (+20/+18)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (71/56);  new base A (5/5);  total (76/61)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 137 total aligned reads displayed.

GCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181454‑3181749
                                                                                                                                                   |                                                                                                                                                    
gCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTc                                                                                                                                                    >  1:1755694/1‑150 (MQ=255)
 cGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACttg                                                                                                                                                   >  1:1869848/1‑148 (MQ=255)
 cGACAGCTATATGCTGATTGGTAGCTGGCGGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                   <  2:1013889/150‑1 (MQ=255)
  gACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGc                                                                                                                                                  <  1:867465/150‑1 (MQ=255)
       cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCACCCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCCCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             >  2:1196305/1‑150 (MQ=255)
       cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             >  1:1647399/1‑150 (MQ=255)
        tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCg                                                                                                                                            <  2:1015174/150‑1 (MQ=255)
        tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCg                                                                                                                                            <  2:308527/150‑1 (MQ=255)
         atatTCTGATGGGTAGCTGGCTGGTGAACGCTCAACCCGCCGGCATTGGCATTCGGGAAGACCGTGCAGTGATCACCCAGGATAGGGCGAGGGTTTATCCACATATTTTTGGTGAATACGCCACGAAACCGGATGGCGAGTGCCATCCgg                                                                                                                                           <  1:1644442/150‑1 (MQ=255)
         atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGTAATTGGCAGTCGTGTAGATCGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                           <  1:1059755/150‑1 (MQ=255)
         atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                           <  1:261356/150‑1 (MQ=255)
          tatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGt                                                                                                                                          >  1:1266812/1‑150 (MQ=255)
          tatGCTGATTGGTAGCTGGCTGGAGACCGATCAAACCGCCGGAATTGGCATTCGGGAAGACCGTTCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGt                                                                                                                                          <  1:886412/150‑1 (MQ=255)
            tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                        >  2:1805/1‑150 (MQ=255)
            tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                        >  2:1597730/1‑150 (MQ=255)
             gCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAt                                                                                                                                       >  1:1103591/1‑150 (MQ=255)
              cTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAtt                                                                                                                                      >  2:868593/1‑150 (MQ=255)
               ggATTGGTAGCTGGCTGGTGAACGATCACCCCGCCGGAATGGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTg                                                                                                                                     <  2:788228/149‑1 (MQ=255)
                  ttGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGCTCACCCAGGCTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCCCGATACCGGATGGCACTCGCCATCCGGTAATTGTTa                                                                                                                                  <  1:393300/150‑1 (MQ=255)
                  ttGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTatcgtggc                                                                                                                                  <  2:1111369/150‑8 (MQ=255)
                    ggTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGc                                                                                                                                >  1:700846/1‑150 (MQ=255)
                    ggTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTGTTATCCACATAGTTTTGTTGAATAAGCCACGATACCGGAGGGCACTCGCCGTCCGGTAATTGTTAGc                                                                                                                                >  2:1712749/1‑150 (MQ=255)
                     gTAGCTGGCTGGTTAACGATCAACCCGCCGGAATTCGCATTCGTGAAGACCGTGCAATGATCACCCATGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                               <  1:53267/150‑1 (MQ=255)
                        gctggctgGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTAt                                                                                                                            >  1:557745/1‑150 (MQ=255)
                        gctggctgGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATAACCCAGGATATGTCTCGGTATTATCCACATATTTTTGTTGAATAAGCAACGATACCGGATGGCTCTCTACATACGGTAATTGTTAGACTAt                                                                                                                            >  1:240310/1‑150 (MQ=255)
                         ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATc                                                                                                                           <  2:1236664/150‑1 (MQ=255)
                           ggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTATCCTATCTg                                                                                                                         >  1:1623803/1‑150 (MQ=255)
                            gctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGc                                                                                                                        >  1:622459/1‑150 (MQ=255)
                               ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCATATCTGCAcc                                                                                                                     >  1:373225/1‑150 (MQ=255)
                                  gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGAc                                                                                                                  >  2:1015534/1‑150 (MQ=255)
                                     cGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGc                                                                                                               >  2:1044096/1‑150 (MQ=255)
                                      gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACTGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCa                                                                                                              <  2:424662/150‑1 (MQ=255)
                                         cAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTAGTTGAATAAGCCACGATACCGGATGGCACTCGACATCCGGTAATTGTAAGCCTATATGCACCGACAGCATaa                                                                                                           >  2:414200/1‑149 (MQ=255)
                                             ccgccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAg                                                                                                       <  1:1326166/150‑1 (MQ=255)
                                             ccgccgGAATTGGCATTCGTGAAGACAGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATTTTTGTTGAATATGCCACGATACCGGAAGACACTCTCCATCCGGTAAATGTTAGCCTATCTAAACCGACATCATACTCAg                                                                                                       >  1:321047/1‑150 (MQ=255)
                                               gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGTATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  2:1579159/150‑1 (MQ=255)
                                               gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  1:793894/150‑1 (MQ=255)
                                               gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  1:912244/150‑1 (MQ=255)
                                                     tttGGCATTCGTGAAGAGCGTGCATTGATCACCCAGGAGATGTATCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGAGGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCa                                                                                               <  1:792674/149‑1 (MQ=255)
                                                         gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  2:753623/150‑1 (MQ=255)
                                                         gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  2:1546949/150‑1 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                          <  2:1045366/150‑1 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                          >  1:1282814/1‑150 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                          <  2:1744634/150‑1 (MQ=255)
                                                           aTTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCta                                                                                         >  1:1139399/1‑150 (MQ=255)
                                                             tCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCtata                                                                                       <  1:1232463/150‑1 (MQ=255)
                                                              cGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAc                                                                                      >  2:1108149/1‑150 (MQ=255)
                                                                tGAAGACCGTGCATTGATCTACCAGGATATGTCTCGGTTTTATACACTTATTTTTGTTGAATAAGCCACGATACCGTATGGCACTGGCCATCCGGTAATTGTTATCCTATCTGCACCGTCAGCATACTCAGGCTGCCGATTTCTATAccg                                                                                    <  2:1352101/150‑2 (MQ=255)
                                                                tGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAccc                                                                                    >  2:1423701/1‑150 (MQ=255)
                                                                   aGACCGTCGAATGATCACCCAGTATATGTCTCGGTTTTATCCACTTATTTTTGTTGAATTAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                 <  2:145079/150‑1 (MQ=255)
                                                                      ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                              >  2:832075/1‑150 (MQ=255)
                                                                       cGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCg                                                                             <  2:1488257/150‑1 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATTTTTTTGTTGAATAAGCCACGATTCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 >  1:1121416/1‑150 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 >  2:1178951/1‑150 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACCTGATGGCACTCGACATCCGGTAATTGTTAGCCTATCTGCACAGACAGCATACTCAGGCTGCCCATTTCTATACCCTAAACCGTAATGGTAATTg                                                                 >  1:1223500/1‑150 (MQ=255)
                                                                                    cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCCCTCGCAATCCGGTAATTGTTAGCCTATCTGCACCGACAGAATACTCAGGCTGCCCAATTCTATACCCTCAACCTGAATGGTAATTgg                                                                >  1:1342382/1‑150 (MQ=255)
                                                                                    cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATAGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTgg                                                                >  2:1292626/1‑150 (MQ=255)
                                                                                     ccAGGTTATGGCTCGGTTTTATCAACATATTTTTGTTGAATAATCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                               <  1:495414/150‑1 (MQ=255)
                                                                                     ccAGGATATGTCTCTGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                               <  1:562582/150‑1 (MQ=255)
                                                                                      cAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCCCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGAATGGTAATTGGCt                                                              >  2:1556420/1‑150 (MQ=255)
                                                                                         gATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCTCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCt                                                           >  1:427716/1‑150 (MQ=255)
                                                                                           tatGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTAGCATATCAGCACAGACAGCAGACTCAGGCTGCCCATTTCTATCCCCTAAACCGGAATGCTAATTGGATCCTGc                                                         >  2:1611591/1‑150 (MQ=255)
                                                                                           tatGTCTCGGTTGTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGAATGGTAATTGGCTCCTGc                                                         <  1:1436786/150‑1 (MQ=255)
                                                                                              gTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTCCCCa                                                      >  1:301186/1‑150 (MQ=255)
                                                                                                ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                    >  1:1084620/1‑150 (MQ=255)
                                                                                                ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTATATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                    >  1:1254538/1‑150 (MQ=255)
                                                                                                   ggTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCa                                                 <  2:1647399/150‑1 (MQ=255)
                                                                                                    gTTTTATCCACATTTTTTTGTTTAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGTCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCAc                                                <  2:275495/150‑1 (MQ=255)
                                                                                                    gTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCAc                                                <  1:1373724/150‑1 (MQ=255)
                                                                                                    gTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCAAc                                                <  1:824055/150‑1 (MQ=255)
                                                                                                      tttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGc                                              >  2:1865239/1‑150 (MQ=255)
                                                                                                             acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  1:511768/150‑1 (MQ=255)
                                                                                                                tatTTTTGTTGAATTAGCCACGATACCGGATGTCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCTGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAAcaca                                    <  2:1295442/150‑1 (MQ=255)
                                                                                                                tatTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCCCCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGACCATCCCACGCACCTAAcaca                                    >  1:1428058/1‑150 (MQ=255)
                                                                                                                   ttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATac                                 <  2:1755694/150‑1 (MQ=255)
                                                                                                                         tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           >  2:87706/1‑150 (MQ=255)
                                                                                                                         tgAATAAGCCACGATACCGGATGCCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           <  2:206986/150‑1 (MQ=255)
                                                                                                                          taaTAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                          <  2:1610075/149‑1 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                          >  2:1578985/1‑150 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                          >  2:853154/1‑150 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                          >  2:1842388/1‑150 (MQ=255)
                                                                                                                            aTAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCAGCACAACCCACGCACCTAACACATACAACAAAGGc                        >  2:1281415/1‑150 (MQ=255)
                                                                                                                               aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                     >  1:868658/1‑150 (MQ=255)
                                                                                                                               aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                     >  2:1319724/1‑150 (MQ=255)
                                                                                                                                gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaat                    >  1:1648907/1‑150 (MQ=255)
                                                                                                                                 ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACACCGGCAaata                   >  1:1014033/1‑150 (MQ=255)
                                                                                                                                     gATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGc               <  2:903171/150‑1 (MQ=255)
                                                                                                                                        aCCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctct            >  2:244835/1‑150 (MQ=255)
                                                                                                                                        aCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATTCCCTCTACCGGAATGGTAATTGGCTCCAGCCAAACCCACGCACCGAACACATACAACAACGGCAAAGAATGctct            >  2:1719412/1‑150 (MQ=255)
                                                                                                                                        aCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctct            >  2:1028513/1‑150 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACCCATACATCCACGGCAAATAATGCTCTg           >  1:687253/1‑150 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACATCAACGGCAAATAATGCTCTg           >  1:1690646/1‑150 (MQ=255)
                                                                                                                                          cGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTgg          <  1:1590303/150‑1 (MQ=255)
                                                                                                                                           ggATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGc         <  2:741748/150‑1 (MQ=255)
                                                                                                                                            gATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCg        <  2:1569914/150‑1 (MQ=255)
                                                                                                                                               ggCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTg     <  2:625873/150‑1 (MQ=255)
                                                                                                                                                gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGGTgg    >  1:851982/1‑150 (MQ=255)
                                                                                                                                                 cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   >  2:554783/1‑150 (MQ=255)
                                                                                                                                                 cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   >  1:179593/1‑150 (MQ=255)
                                                                                                                                                  aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt  >  1:1469158/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181454‑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: