Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,759,788 T→C M246V (ATG→GTG)  sufS ← cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,759,7880TC99.1% 396.8 / ‑5.2 108M246V (ATG→GTG) sufScysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (61/46);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAG  >  NC_000913/1759640‑1759932
                                                                                                                                                    |                                                                                                                                                
gAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAcct                                                                                                                                                <  1:35450/151‑1 (MQ=255)
 aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAcct                                                                                                                                                <  2:317794/150‑1 (MQ=255)
 aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATGGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAcctc                                                                                                                                               <  3:212253/151‑1 (MQ=255)
 aaaCATACTCCAGCGCCGCGCAAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAcct                                                                                                                                                <  3:995033/150‑1 (MQ=255)
  aaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTcc                                                                                                                                              >  2:1312841/1‑151 (MQ=255)
  aaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTcc                                                                                                                                              <  3:1799/151‑1 (MQ=255)
   aCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAcctcct                                                                                                                                             >  2:536187/1‑151 (MQ=255)
       aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGca                                                                                                                                          >  3:322412/1‑150 (MQ=255)
        cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTCCCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGcaa                                                                                                                                         <  2:320322/150‑1 (MQ=255)
          ccAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGcaacaa                                                                                                                                      <  2:1225621/151‑1 (MQ=255)
          ccAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGcaaca                                                                                                                                       <  4:373021/150‑1 (MQ=255)
             gcgccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGc                                                                                                                                   <  4:345482/151‑1 (MQ=255)
              cgccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGct                                                                                                                                  <  1:170603/151‑1 (MQ=255)
               gccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGctt                                                                                                                                 <  1:109149/151‑1 (MQ=255)
               gccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGctt                                                                                                                                 >  3:724234/1‑151 (MQ=255)
               gccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGctt                                                                                                                                 <  2:626595/151‑1 (MQ=255)
                ccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGctt                                                                                                                                 <  3:201121/150‑1 (MQ=255)
                 cgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGcttct                                                                                                                               >  3:1079774/1‑151 (MQ=255)
                      cAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCGAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTcac                                                                                                                          >  3:607957/1‑151 (MQ=255)
                      cAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTcac                                                                                                                          >  4:2594/1‑151 (MQ=255)
                        aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACAt                                                                                                                        >  2:857251/1‑151 (MQ=255)
                         gACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACAt                                                                                                                        >  1:435527/1‑150 (MQ=255)
                          aCCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATa                                                                                                                       >  2:68180/1‑150 (MQ=255)
                           ccAATGATGCCCCCGGTATTGGGTGTACCGTCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTTCCTTCACTCAGCCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACCAGGCTTCTTTCACATaa                                                                                                                      <  3:1206026/150‑1 (MQ=255)
                           ccAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATaaa                                                                                                                     <  2:1107687/151‑1 (MQ=255)
                           ccAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATaa                                                                                                                      >  2:769388/1‑150 (MQ=255)
                             aaTGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATaa                                                                                                                      >  2:1279450/1‑148 (MQ=255)
                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGa                                                                                                                   >  2:807684/1‑150 (MQ=255)
                                   gCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGcct                                                                                                             <  3:102325/151‑2 (MQ=255)
                                   gCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGcc                                                                                                              <  2:1279415/150‑1 (MQ=255)
                                   gCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGcc                                                                                                              >  1:162890/1‑150 (MQ=255)
                                   gCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCa                                                                                                             <  2:163201/151‑1 (MQ=255)
                                    cccccGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCaa                                                                                                            <  3:395999/151‑1 (MQ=255)
                                     ccccGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAAt                                                                                                           >  3:244025/1‑151 (MQ=255)
                                      cccGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAAtt                                                                                                          <  2:45057/151‑1 (MQ=255)
                                      cccGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAAtt                                                                                                          >  3:1303175/1‑151 (MQ=255)
                                         ggTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCg                                                                                                       >  2:842340/1‑151 (MQ=255)
                                            aTTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTg                                                                                                    >  1:64959/1‑151 (MQ=255)
                                              tGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTggg                                                                                                  >  2:78815/1‑151 (MQ=255)
                                               gggTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTggg                                                                                                  <  2:1362975/150‑1 (MQ=255)
                                                   gtACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCa                                                                                              <  4:448215/150‑1 (MQ=255)
                                                     aCCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATAc                                                                                           >  2:945736/1‑151 (MQ=255)
                                                        ggCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGt                                                                                        >  1:115003/1‑151 (MQ=255)
                                                           ttCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGttt                                                                                      >  2:1326266/1‑150 (MQ=255)
                                                                 ccGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCg                                                                               >  3:1237467/1‑151 (MQ=255)
                                                                  cGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCgg                                                                              <  4:291738/151‑1 (MQ=255)
                                                                   gCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGga                                                                             <  2:280423/151‑1 (MQ=255)
                                                                    ccATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCgg                                                                              >  3:629244/1‑149 (MQ=255)
                                                                     cATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCgg                                                                              >  2:26579/1‑148 (MQ=255)
                                                                     cATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGgaga                                                                           <  2:810733/151‑1 (MQ=255)
                                                                      aTGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGgaga                                                                           >  1:192615/1‑150 (MQ=255)
                                                                      aTGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGgaga                                                                           >  2:354799/1‑150 (MQ=255)
                                                                           gCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGt                                                                     <  2:1103054/151‑1 (MQ=255)
                                                                              ttGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGACCACGTa                                                                    <  1:73058/149‑1 (MQ=255)
                                                                               tGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTaaa                                                                  <  3:1279516/150‑1 (MQ=255)
                                                                                   ccAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCg                                                             <  3:124131/151‑1 (MQ=255)
                                                                                    cAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCg                                                             >  2:910688/1‑150 (MQ=255)
                                                                                    cAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGc                                                            >  4:448902/1‑151 (MQ=255)
                                                                                        tAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAAt                                                         >  3:860895/1‑150 (MQ=255)
                                                                                        tAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATc                                                        >  3:948912/1‑151 (MQ=255)
                                                                                         aGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATcc                                                       >  2:93091/1‑151 (MQ=255)
                                                                                         aGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATcc                                                       >  4:137605/1‑151 (MQ=255)
                                                                                         aGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATc                                                        <  2:322382/150‑1 (MQ=255)
                                                                                          gTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                      <  3:1312911/151‑1 (MQ=255)
                                                                                          gTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                      <  2:1261744/151‑1 (MQ=255)
                                                                                           tGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATcc                                                       >  3:1027807/1‑149 (MQ=255)
                                                                                           tGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATcc                                                       >  2:124278/1‑149 (MQ=255)
                                                                                            gCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAgc                                                    <  2:874401/151‑1 (MQ=255)
                                                                                             ccTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAgcg                                                   >  1:294065/1‑151 (MQ=255)
                                                                                                 cACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCt                                                <  4:115002/150‑1 (MQ=255)
                                                                                                 cACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTg                                               <  3:807719/151‑1 (MQ=255)
                                                                                                      aGGCTGACGGTGACGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACAt                                          >  3:289786/1‑151 (MQ=255)
                                                                                                           gACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                      <  1:109576/150‑1 (MQ=255)
                                                                                                           gACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                      <  1:137606/150‑1 (MQ=255)
                                                                                                            aCGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAcc                                     <  2:705979/150‑1 (MQ=255)
                                                                                                                    gATCATAGAACCGCCCCCTTCCCACGGTGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGc                            >  2:1175128/1‑151 (MQ=255)
                                                                                                                    gATCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGc                            >  3:796134/1‑151 (MQ=255)
                                                                                                                     aTCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGc                            >  1:261784/1‑150 (MQ=255)
                                                                                                                      tCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCAt                          >  2:514808/1‑151 (MQ=255)
                                                                                                                      tCATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCAt                          >  2:322701/1‑151 (MQ=255)
                                                                                                                       cATAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATc                         >  2:994765/1‑151 (MQ=255)
                                                                                                                         tAGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCAc                       >  2:52820/1‑151 (MQ=255)
                                                                                                                          aGAACCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCAc                       <  1:452249/150‑1 (MQ=255)
                                                                                                                            aaCCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGc                    >  1:280308/1‑151 (MQ=255)
                                                                                                                            aaCCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGc                    >  4:258337/1‑151 (MQ=255)
                                                                                                                             aCCGCCCCCTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGcc                   >  3:559521/1‑151 (MQ=255)
                                                                                                                                   cccTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcg             <  3:910722/151‑1 (MQ=255)
                                                                                                                                   cccTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcg             >  2:1118358/1‑151 (MQ=255)
                                                                                                                                    ccTTCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcgc            >  2:1228894/1‑151 (MQ=255)
                                                                                                                                      ttCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcgc            <  2:607930/149‑1 (MQ=255)
                                                                                                                                       tCCCACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCAt         >  3:485032/1‑151 (MQ=255)
                                                                                                                                        cccACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATc        >  3:658238/1‑151 (MQ=255)
                                                                                                                                         ccACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc       <  3:847945/151‑1 (MQ=255)
                                                                                                                                         ccACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc       >  2:694989/1‑151 (MQ=255)
                                                                                                                                          cACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc       >  3:744103/1‑150 (MQ=255)
                                                                                                                                          cACGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc       >  2:384538/1‑150 (MQ=255)
                                                                                                                                           aCGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg               >  4:208305/1‑141 (MQ=255)
                                                                                                                                           aCGGCGGCACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg               <  1:208307/141‑1 (MQ=255)
                                                                                                                                              gcggcACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAg  <  3:1326336/151‑1 (MQ=255)
                                                                                                                                              gcggcACCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAg  <  3:124293/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAG  >  NC_000913/1759640‑1759932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: