Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,738 A→G 6.1% R282R (CGT→CGC)  cysM ← cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,538,7380AG 6.1% 241.7 / 5.3 82R282R (CGT→CGC) cysMcysteine synthase B
Reads supporting (aligned to +/- strand):  ref base A (38/39);  new base G (2/3);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATG  >  NC_000913/2538589‑2538887
                                                                                                                                                     |                                                                                                                                                     
tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCa                                                                                                                                                       >  1:1045376/1‑150 (MQ=255)
    gCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGAt                                                                                                                                                   >  1:163816/1‑150 (MQ=255)
    gCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGAt                                                                                                                                                   <  1:1522268/150‑1 (MQ=255)
     cAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATc                                                                                                                                                  <  1:1461100/150‑1 (MQ=255)
      aGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCg                                                                                                                                                 <  1:728560/150‑1 (MQ=255)
       gAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGc                                                                                                                                                <  1:1531523/150‑1 (MQ=255)
         aaCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGcta                                                                                                                                              <  1:543252/150‑3 (MQ=255)
         aaCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAg                                                                                                                                              <  2:1253703/150‑1 (MQ=255)
          aCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAga                                                                                                                                             >  1:670500/1‑150 (MQ=255)
           ccTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGctacc                                                                                                                                            <  1:1004039/150‑5 (MQ=255)
           ccTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGTGGCGATCGctacc                                                                                                                                            <  1:575781/150‑5 (MQ=255)
                tAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCg                                                                                                                                       <  1:515143/150‑1 (MQ=255)
                  ggCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGcc                                                                                                                                     >  1:862526/1‑150 (MQ=255)
                  ggCCTGAAAATCCTAGAGCATCAGGCAGTTTTGCGTTTGTCATAAGTCTCCGATGCTATTAATCCTTAAATCCCCCCCCCCTGGCTAAAATTCACTTACCCAAACACACCCGTTGAAAGCTAGAGATACCCAACATCGCAGATGATCGcc                                                                                                                                     >  2:827815/1‑150 (MQ=255)
                   gCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGcca                                                                                                                                    >  1:701853/1‑150 (MQ=255)
                    ccTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccac                                                                                                                                   >  2:1032187/1‑150 (MQ=255)
                    ccTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccac                                                                                                                                   <  1:360480/150‑1 (MQ=255)
                     cTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccacc                                                                                                                                  >  2:211267/1‑150 (MQ=255)
                      tGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccacca                                                                                                                                 <  1:896192/150‑1 (MQ=255)
                         tAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccaccacca                                                                                                                              <  1:1061469/150‑1 (MQ=255)
                            gCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGCTCGCCACGATCGCAGATGATCGCCACCACCACcg                                                                                                                           >  2:1289862/1‑150 (MQ=255)
                            gCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACcg                                                                                                                           >  1:87996/1‑150 (MQ=255)
                             cGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACcgc                                                                                                                          >  2:1575510/1‑150 (MQ=255)
                                aGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTc                                                                                                                       <  1:1054140/150‑1 (MQ=255)
                                    cATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAggg                                                                                                                   <  2:518622/150‑1 (MQ=255)
                                      tCAGGCAGTTTTGCGTTTTTCATAAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACACCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGtt                                                                                                                 >  2:1095488/1‑150 (MQ=255)
                                       cAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGCTCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGtaa                                                                                                                >  2:1538159/1‑148 (MQ=255)
                                         ggCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGc                                                                                                              <  1:84541/150‑1 (MQ=255)
                                              ttttGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTg                                                                                                         <  2:318601/150‑1 (MQ=255)
                                               tttGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGc                                                                                                        <  1:827815/150‑1 (MQ=255)
                                                       tGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCACTTCCCCAAACCCCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCa                                                                                                >  2:1410199/1‑150 (MQ=255)
                                                          catcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTg                                                                                             <  2:1020768/150‑1 (MQ=255)
                                                          catcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTg                                                                                             <  2:1081218/150‑1 (MQ=255)
                                                            tcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCACTTACCCAAACACCCCGGTAGAAAGGTAGCGATCCCCACGATCGCAGATGATCGCCACAACCACCGCGTCAGCGATAGCTTTTGCCACCCGCAGTGCt                                                                                           >  2:172121/1‑150 (MQ=255)
                                                                tctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCgg                                                                                       <  1:1284834/150‑1 (MQ=255)
                                                                tctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCgg                                                                                       <  2:701853/150‑1 (MQ=255)
                                                                 ctcCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGc                                                                                      <  2:1237201/150‑1 (MQ=255)
                                                                   ccGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCaa                                                                                    >  2:1041449/1‑150 (MQ=255)
                                                                    cGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCTCCACGATCGCAGATGATCGCCACCACCACCGCCTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAAc                                                                                   >  1:198501/1‑150 (MQ=255)
                                                                     gATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAAcc                                                                                  <  2:494220/150‑1 (MQ=255)
                                                                       tGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACcgc                                                                                <  1:1272765/150‑1 (MQ=255)
                                                                          tATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCCCCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgcc                                                                             >  2:1454660/1‑150 (MQ=255)
                                                                           aTTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgccg                                                                            >  1:392684/1‑150 (MQ=255)
                                                                           aTTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgccg                                                                            >  2:135770/1‑150 (MQ=255)
                                                                              aaTCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGccgccg                                                                         >  2:95077/1‑150 (MQ=255)
                                                                                 ccTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCACAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGAACGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGac                                                                      >  2:887769/1‑149 (MQ=255)
                                                                                  cTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGc                                                                     <  2:261853/150‑1 (MQ=255)
                                                                                  cTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGc                                                                     <  2:98646/150‑1 (MQ=255)
                                                                                      aaTCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGAc                                                                 <  2:467682/150‑1 (MQ=255)
                                                                                        tCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCAGAGCTGACGc                                                               >  1:987541/1‑150 (MQ=255)
                                                                                          cccGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCca                                                             >  2:700268/1‑150 (MQ=255)
                                                                                          cccGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCca                                                             >  2:1670880/1‑150 (MQ=255)
                                                                                           ccGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCcac                                                            >  1:1428358/1‑150 (MQ=255)
                                                                                            cGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCcaca                                                           <  2:470096/150‑1 (MQ=255)
                                                                                            cGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCcaca                                                           <  1:1245949/150‑1 (MQ=255)
                                                                                             gCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCCGGGTTAGCTTTTGCGACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAg                                                          >  2:894995/1‑150 (MQ=255)
                                                                                              cccccTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGa                                                         <  2:392684/150‑1 (MQ=255)
                                                                                               ccccTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGaa                                                        <  1:340912/150‑1 (MQ=255)
                                                                                                 ccTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAata                                                      >  2:1330573/1‑150 (MQ=255)
                                                                                                   tGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAtt                                                    <  1:135770/150‑1 (MQ=255)
                                                                                                    ggCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttc                                                   <  2:1045376/150‑1 (MQ=255)
                                                                                                     gCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcc                                                  <  1:1670880/150‑1 (MQ=255)
                                                                                                      cTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcct                                                 >  2:454877/1‑150 (MQ=255)
                                                                                                      cTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcct                                                 >  2:41663/1‑150 (MQ=255)
                                                                                                          aaTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTccc                                             <  1:173051/150‑1 (MQ=255)
                                                                                                                   ccccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCa                                    <  1:1095488/150‑1 (MQ=255)
                                                                                                                    cccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAg                                   >  1:812927/1‑150 (MQ=255)
                                                                                                                    cccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAg                                   >  2:889510/1‑150 (MQ=255)
                                                                                                                       aaaCACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTc                                >  2:283928/1‑150 (MQ=255)
                                                                                                                         acacCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTcgc                              <  1:1032187/150‑1 (MQ=255)
                                                                                                                           acCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCgcgc                            <  2:1043027/150‑1 (MQ=255)
                                                                                                                            ccccGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCa                           <  2:1565619/150‑1 (MQ=255)
                                                                                                                                 gTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGgtg                      >  1:1102939/1‑150 (MQ=255)
                                                                                                                                      aaGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTc                 >  2:688062/1‑150 (MQ=255)
                                                                                                                                       catctgCGATCGTGGCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTcc                >  1:513627/6‑150 (MQ=255)
                                                                                                                                        ggTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCg               <  2:1413287/150‑1 (MQ=255)
                                                                                                                                         tctgCGATCGTGGCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGc              >  2:198311/4‑150 (MQ=255)
                                                                                                                                         gTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGc              <  1:700268/150‑1 (MQ=255)
                                                                                                                                          tAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCa             >  2:384267/1‑150 (MQ=255)
                                                                                                                                           aGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCAt            >  1:1303891/1‑150 (MQ=255)
                                                                                                                                            gCGATCGCCACGATCGCAGATGATCGCCACCACCCCCGCGTCAGGGTTAGCTTTTGCCCCCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATc           >  2:1642840/1‑150 (MQ=255)
                                                                                                                                                tCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCgcgc       <  1:1575510/150‑1 (MQ=255)
                                                                                                                                                   ccACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtga    >  2:323736/1‑150 (MQ=255)
                                                                                                                                                     aCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGCCGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtgatg  >  2:1332081/1‑150 (MQ=255)
                                                                                                                                                     aCGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCtgatg  <  2:87996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATG  >  NC_000913/2538589‑2538887

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