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

Only 100 of 123 total aligned reads displayed.

GGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181453‑3181749
                                                                                                                                                     |                                                                                                                                                   
ggCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACt                                                                                                                                                     >  2:908404/1‑150 (MQ=255)
  cGACAGCTATATGCTGATTGGTAGCTTGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCATGAAGACCGTGCATTGATCACCCAGGATATGTCTCGTTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                   <  1:466125/150‑1 (MQ=255)
   gACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGc                                                                                                                                                  <  2:53760/150‑1 (MQ=255)
    aCAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                 >  1:1268995/1‑150 (MQ=255)
     cAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCa                                                                                                                                                >  1:404186/1‑150 (MQ=255)
       gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATc                                                                                                                                              <  1:880380/150‑1 (MQ=255)
        cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCCACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCTCCCAGGCTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                             >  1:90834/1‑150 (MQ=255)
        cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTTATCACACAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACCGGATGGCACTCGCCATcc                                                                                                                                             >  2:134982/1‑150 (MQ=255)
          atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                           <  2:648026/150‑1 (MQ=255)
          atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                           >  2:402452/1‑150 (MQ=255)
            atGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTa                                                                                                                                         >  1:1295436/1‑150 (MQ=255)
            atGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCAATCCGGTa                                                                                                                                         >  2:353520/1‑150 (MQ=255)
             tGCTGATTGGTAGCTGTCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTAATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTat                                                                                                                                        <  1:715393/150‑2 (MQ=255)
             tGCTGATTGGTAGCTGGCTGGTGAACTATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGAATTGATCACCCAGGTTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATACGGTat                                                                                                                                        <  2:545541/150‑2 (MQ=255)
              gCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAAAAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAt                                                                                                                                       >  2:1133826/1‑150 (MQ=255)
                  aTTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGtt                                                                                                                                   <  2:1665139/150‑1 (MQ=255)
                      gTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                               <  1:538021/150‑1 (MQ=255)
                        aGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTa                                                                                                                             >  1:17932/1‑150 (MQ=255)
                         gctggctgGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTAt                                                                                                                            >  1:139011/1‑150 (MQ=255)
                           tggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCt                                                                                                                          <  1:180389/150‑1 (MQ=255)
                             gctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACCATACCGGATGGCACTCGCCATCAGGTAATTATTAGCCTATCTGc                                                                                                                        >  2:733081/1‑150 (MQ=255)
                               tggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAc                                                                                                                      <  1:716033/150‑1 (MQ=255)
                               tggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAc                                                                                                                      >  2:219210/1‑150 (MQ=255)
                               tggTGAACGATCAACCCGCCGGAATTGCCATTCGTGAAGACCGTGCATTGATCCCACAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTAGCAGATCTGAAc                                                                                                                      >  2:306577/1‑150 (MQ=255)
                                ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                     <  2:1447552/150‑1 (MQ=255)
                                ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                     >  1:180249/1‑150 (MQ=255)
                                 gTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCg                                                                                                                    >  1:1674296/1‑150 (MQ=255)
                                      cGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGc                                                                                                               >  1:1144375/1‑150 (MQ=255)
                                       gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACATCa                                                                                                              >  1:80663/1‑150 (MQ=255)
                                        aTCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCAt                                                                                                             <  2:1620917/150‑1 (MQ=255)
                                          cAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATAc                                                                                                           <  2:67270/150‑1 (MQ=255)
                                            aCCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTc                                                                                                         <  2:1340303/150‑1 (MQ=255)
                                                gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAAAAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  1:1133826/150‑1 (MQ=255)
                                                 ccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCt                                                                                                    <  2:14232/150‑1 (MQ=255)
                                                   ggAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGc                                                                                                  <  1:998865/150‑1 (MQ=255)
                                                    gAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGcc                                                                                                 <  1:1418275/150‑1 (MQ=255)
                                                     aaTTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGccc                                                                                                >  1:69996/1‑150 (MQ=255)
                                                      aTTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCa                                                                                               >  1:735779/1‑150 (MQ=255)
                                                       ttGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAt                                                                                              <  1:621410/150‑1 (MQ=255)
                                                        tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCGACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCTCCGACAGCATACTCAGGCGGCCCAtt                                                                                             >  1:204969/1‑150 (MQ=255)
                                                         ggCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAttt                                                                                            <  2:1569174/150‑1 (MQ=255)
                                                          gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  1:1214764/150‑1 (MQ=255)
                                                          gCATACGTTAAGACCCTGCATTGAGCCCCCAGGTTATCTCTCGGTTAAATCCACATATTTTTGTTGAATAAGGCACGAGACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  2:106858/150‑1 (MQ=255)
                                                              tCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCtata                                                                                       <  1:1391268/150‑1 (MQ=255)
                                                                gTGAAGACCGAGCATTGATCACCGAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCTCCGACAGCATACTCAGGCTGCCCATTTCTATAcc                                                                                     <  1:374472/150‑1 (MQ=255)
                                                                    aGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCCTCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                 >  2:322624/1‑150 (MQ=255)
                                                                        cGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCg                                                                             >  2:1324164/1‑150 (MQ=255)
                                                                             aTTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATg                                                                        <  1:1130396/150‑1 (MQ=255)
                                                                                gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                     <  2:1274925/150‑1 (MQ=255)
                                                                                gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                     >  2:1402217/1‑150 (MQ=255)
                                                                                   cACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                  <  1:402452/150‑1 (MQ=255)
                                                                                   cACCAAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                  <  2:978145/150‑1 (MQ=255)
                                                                                    aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 <  1:1093289/150‑1 (MQ=255)
                                                                                    aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCAAGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 <  1:614929/150‑1 (MQ=255)
                                                                                     cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTACTATACCCTCAACCGGAATGGTAATTgg                                                                <  1:1715872/150‑1 (MQ=255)
                                                                                      ccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                               >  2:317318/1‑150 (MQ=255)
                                                                                       cAGTATATGTCTCTGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGATGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCt                                                              <  1:110929/150‑1 (MQ=255)
                                                                                            tatGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGc                                                         <  2:368021/150‑1 (MQ=255)
                                                                                            tatGTCTAGTTTTTATCCACATATTTTTGTTTAATAAGCCACGATACCGGATGGCACTCGGCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGc                                                         <  1:1423806/150‑1 (MQ=255)
                                                                                               gTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCa                                                      <  2:1536377/150‑1 (MQ=255)
                                                                                                 ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                    <  1:798897/150‑1 (MQ=255)
                                                                                                    ggTTTTATCCTCATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCTGCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCTGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATccct                                                 <  2:1390660/150‑2 (MQ=255)
                                                                                                    ggTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCa                                                 <  1:1380307/150‑1 (MQ=255)
                                                                                                        ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCTCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGCCCTTCCCACGCa                                             <  2:1167112/150‑1 (MQ=255)
                                                                                                        ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTATTTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAGTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGcc                                             >  1:956073/1‑149 (MQ=255)
                                                                                                           tCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCCCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCCTACTCACGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCCTCCCACGCGCCt                                          >  2:917193/1‑150 (MQ=255)
                                                                                                           tCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCt                                          >  2:692714/1‑150 (MQ=255)
                                                                                                             cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTaa                                        >  2:603642/1‑150 (MQ=255)
                                                                                                              acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  1:133503/150‑1 (MQ=255)
                                                                                                              acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  2:1695104/150‑1 (MQ=255)
                                                                                                               aatatTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAaca                                      <  1:1442142/149‑1 (MQ=255)
                                                                                                                    ttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATac                                 <  2:180249/150‑1 (MQ=255)
                                                                                                                    ttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATac                                 <  2:139011/150‑1 (MQ=255)
                                                                                                                       tgttgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATacaac                              >  2:656090/1‑150 (MQ=255)
                                                                                                                        tttgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCTCATCCCACGCACCTAACACATacaaca                             <  1:563485/149‑1 (MQ=255)
                                                                                                                         ttgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATAcaacaa                            >  1:775189/1‑150 (MQ=255)
                                                                                                                          tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           >  1:139856/1‑150 (MQ=255)
                                                                                                                          tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           <  2:80663/150‑1 (MQ=255)
                                                                                                                          tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           <  2:560122/150‑1 (MQ=255)
                                                                                                                                aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCCCCTAACACATACAACAACGGCaaa                     >  2:674746/1‑150 (MQ=255)
                                                                                                                                aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                     <  2:1295436/150‑1 (MQ=255)
                                                                                                                                 gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCACCCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaat                    >  1:640383/1‑150 (MQ=255)
                                                                                                                                 gCCACGATACCGGATGGCACTCGCCATCCGGAAATTGTTAGCCTATCTGCACCGACCGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACCTACAACACCGGCaact                    >  1:158380/1‑148 (MQ=255)
                                                                                                                                    caattACCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaataat                 >  2:552583/5‑150 (MQ=255)
                                                                                                                                    aCGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaataat                 <  1:177181/150‑1 (MQ=255)
                                                                                                                                    aCGATACCGGATGGCACTAGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAaataat                 <  2:723129/150‑1 (MQ=255)
                                                                                                                                     aattACCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATg                >  2:1007692/4‑150 (MQ=255)
                                                                                                                                       aTACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGct              >  2:442257/1‑150 (MQ=255)
                                                                                                                                         aCCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctct            >  1:1643183/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  1:1437424/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCACCCCACGCACCTACCACCTACAACAACGGCAAATAATGCTCTg           >  1:1211576/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  1:494995/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  1:717131/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCACTCGCCATCCGGTAATTGTGAGCCTAGCTGCACCGACAGCATACGCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTACCACATACAACAACGGCAAATAATGCTCTg           <  2:1062051/150‑1 (MQ=255)
                                                                                                                                               tGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGACCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGtt      <  1:322624/150‑1 (MQ=255)
                                                                                                                                                 gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg    <  2:1276580/150‑1 (MQ=255)
                                                                                                                                                  cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   >  1:552437/1‑150 (MQ=255)
                                                                                                                                                  cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   <  2:1601049/150‑1 (MQ=255)
                                                                                                                                                  cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   <  1:692714/150‑1 (MQ=255)
                                                                                                                                                   aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt  <  2:69996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181453‑3181749

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