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

Only 100 of 125 total aligned reads displayed.

CGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181452‑3181749
                                                                                                                                                     |                                                                                                                                                    
cGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCAc                                                                                                                                                      >  1:765362/1‑150 (MQ=255)
cGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTGTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCAc                                                                                                                                                      >  2:110707/1‑150 (MQ=255)
 ggCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACt                                                                                                                                                     >  2:908404/1‑150 (MQ=255)
   cGACAGCTATATGCTGATTGGTAGCTTGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCATGAAGACCGTGCATTGATCACCCAGGATATGTCTCGTTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                   <  1:466125/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)
          tataTGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAAATGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCAACATATTTTTGTTGAATAAGAAACGATACCGGATGGCACTCGCCATGCg                                                                                                                                            >  1:1180637/1‑150 (MQ=255)
           atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                           <  2:648026/150‑1 (MQ=255)
           atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                           >  2:402452/1‑150 (MQ=255)
              tGCTGATTGGTAGCTGTCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTAATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTat                                                                                                                                        <  1:715393/150‑2 (MQ=255)
              tGCTGATTGGTAGCTGGCTGGTGAACTATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGAATTGATCACCCAGGTTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATACGGTat                                                                                                                                        <  2:545541/150‑2 (MQ=255)
                   aTTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGtt                                                                                                                                   <  2:1665139/150‑1 (MQ=255)
                       gTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                               <  1:538021/150‑1 (MQ=255)
                        tAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCt                                                                                                                              >  1:560122/1‑150 (MQ=255)
                         aGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTa                                                                                                                             >  1:17932/1‑150 (MQ=255)
                         aGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTa                                                                                                                             >  1:1311751/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)
                                tggTGAACGATCAACCCGCCGGAATTGCCATTCGTGAAGACCGTGCATTGATCCCACAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTAGCAGATCTGAAc                                                                                                                      >  2:306577/1‑150 (MQ=255)
                                        gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACATCa                                                                                                              >  1:80663/1‑150 (MQ=255)
                                        gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAACGACAGCa                                                                                                              >  2:343447/1‑150 (MQ=255)
                                           cAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATAc                                                                                                           <  2:67270/150‑1 (MQ=255)
                                             aCCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTc                                                                                                         <  2:1340303/150‑1 (MQ=255)
                                                cgccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAgg                                                                                                      <  2:1268995/150‑1 (MQ=255)
                                                 gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAAAAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                     <  1:1133826/150‑1 (MQ=255)
                                                 gccgGAATTGGCAATCGAGTAGACCGTGCATTGATCACCCAGGATATGTATCGGTTTAATCCACATATTTTTGTTGAATATGCCTCGATACCGGATGGCACTCGCCTTCGGGTAATTGTTAGCCTATATGCACCGACAGCATACTCAGGc                                                                                                     <  2:1408534/150‑1 (MQ=255)
                                                  ccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCt                                                                                                    <  2:14232/150‑1 (MQ=255)
                                                    ggAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGc                                                                                                  <  1:998865/150‑1 (MQ=255)
                                                     gAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGcc                                                                                                 <  2:1727486/150‑1 (MQ=255)
                                                     gAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGcc                                                                                                 <  1:1418275/150‑1 (MQ=255)
                                                        ttGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAt                                                                                              <  1:621410/150‑1 (MQ=255)
                                                         tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCGACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCTCCGACAGCATACTCAGGCGGCCCAtt                                                                                             >  1:204969/1‑150 (MQ=255)
                                                         tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACAGGATGGCACTCGCCATCAGGTAATTGTTAGCATATCTGCACAGAAAGCATACTCAGGCTGCCCAta                                                                                             >  1:72734/1‑149 (MQ=255)
                                                          ggCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAttt                                                                                            <  2:1569174/150‑1 (MQ=255)
                                                           gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           >  1:172670/1‑150 (MQ=255)
                                                           gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  1:1214764/150‑1 (MQ=255)
                                                           gCATACGTTAAGACCCTGCATTGAGCCCCCAGGTTATCTCTCGGTTAAATCCACATATTTTTGTTGAATAAGGCACGAGACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                           <  2:106858/150‑1 (MQ=255)
                                                              ttCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCtat                                                                                        >  1:430004/1‑150 (MQ=255)
                                                                     aGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCCTCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                 >  2:322624/1‑150 (MQ=255)
                                                                         cGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCg                                                                             >  2:1324164/1‑150 (MQ=255)
                                                                             cATTGATCACCCAGGATATGTCTCGGTTTTATCCACATAATTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAAt                                                                         <  2:1013910/150‑1 (MQ=255)
                                                                              aTTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATg                                                                        <  1:1130396/150‑1 (MQ=255)
                                                                                 gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                     >  2:1402217/1‑150 (MQ=255)
                                                                                 gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                     <  2:1274925/150‑1 (MQ=255)
                                                                                    cACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                  <  1:402452/150‑1 (MQ=255)
                                                                                    cACCAAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                  <  2:978145/150‑1 (MQ=255)
                                                                                     aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                 <  1:1093289/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)
                                                                                          ggATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTCGCtac                                                            >  1:1453890/1‑148 (MQ=255)
                                                                                             tatGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGc                                                         <  2:368021/150‑1 (MQ=255)
                                                                                             tatGTCTAGTTTTTATCCACATATTTTTGTTTAATAAGCCACGATACCGGATGGCACTCGGCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGc                                                         <  1:1423806/150‑1 (MQ=255)
                                                                                                gTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCa                                                      <  2:1536377/150‑1 (MQ=255)
                                                                                                gTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTACTGCCCa                                                      >  1:696275/1‑150 (MQ=255)
                                                                                                  ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                    <  1:798897/150‑1 (MQ=255)
                                                                                                     ggTTTTATCCTCATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCTGCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCTGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATccct                                                 <  2:1390660/150‑2 (MQ=255)
                                                                                                     ggTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCa                                                 <  1:1380307/150‑1 (MQ=255)
                                                                                                      tttttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCAc                                                <  2:578692/149‑1 (MQ=255)
                                                                                                         ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCTCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGTAATGGTAATTGGCTCCTGCCCTTCCCACGCa                                             <  2:1167112/150‑1 (MQ=255)
                                                                                                         ttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTATTTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAGTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGcc                                             >  1:956073/1‑149 (MQ=255)
                                                                                                           aTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAcc                                           >  2:521804/1‑150 (MQ=255)
                                                                                                            tCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCCCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCCTACTCACGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCCTCCCACGCGCCt                                          >  2:917193/1‑150 (MQ=255)
                                                                                                            tCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCt                                          >  2:692714/1‑150 (MQ=255)
                                                                                                               acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  2:1695104/150‑1 (MQ=255)
                                                                                                               acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  1:317318/150‑1 (MQ=255)
                                                                                                               acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                       <  1:133503/150‑1 (MQ=255)
                                                                                                                aatatTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAaca                                      <  1:1442142/149‑1 (MQ=255)
                                                                                                                     ttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATac                                 <  2:180249/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:560122/150‑1 (MQ=255)
                                                                                                                           tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                           <  2:80663/150‑1 (MQ=255)
                                                                                                                             aaTAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACgg                         <  2:1592794/150‑1 (MQ=255)
                                                                                                                                 aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCCCCTAACACATACAACAACGGCaaa                     >  2:674746/1‑150 (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)
                                                                                                                                      aattACCGGATGGCGAGTGCCATCCGGTATTTGTTAGCCTATCTGCACCGACCGCATACTCCGGCTGCCCAATTCTATACCCTCAACCGGAATGGTAAGTGGCTCCTGCCCATCCCACGCACCAAACACATACACCCAAGGCAaataaag                >  2:1325461/4‑148 (MQ=255)
                                                                                                                                        aTACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGct              >  2:442257/1‑150 (MQ=255)
                                                                                                                                          aCCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctct            >  1:1643183/1‑150 (MQ=255)
                                                                                                                                           ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCACCCCACGCACCTACCACCTACAACAACGGCAAATAATGCTCTg           >  1:1211576/1‑150 (MQ=255)
                                                                                                                                           ccGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg           >  1:717131/1‑150 (MQ=255)
                                                                                                                                           ccGGATGGCACTCGCCATCCGGTAATTGTGAGCCTAGCTGCACCGACAGCATACGCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTACCACATACAACAACGGCAAATAATGCTCTg           <  2:1062051/150‑1 (MQ=255)
                                                                                                                                            cGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTgg          <  1:1402217/150‑1 (MQ=255)
                                                                                                                                                tGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGACCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGtt      <  1:322624/150‑1 (MQ=255)
                                                                                                                                                  gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg    <  2:1276580/150‑1 (MQ=255)
                                                                                                                                                   cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   <  2:1601049/150‑1 (MQ=255)
                                                                                                                                                   cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   >  1:552437/1‑150 (MQ=255)
                                                                                                                                                   cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTggg   <  1:692714/150‑1 (MQ=255)
                                                                                                                                                    aCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt  <  2:69996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGT  >  NC_000913/3181452‑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: