Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,736 C→T 6.1% G283D (GGC→GAC) ‡ cysM ← cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,538,7360CT 6.1% 262.5 / 4.4 82G283D (GGC→GAC) ‡cysMcysteine synthase B
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base T (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 = 9.94e-01

TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGA  >  NC_000913/2538589‑2538885
                                                                                                                                                   |                                                                                                                                                     
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:1081218/150‑1 (MQ=255)
                                                          catcaGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTg                                                                                           <  2:1020768/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                                                                          >  2:135770/1‑150 (MQ=255)
                                                                           aTTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGcgccg                                                                          >  1:392684/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:41663/1‑150 (MQ=255)
                                                                                                      cTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATAttcct                                               >  2:454877/1‑150 (MQ=255)
                                                                                                          aaTGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTccc                                           <  1:173051/150‑1 (MQ=255)
                                                                                                                   ccccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCa                                  <  1:1095488/150‑1 (MQ=255)
                                                                                                                    cccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAg                                 >  2:889510/1‑150 (MQ=255)
                                                                                                                    cccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAg                                 >  1:812927/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)
                                                                                                                                                   |                                                                                                                                                     
TATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGA  >  NC_000913/2538589‑2538885

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: ‑