Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,385,837 T→A 5.5% intergenic (+36/+20) yhcN → / ← yhcO cadmium and peroxide resistance protein, stress‑induced/putative barnase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,385,8370TA 5.5% 430.6 / 5.6 158intergenic (+36/+20)yhcN/yhcOcadmium and peroxide resistance protein, stress‑induced/putative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base T (79/68);  new base A (3/7);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGC  >  NC_000913/3385690‑3385986
                                                                                                                                                   |                                                                                                                                                      
  tGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCgg                                                                                                                                                   >  1:1288800/1‑150 (MQ=255)
   gCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                                                                                                                                                  >  1:1417374/1‑150 (MQ=255)
   gCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                                                                                                                                                  >  2:1590198/1‑150 (MQ=255)
   gCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                                                                                                                                                  <  1:1568990/150‑1 (MQ=255)
      tGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGc                                                                                                                                               <  1:1249375/150‑1 (MQ=255)
      tGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGAACCGGGGGc                                                                                                                                               <  1:126466/150‑1 (MQ=255)
      ggAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCGCGTCGCGGGGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCCTCGTCTTGTCCGCCGATATTGCCCCCGAACCGGGGGc                                                                                                                                               <  1:1216155/149‑1 (MQ=255)
       gAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCt                                                                                                                                              >  1:1676097/1‑150 (MQ=255)
       gAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCt                                                                                                                                              >  1:1694644/1‑150 (MQ=255)
          aTGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATACACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttt                                                                                                                                           >  1:339265/1‑150 (MQ=255)
          aTGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGAACCGGGGGcaata                                                                                                                                           <  2:956468/150‑5 (MQ=255)
           tGCTGAACAAAAAAGCGGAAGAGAAAGGCTCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGAACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGAACCGGGGGcaatat                                                                                                                                          <  1:1522933/150‑6 (MQ=255)
           tGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGAACCGGGGGcaatat                                                                                                                                          <  1:1340152/150‑6 (MQ=255)
            gCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATTAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttttt                                                                                                                                         <  2:697355/150‑1 (MQ=255)
             cTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCttttttt                                                                                                                                        >  1:239952/1‑150 (MQ=255)
              tGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGCCACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttttttt                                                                                                                                       >  2:956606/1‑150 (MQ=255)
                aaCAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTgc                                                                                                                                     <  2:935633/150‑1 (MQ=255)
                 aCAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTgcg                                                                                                                                    <  2:457405/150‑1 (MQ=255)
                  cAAAAAAGCGGAAGAGAAAGTCGCACCGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACCACTGGCACGCTACGGCTGATCAGGACAAATTAACCCTCATCGGCTTGGCCGACGATATTGCCGGCGGAGCGGGGGCTTTTTTTTgcgt                                                                                                                                   >  1:106995/1‑149 (MQ=255)
                  cAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTgcgc                                                                                                                                   >  1:64418/1‑150 (MQ=255)
                    aaaaaGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTa                                                                                                                                 >  1:1189212/1‑150 (MQ=255)
                        aGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGa                                                                                                                             <  1:1000574/150‑1 (MQ=255)
                          cggtacagaAAGGCGCTACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACATATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACg                                                                                                                           <  1:188768/144‑1 (MQ=255)
                              agagAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCgggg                                                                                                                                                 <  2:1527420/124‑1 (MQ=255)
                              agagAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCgggg                                                                                                                                                 >  1:1527420/1‑124 (MQ=255)
                               gagaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAAGTGTTCAAATAAACCCTCATCGTCGTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACTAAcct                                                                                                                      >  1:220868/1‑148 (MQ=255)
                                agaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACAtt                                                                                                                     >  1:1634382/1‑150 (MQ=255)
                                 gaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTa                                                                                                                    >  1:777009/1‑150 (MQ=255)
                                 gaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTa                                                                                                                    <  1:1443436/150‑1 (MQ=255)
                                 gaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTa                                                                                                                    >  2:777488/1‑150 (MQ=255)
                                    aGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTaaaa                                                                                                                 <  1:954631/150‑1 (MQ=255)
                                    aGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACAATCATCGTCTTGTCCGACGATATTGCCCCCCGTTCGGGCGCTAAATTTTCCGCTAATGCCGAACATTaaaa                                                                                                                 >  2:1243361/1‑150 (MQ=255)
                                       cgcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGc                                                                                                              <  1:437523/150‑1 (MQ=255)
                                        gcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCa                                                                                                             <  2:1479705/150‑1 (MQ=255)
                                        gcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCa                                                                                                             <  2:777009/150‑1 (MQ=255)
                                        gcAACGGCCTACAAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCa                                                                                                             <  2:347764/150‑1 (MQ=255)
                                                 tACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCtt                                                                                                    >  2:1158690/1‑150 (MQ=255)
                                                  aCCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTCGCTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                                   >  2:111271/1‑150 (MQ=255)
                                                  aCCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                                   >  1:1005106/1‑150 (MQ=255)
                                                   ccAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                  <  2:1339832/150‑1 (MQ=255)
                                                    cAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGTGCTAATGACGAACATTAAAACGCAAATGCCCTTCCa                                                                                                 <  1:457284/150‑1 (MQ=255)
                                                    cAGATAACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCa                                                                                                 <  1:1424858/150‑1 (MQ=255)
                                                     aGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAg                                                                                                <  1:1460417/150‑1 (MQ=255)
                                                         tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                                            >  1:497915/1‑150 (MQ=255)
                                                         tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                                            <  2:1636146/150‑1 (MQ=255)
                                                         tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                                            <  2:1417374/150‑1 (MQ=255)
                                                          aCTGAAGCTCCTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCtt                                                                                           >  2:547758/1‑150 (MQ=255)
                                                            tGAAGCTCGTACCGGTGACAACTGCCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTcc                                                                                         <  1:334482/150‑1 (MQ=255)
                                                             gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTccg                                                                                        >  1:780091/1‑149 (MQ=255)
                                                             gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGACCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCct                                                                                        >  1:959322/1‑150 (MQ=255)
                                                              aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCct                                                                                        <  2:505480/150‑1 (MQ=255)
                                                              aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctc                                                                                       >  2:759960/1‑150 (MQ=255)
                                                               aGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGATCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctct                                                                                      <  1:613262/150‑1 (MQ=255)
                                                               aGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctct                                                                                      <  2:1562220/150‑1 (MQ=255)
                                                                gCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTg                                                                                     <  2:984697/150‑1 (MQ=255)
                                                                gCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTg                                                                                     >  1:1373787/1‑150 (MQ=255)
                                                                   cGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCt                                                                                  <  2:1615496/150‑1 (MQ=255)
                                                                   cGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGGTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCt                                                                                  <  1:1677391/150‑1 (MQ=255)
                                                                     tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCtca                                                                                <  1:1590198/150‑1 (MQ=255)
                                                                            gACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACa                                                                         <  1:1660427/150‑1 (MQ=255)
                                                                              cacCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGc                                                                       <  2:1446739/150‑1 (MQ=255)
                                                                                ccTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCaa                                                                     >  1:1032013/1‑150 (MQ=255)
                                                                                 cTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAat                                                                    <  1:582526/150‑1 (MQ=255)
                                                                                 cTGACACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAat                                                                    >  2:1492275/1‑150 (MQ=255)
                                                                                  tGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAata                                                                   >  2:1126528/1‑150 (MQ=255)
                                                                                     cACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTa                                                                <  2:375048/150‑1 (MQ=255)
                                                                                          tACGGCTGAACTGTACTAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTGTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGATCTTCCTCTGCCTCATCAAACAGCAATATTAActcc                                                           >  2:1655342/1‑147 (MQ=255)
                                                                                             ggCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCa                                                        <  2:664011/150‑1 (MQ=255)
                                                                                                   aCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCACCCGGTTCGGGGGCTGATTTTTGCGCTAATGCCGAACATTAAAACGCAACTGCACTTCCAGCTCTTCATCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACgc                                                  >  2:1319636/1‑149 (MQ=255)
                                                                                                         cAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGAcgc                                            >  1:354150/1‑150 (MQ=255)
                                                                                                         cAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGAcgc                                            >  1:933339/1‑150 (MQ=255)
                                                                                                             tAAACCCTCATCGTCTTGTCCGACGATAATGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAACACGAAAATGCCCTTCCAGCTCTTGCTCTGCCTCATCATAAAGCAATATTAAAGAGCCAAAACGGCGCcgaattt                                        >  1:106046/1‑145 (MQ=255)
                                                                                                              aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGGCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                       <  1:899433/150‑1 (MQ=255)
                                                                                                              aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                       >  2:445391/1‑150 (MQ=255)
                                                                                                                aCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtct                                     <  2:581334/150‑1 (MQ=255)
                                                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctc                                    >  2:871778/1‑150 (MQ=255)
                                                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATAAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctc                                    <  2:1096701/150‑1 (MQ=255)
                                                                                                                   cTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctctc                                  >  1:189363/1‑150 (MQ=255)
                                                                                                                      aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCa                               <  1:577301/150‑1 (MQ=255)
                                                                                                                         gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGAt                            <  2:1032013/150‑1 (MQ=255)
                                                                                                                         gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTGTTCTCTCCCAGATg                           >  2:597151/1‑150 (MQ=255)
                                                                                                                            ttGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGaa                         <  2:663465/150‑1 (MQ=255)
                                                                                                                              gTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACa                       >  1:1314442/1‑150 (MQ=255)
                                                                                                                               tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACaa                      >  2:838896/1‑150 (MQ=255)
                                                                                                                               tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACaa                      >  2:491928/1‑150 (MQ=255)
                                                                                                                               tCCGACGATATTGCCCCCGGTTAGGGGGCTATTATTTGCGCTAATGACTAACATTAAAACTCAAATGCACTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCAAAAACGGCGACTCTGATTCTCTCCCAGATGAACaa                      >  1:102556/1‑150 (MQ=255)
                                                                                                                                  gacgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGATAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGGTTTCTCTCCCAGCTGAACAAAtt                   >  2:1478645/1‑150 (MQ=255)
                                                                                                                                   acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTc                  <  1:995802/150‑1 (MQ=255)
                                                                                                                                   acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTc                  <  2:74149/150‑1 (MQ=255)
                                                                                                                                    aaaaaaagCCCCCGAACCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCa                 >  2:831375/8‑150 (MQ=255)
                                                                                                                                    aaaaaaagCCCCCGAACCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCa                 >  2:233790/8‑150 (MQ=255)
                                                                                                                                      atatTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAAt               <  2:423485/150‑1 (MQ=255)
                                                                                                                                        atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCATACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAAtct             <  2:64418/150‑1 (MQ=255)
                                                                                                                                        atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAAtct             >  1:818385/1‑150 (MQ=255)
                                                                                                                                          tGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCa           >  2:1386351/1‑150 (MQ=255)
                                                                                                                                           gCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTGCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAACCGGCGACGCGATTTCTCTCCCAGATGAACAAATTCAATCTCaa          >  1:1340198/1‑150 (MQ=255)
                                                                                                                                              cccGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTATCCCAGATGAACAAATACAATATCAAGTg       >  2:563920/1‑150 (MQ=255)
                                                                                                                                               ccGGTTCTGGGGCTTTTTTTTGCGCTAATGACCAACATTAAACCGCAAATGCCCGACCCGCACTTCCTCTGCCTCATCAAACATAATTAATAACGCGCCAAAACGGCGACGCGGTTTCTCTCCCAGATGAACAACTTCAATCTCAAGTgg      <  2:910937/150‑1 (MQ=255)
                                                                                                                                                 ggTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTa    >  1:945976/1‑150 (MQ=255)
                                                                                                                                                   ttCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGc  <  1:1624731/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGC  >  NC_000913/3385690‑3385986

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