Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,258,251 T→A 5.4% intergenic (‑241/+34) yhaO ← / ← tdcG putative transporter/L‑serine dehydratase 3, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,258,2510TA 5.4% 319.9 / 7.7 94intergenic (‑241/+34)yhaO/tdcGputative transporter/L‑serine dehydratase 3, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (42/46);  new base A (3/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACG  >  NC_000913/3258103‑3258398
                                                                                                                                                     |                                                                                                                                                   
gAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGATGCACAt                                                                                                                                                     >  1:715011/1‑149 (MQ=255)
gAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAAGCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAGATTAATCGCACGAGAATATTCATTCTCTATTTGCGGCAGATCACAACAAAAACGTGCACAt                                                                                                                                                     >  4:900100/1‑149 (MQ=255)
 aaTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAtt                                                                                                                                                    >  2:564277/1‑149 (MQ=255)
    aaaGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtgt                                                                                                                                                 <  1:635147/149‑1 (MQ=255)
    aaaGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGACGGATAAGAATTTTCTACTATATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGCGCACATTtgt                                                                                                                                                 <  3:900100/149‑1 (MQ=255)
      aGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtgtga                                                                                                                                               >  4:411801/1‑148 (MQ=255)
          aTGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAAATGTGCAcct                                                                                                                                           <  1:417502/149‑2 (MQ=255)
              tACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAgg                                                                                                                                       >  2:787309/1‑149 (MQ=255)
               aCGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGa                                                                                                                                      <  2:558120/149‑1 (MQ=255)
                cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGAt                                                                                                                                     <  1:951981/149‑1 (MQ=255)
                     gAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAg                                                                                                                                >  4:453109/1‑149 (MQ=255)
                       ggTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCt                                                                                                                              <  1:694220/149‑1 (MQ=255)
                        gTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTg                                                                                                                             >  3:982421/1‑149 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            <  1:314855/149‑1 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            <  3:324511/149‑1 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            >  1:27054/1‑149 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            <  2:958892/149‑1 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            <  4:581894/149‑1 (MQ=255)
                         tCAAAGAAAATATTCACAATTAAGACTTTGAAATAAAGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGa                                                                                                                            <  4:39821/149‑1 (MQ=255)
                          cAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGAc                                                                                                                           >  2:6458/1‑149 (MQ=255)
                              gaAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCa                                                                                                                       >  2:987196/1‑149 (MQ=255)
                               aaaaTCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCaa                                                                                                                      >  1:526422/1‑149 (MQ=255)
                                   tCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTc                                                                                                                  <  2:907375/149‑1 (MQ=255)
                                    cTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCa                                                                                                                 <  1:136746/149‑1 (MQ=255)
                                      tCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGc                                                                                                               <  2:373530/149‑1 (MQ=255)
                                       cacaATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGcc                                                                                                              >  2:585582/1‑149 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtgt                                                                                                                                                 >  3:1074140/1‑109 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtgt                                                                                                                                                 <  4:1074140/109‑1 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                                                         >  1:595195/1‑149 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                                                         >  3:505536/1‑149 (MQ=255)
                                                gACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCAc                                                                                                     <  4:560804/149‑1 (MQ=255)
                                                 aCTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACt                                                                                                    <  1:987196/149‑1 (MQ=255)
                                                   tttGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACttt                                                                                                  >  1:162816/1‑149 (MQ=255)
                                                   tttGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACttt                                                                                                  >  1:926616/1‑149 (MQ=255)
                                                     tGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTaa                                                                                                <  4:1071526/149‑1 (MQ=255)
                                                      gAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               <  2:1056501/149‑1 (MQ=255)
                                                      gAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               <  4:491831/149‑1 (MQ=255)
                                                       aaaTAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               <  1:532930/149‑1 (MQ=255)
                                                       aaaTAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATg                                                                                              >  4:166241/1‑149 (MQ=255)
                                                       aaaTAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATg                                                                                              >  2:410711/1‑149 (MQ=255)
                                                          tAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGcc                                                                                           <  3:735009/149‑1 (MQ=255)
                                                                     gATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                                                >  3:540899/1‑149 (MQ=255)
                                                                       tAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTg                                                                              >  3:379720/1‑149 (MQ=255)
                                                                         aGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGAt                                                                            <  2:27054/149‑1 (MQ=255)
                                                                          gAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATg                                                                           >  2:144556/1‑149 (MQ=255)
                                                                           aaTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGt                                                                          <  3:15639/149‑1 (MQ=255)
                                                                              tttCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTc                                                                       <  1:163978/149‑1 (MQ=255)
                                                                                  tactaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCgg                                                                   >  4:1080432/1‑149 (MQ=255)
                                                                                  tactaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCgg                                                                   >  2:799225/1‑149 (MQ=255)
                                                                                   actaAATTAATCGCAGGAGATTATTCATTCTCTATATGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGAAAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGt                                                                  >  2:464854/1‑149 (MQ=255)
                                                                                     taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCaa                                                                                                                      <  4:526569/95‑1 (MQ=255)
                                                                                     taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCaa                                                                                                                      >  3:526569/1‑95 (MQ=255)
                                                                                      aaaTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAtt                                                               >  4:823282/1‑149 (MQ=255)
                                                                                        aTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTa                                                             >  2:965394/1‑149 (MQ=255)
                                                                                        aTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTa                                                             <  3:339986/149‑1 (MQ=255)
                                                                                          tAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATc                                                           <  1:965394/149‑1 (MQ=255)
                                                                                              cGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTc                                                       <  3:1001564/149‑1 (MQ=255)
                                                                                                 aGGAGAATATTCATTCTCTATTTGCGGCAGGTCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCata                                                    >  2:891914/1‑149 (MQ=255)
                                                                                                 aGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               <  2:172525/106‑1 (MQ=255)
                                                                                                 aGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               >  1:172525/1‑106 (MQ=255)
                                                                                                  ggAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCatat                                                   >  1:1021595/1‑149 (MQ=255)
                                                                                                   gagaATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATc                                                  <  1:144556/149‑1 (MQ=255)
                                                                                                      aaTATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCttt                                               <  2:162816/149‑1 (MQ=255)
                                                                                                          ttcattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg                                           >  2:1024665/1‑149 (MQ=255)
                                                                                                            cattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGt                                         <  1:678231/149‑1 (MQ=255)
                                                                                                             attcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGtt                                        >  4:388452/1‑149 (MQ=255)
                                                                                                              ttcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCCCACAAGGATGAAACCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGttt                                       <  2:969546/149‑1 (MQ=255)
                                                                                                              ttcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGttt                                       >  3:70513/1‑149 (MQ=255)
                                                                                                                ctctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCa                                     <  2:171713/149‑1 (MQ=255)
                                                                                                                ctctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCa                                     <  1:668257/149‑1 (MQ=255)
                                                                                                                 tctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAgg                                                                                                                                       >  2:920009/1‑50 (MQ=255)
                                                                                                                 tctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAgg                                                                                                                                       <  1:920009/50‑1 (MQ=255)
                                                                                                                 tctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt                                    >  3:424675/1‑149 (MQ=255)
                                                                                                                 tctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt                                    >  2:1098689/1‑149 (MQ=255)
                                                                                                                   tATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGtt                                                        <  4:829965/127‑1 (MQ=255)
                                                                                                                   tATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGtt                                                        >  3:829965/1‑127 (MQ=255)
                                                                                                                     tttGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACAt                                >  2:810757/1‑149 (MQ=255)
                                                                                                                      ttGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGCTGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATc                               <  2:106578/149‑1 (MQ=255)
                                                                                                                              gATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGAt                       <  3:126197/149‑1 (MQ=255)
                                                                                                                              gATCACAAAAAAAAGGTGCACAGTTGTGCACCCAAGGATGAAAGCTGACAGAAATGTCATCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCTCGGTATTGATCGTTCATATCTTTGCCGGTTTCATACATCGTCACGAt                       >  4:674038/1‑149 (MQ=255)
                                                                                                                               aTCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATc                      >  3:557774/1‑149 (MQ=255)
                                                                                                                                 cacaAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCAc                    >  2:483816/1‑149 (MQ=255)
                                                                                                                                   caAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACtt                  <  1:799225/149‑1 (MQ=255)
                                                                                                                                     aaaaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACttt                 >  2:877746/1‑149 (MQ=255)
                                                                                                                                     aaaaaaaaGGGGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACttt                 <  1:825062/149‑1 (MQ=255)
                                                                                                                                       aaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCgg                                                                   >  1:981315/1‑97 (MQ=255)
                                                                                                                                       aaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCgg                                                                   <  2:981315/97‑1 (MQ=255)
                                                                                                                                        aaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATc              <  1:483816/149‑1 (MQ=255)
                                                                                                                                          ttgggTGCACAAATGTGCACCTAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGa            >  1:984591/4‑149 (MQ=255)
                                                                                                                                          ttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtc                            >  4:263913/4‑133 (MQ=255)
                                                                                                                                          ttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtc                            <  3:263913/130‑1 (MQ=255)
                                                                                                                                           tgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAg           >  3:1034806/3‑149 (MQ=255)
                                                                                                                                            aGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGt          >  1:17870/1‑149 (MQ=255)
                                                                                                                                            aGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGt          <  3:520217/149‑1 (MQ=255)
                                                                                                                                                gCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGaaa      <  4:379720/149‑1 (MQ=255)
                                                                                                                                                   caTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAacac   <  3:246877/149‑1 (MQ=255)
                                                                                                                                                    aTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACAcg  >  2:906566/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACG  >  NC_000913/3258103‑3258398

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