Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,092,779 Δ1 bp 11.4% coding (902/951 nt) gshB → glutathione synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,092,7790T.11.4% 638.9 / 72.4 184coding (902/951 nt)gshBglutathione synthetase
Reads supporting (aligned to +/- strand):  ref base T (79/84);  new base . (10/11);  total (89/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 189 total aligned reads displayed.

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926
                                                                                                                                                     |                                                                                                                                                   
tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAATGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  3:952360/1‑149 (MQ=255)
tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  4:491273/1‑149 (MQ=255)
tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGAt                                                                                                                                                     >  1:651259/1‑149 (MQ=255)
        cGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaa                                                                                                                                             <  1:1289015/149‑1 (MQ=255)
         gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                            <  3:147104/149‑2 (MQ=255)
           ccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGt                                                                                                                                          <  3:491273/149‑1 (MQ=255)
               cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      <  2:651259/149‑1 (MQ=255)
                       gCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCa                                                                                                                              >  4:1173446/1‑149 (MQ=255)
                                gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                                     >  2:792157/1‑149 (MQ=255)
                                gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                                     <  1:792157/149‑1 (MQ=255)
                                   agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                                  >  1:434810/1‑149 (MQ=255)
                                    gaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                            <  4:490942/122‑2 (MQ=255)
                                    gaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGA‑CACCGGaaa                                                                                                                                            >  3:490942/1‑121 (MQ=255)
                                     aaaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGaa                                                                                                                         >  1:1411649/1‑140 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCGGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               >  3:919535/1‑149 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               <  2:1192900/149‑1 (MQ=255)
                                           ggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                          <  2:1453826/149‑1 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                                  >  2:764205/1‑140 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                                  <  1:764205/140‑1 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagc                                                                                                         >  2:1290943/1‑149 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                             <  4:577180/145‑1 (MQ=255)
                                                          gTCTGGATATCATCGGCGACCGTCTGACAGAGATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  4:899339/1‑149 (MQ=255)
                                                           tCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                           >  4:106227/1‑148 (MQ=255)
                                                           tCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAg                                                                                          <  2:269257/149‑1 (MQ=255)
                                                               gATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCgag                                                                                      >  2:382869/1‑149 (MQ=255)
                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  2:1177636/149‑1 (MQ=255)
                                                                   tcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                           <  1:1257316/124‑1 (MQ=255)
                                                                   tcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                           >  2:1257316/1‑124 (MQ=255)
                                                                    catcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGa                                                                                 >  2:310657/1‑149 (MQ=255)
                                                                        ggCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                             >  1:1182546/1‑149 (MQ=255)
                                                                        ggCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtctc                                                                             >  1:524531/1‑149 (MQ=255)
                                                                         gCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGaa                                                                                                                                             >  3:796385/1‑84 (MQ=255)
                                                                         gCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGaa                                                                                                                                             <  4:796385/84‑1 (MQ=255)
                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      <  2:1367604/89‑1 (MQ=255)
                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      <  3:987073/89‑1 (MQ=255)
                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGAt                                                                                                                                  >  2:700829/1‑93 (MQ=255)
                                                                            aCCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                        >  3:772200/1‑118 (MQ=255)
                                                                            aCCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                        <  4:772200/118‑1 (MQ=255)
                                                                             ccGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        <  4:979370/149‑1 (MQ=255)
                                                                               gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                      <  2:524531/149‑1 (MQ=255)
                                                                                    actgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  2:1084045/130‑1 (MQ=255)
                                                                                    actgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    >  1:1084045/1‑130 (MQ=255)
                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAc                                                                                                   >  1:147221/1‑114 (MQ=255)
                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  3:591761/149‑1 (MQ=255)
                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  3:661606/149‑1 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGaa                                                                                                                                             <  4:604362/71‑1 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAAccc                                                                                                 >  3:1198163/1‑115 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    >  1:100019/1‑128 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  4:515506/1‑149 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  2:155480/1‑149 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  2:933577/1‑149 (MQ=255)
                                                                                       gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              <  1:310657/149‑1 (MQ=255)
                                                                                       gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              >  1:884913/1‑149 (MQ=255)
                                                                                         aaTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                                     <  1:185126/92‑1 (MQ=255)
                                                                                         aaTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGa                                                            <  3:899339/149‑1 (MQ=255)
                                                                                          aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGt                                                                                                                                          <  3:1401744/70‑1 (MQ=255)
                                                                                                 tCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        >  1:1299381/1‑129 (MQ=255)
                                                                                                  cACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                >  3:1108001/1‑136 (MQ=255)
                                                                                                   aCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCctt                                                  >  1:284815/1‑149 (MQ=255)
                                                                                                     cAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAAt                                                                                                                                      >  4:203778/1‑63 (MQ=255)
                                                                                                        cccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  4:574781/1‑131 (MQ=255)
                                                                                                        cccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               <  3:574781/131‑1 (MQ=255)
                                                                                                         ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                            <  1:382869/149‑1 (MQ=255)
                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                         <  1:394126/119‑1 (MQ=255)
                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                         >  2:394126/1‑119 (MQ=255)
                                                                                                            aCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACg                                         >  3:581948/1‑149 (MQ=255)
                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAg                                                                                   >  1:763153/1‑105 (MQ=255)
                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                       <  1:933577/149‑1 (MQ=255)
                                                                                                                    tCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTgg                                 <  2:884913/149‑1 (MQ=255)
                                                                                                                      gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                                <  1:343560/100‑1 (MQ=255)
                                                                                                                      gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                                >  2:343560/1‑100 (MQ=255)
                                                                                                                      gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTgg                                 <  4:591689/147‑1 (MQ=255)
                                                                                                                      gTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgc                               >  3:1313997/1‑149 (MQ=255)
                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  3:753154/1‑105 (MQ=255)
                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:753154/105‑1 (MQ=255)
                                                                                                                           aTTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGt                          <  2:434810/149‑1 (MQ=255)
                                                                                                                                aGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtt                     >  2:1015856/1‑149 (MQ=255)
                                                                                                                                 gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGc                                                                                         >  1:690312/1‑80 (MQ=255)
                                                                                                                                 gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        >  3:1399297/1‑97 (MQ=255)
                                                                                                                                 gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        <  4:1399297/97‑1 (MQ=255)
                                                                                                                                 gCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttt                    >  4:166362/1‑149 (MQ=255)
                                                                                                                                     agTTTCCGGTG‑TCAATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCAAGAAGGATCTCGTTGAGACTCTGa                                                                >  4:944925/1‑101 (MQ=255)
                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCAcc                                                                                              <  4:378136/69‑1 (MQ=255)
                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCAcc                                                                                              >  3:378136/1‑69 (MQ=255)
                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              <  4:274397/101‑1 (MQ=255)
                                                                                                                                       tttCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              >  3:274397/1‑101 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAGCCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa              <  3:1051251/148‑1 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                                >  3:196538/1‑82 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        >  4:422083/1‑90 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                        <  3:422083/90‑1 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCctt                                                  >  4:1218520/1‑112 (MQ=255)
                                                                                                                                        ttCCGGTGATCGA‑CACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt                 >  3:500102/1‑145 (MQ=255)
                                                                                                                                          ccGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCa           <  2:292124/149‑1 (MQ=255)
                                                                                                                                           cGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAg          >  1:962265/1‑149 (MQ=255)
                                                                                                                                           cGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAg          <  1:155480/149‑1 (MQ=255)
                                                                                                                                             gtg‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  4:1441953/1‑83 (MQ=255)
                                                                                                                                              tg‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               <  2:763505/93‑1 (MQ=255)
                                                                                                                                               g‑tCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGaaa      <  1:1015856/149‑1 (MQ=255)
                                                                                                                                                    aTCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  <  4:1313997/149‑1 (MQ=255)
                                                                                                                                                    aTCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  >  3:24416/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTG‑TCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑