Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,952 T→A 5.7% Q316L (CAG→CTG)  cysJ ← sulfite reductase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,9520TA 5.7% 278.7 / 5.5 88Q316L (CAG→CTG) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (36/47);  new base A (2/3);  total (38/50)
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.52e-01

CACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCC  >  NC_000913/2890807‑2891101
                                                                                                                                                 |                                                                                                                                                     
cACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAg                                                                                                                                                   >  1:598323/1‑150 (MQ=255)
cACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGCGCTACAGtgg                                                                                                                                                   >  1:911744/1‑148 (MQ=255)
     tGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCtt                                                                                                                                              >  2:1341019/1‑150 (MQ=255)
       gcaacacTCGGCGGCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCCGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGGCAGTTCGAAGCGCTACAGTGGCACTTCg                                                                                                                                            <  1:1315207/149‑1 (MQ=255)
        caacaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGt                                                                                                                                           <  1:1720555/150‑1 (MQ=255)
        caacaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGt                                                                                                                                           >  2:1733326/1‑150 (MQ=255)
        caacaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCACAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTTGt                                                                                                                                           >  2:1387796/1‑150 (MQ=255)
          acaaTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTTGTGTAATTCTCAACAATTTTTGCGGTTTTTACTGTCATTTCGAAGTGCCACTTTAGCGATTCTTTc                                                                                                                                         <  2:856053/150‑1 (MQ=255)
             aTCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGCGCTACAGTGGCACTTCgaactga                                                                                                                                      <  2:911744/150‑7 (MQ=255)
                 gCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCa                                                                                                                                  >  2:813786/1‑150 (MQ=255)
                        cgcgGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGtttt                                                                                                                           >  2:2988/1‑150 (MQ=255)
                          ctgcgTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGc                                                                                                                         <  2:1168904/148‑1 (MQ=255)
                           ggcgTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGcc                                                                                                                        <  2:405340/150‑1 (MQ=255)
                            gcgTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGccc                                                                                                                       <  2:1821859/150‑1 (MQ=255)
                             cgTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCt                                                                                                                      >  1:126035/1‑150 (MQ=255)
                             cgTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCt                                                                                                                      >  2:1203953/1‑150 (MQ=255)
                              gTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTc                                                                                                                     <  1:1098919/150‑1 (MQ=255)
                                 aTGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGAc                                                                                                                  <  1:1821493/150‑1 (MQ=255)
                                 aTGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGAc                                                                                                                  <  1:143948/150‑1 (MQ=255)
                                 aTGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGAc                                                                                                                  <  2:979018/150‑1 (MQ=255)
                                  tGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACg                                                                                                                 >  1:1622097/1‑150 (MQ=255)
                                  tGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACg                                                                                                                 <  1:813786/150‑1 (MQ=255)
                                   gCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACgg                                                                                                                >  2:1358962/1‑150 (MQ=255)
                                     tGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTg                                                                                                              >  2:1670635/1‑150 (MQ=255)
                                      gTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGa                                                                                                             >  2:938973/1‑150 (MQ=255)
                                         aCTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAg                                                                                                          <  1:1341019/150‑1 (MQ=255)
                                              gCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCa                                                                                                     >  2:1457533/1‑150 (MQ=255)
                                                  tATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGTCCTCGACGGTGACAGGTTCATCGc                                                                                                 >  1:1593410/1‑150 (MQ=255)
                                                  tATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGc                                                                                                 >  1:400959/1‑150 (MQ=255)
                                                    tCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCt                                                                                               >  2:677952/1‑150 (MQ=255)
                                                     cGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCtt                                                                                              <  1:780066/150‑1 (MQ=255)
                                                     cGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCtt                                                                                              <  1:677952/150‑1 (MQ=255)
                                                       cccACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTc                                                                                            <  2:738415/150‑1 (MQ=255)
                                                         caccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAg                                                                                          <  2:898994/150‑1 (MQ=255)
                                                         caccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAg                                                                                          <  1:423598/150‑1 (MQ=255)
                                                           ccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGcc                                                                                        <  2:604317/150‑1 (MQ=255)
                                                           ccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGcc                                                                                        <  2:419027/150‑1 (MQ=255)
                                                             aGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCcac                                                                                      <  1:182431/150‑1 (MQ=255)
                                                                 gcagcaGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCAcagca                                                                                  <  1:508718/150‑1 (MQ=255)
                                                                 gcagcaGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCAcagca                                                                                  >  2:1181377/1‑150 (MQ=255)
                                                                  cagcagTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCTCCTTTCAGCCAcagcag                                                                                 <  2:613309/150‑1 (MQ=255)
                                                                  cagcagTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCAcagcag                                                                                 <  2:550687/150‑1 (MQ=255)
                                                                           ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCcac                                                                                      >  1:1166825/1‑136 (MQ=255)
                                                                           ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCcac                                                                                      <  2:1166825/136‑1 (MQ=255)
                                                                            tCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGCGCTACAGtgg                                                                                                                                                   <  1:401578/74‑3 (MQ=255)
                                                                            tCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGCGCTACAGtgg                                                                                                                                                   >  2:401578/1‑72 (MQ=255)
                                                                             cACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTGCATCGCCTTTCAGCCACAGCAGTTCGACAAGtt                                                                      >  1:1244311/1‑150 (MQ=255)
                                                                              aCTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTc                                                                     >  2:601/1‑150 (MQ=255)
                                                                                 gCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCttt                                                                  <  2:1106468/150‑1 (MQ=255)
                                                                                  cGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTc                                                                 >  2:1082843/1‑150 (MQ=255)
                                                                                  cGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTc                                                                 <  1:1181377/150‑1 (MQ=255)
                                                                                  cGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTc                                                                 <  2:1244311/150‑1 (MQ=255)
                                                                                  cGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTc                                                                 >  1:845711/1‑150 (MQ=255)
                                                                                   gggTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCa                                                                <  2:313194/150‑1 (MQ=255)
                                                                                     gTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCAcc                                                              <  1:1733326/150‑1 (MQ=255)
                                                                                       aaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAg                                                            <  2:400959/150‑1 (MQ=255)
                                                                                             ggCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCgg                                                      <  1:1671809/150‑1 (MQ=255)
                                                                                              gCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGa                                                     <  1:136780/150‑1 (MQ=255)
                                                                                                gTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATc                                                   >  1:813349/1‑150 (MQ=255)
                                                                                                gTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATc                                                   <  1:1457533/150‑1 (MQ=255)
                                                                                                   aTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGtt                                                <  1:2988/150‑1 (MQ=255)
                                                                                                      ctcAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTg                                             <  1:1408851/150‑1 (MQ=255)
                                                                                                        caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAt                                           >  1:776976/1‑150 (MQ=255)
                                                                                                        caacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGAt                                           >  2:1907627/1‑150 (MQ=255)
                                                                                                          acaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAc                                         >  1:1823073/1‑150 (MQ=255)
                                                                                                           caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATAcc                                        >  1:1753723/1‑150 (MQ=255)
                                                                                                             aTGTTGGCGGTGTTGCCGGTCAGTTCGAAGTGGCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGCCGGTGACAGGTTCCTCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGGCGGATCGTTCTGATACCAg                                      >  2:1871341/1‑150 (MQ=255)
                                                                                                              tGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGTCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGa                                     <  2:1593410/150‑1 (MQ=255)
                                                                                                              tGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGa                                     <  1:1358962/150‑1 (MQ=255)
                                                                                                               gTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGAc                                    <  1:938973/150‑1 (MQ=255)
                                                                                                                 tGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                  >  2:979821/1‑150 (MQ=255)
                                                                                                                 tGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGc                                  <  1:1670635/150‑1 (MQ=255)
                                                                                                                   gCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGccc                                <  1:852482/150‑1 (MQ=255)
                                                                                                                     ggTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGGTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAg                              <  1:1471816/150‑1 (MQ=255)
                                                                                                                      gtgtTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAgc                             >  1:1057390/1‑150 (MQ=255)
                                                                                                                         ttGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGcgcg                          <  2:598323/150‑1 (MQ=255)
                                                                                                                          tGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                         >  1:663036/1‑150 (MQ=255)
                                                                                                                          tGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGt                         <  2:845711/150‑1 (MQ=255)
                                                                                                                           gACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTc                        >  2:1226199/1‑150 (MQ=255)
                                                                                                                                  aGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                 <  2:776976/150‑1 (MQ=255)
                                                                                                                                  aGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGc                 <  1:1563117/150‑1 (MQ=255)
                                                                                                                                      cGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGt             <  1:1203953/150‑1 (MQ=255)
                                                                                                                                          gTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACg         <  1:670882/150‑1 (MQ=255)
                                                                                                                                            gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa       >  1:357579/1‑150 (MQ=255)
                                                                                                                                            gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa       <  1:601/150‑1 (MQ=255)
                                                                                                                                            gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCa       >  1:1135364/1‑150 (MQ=255)
                                                                                                                                             ccACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCAt      >  2:304329/1‑150 (MQ=255)
                                                                                                                                              cACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATg     >  1:1414060/1‑150 (MQ=255)
                                                                                                                                                 tGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGccc  >  2:796946/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCC  >  NC_000913/2890807‑2891101

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: