Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,092,779 Δ1 bp 6.5% 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. 6.5% 324.4 / 14.8 92coding (902/951 nt)gshBglutathione synthetase
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base . (4/2);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926
                                                                                                                                                     |                                                                                                                                                   
tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGAt                                                                                                                                                     <  3:228890/149‑1 (MQ=255)
tGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGAt                                                                                                                                                     >  2:432428/1‑149 (MQ=255)
  gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCa                                                                                                                                                   >  3:449851/1‑149 (MQ=255)
  gAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTCGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCa                                                                                                                                                   >  2:407410/1‑149 (MQ=255)
        cGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGaa                                                                                                                                             <  2:420124/149‑1 (MQ=255)
         gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGACACCGGaaa                                                                                                                                            <  4:10628/149‑2 (MQ=255)
           ccGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGt                                                                                                                                          <  1:63797/149‑1 (MQ=255)
                  aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGa                                                                                                                                   >  3:338620/1‑149 (MQ=255)
                   tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGAt                                                                                                                                  <  3:66357/149‑1 (MQ=255)
                       gCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCa                                                                                                                              <  3:14417/149‑1 (MQ=255)
                         cGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATc                                                                                                                            >  3:139191/1‑149 (MQ=255)
                            cGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGaa                                                                                                                         <  4:5492/149‑1 (MQ=255)
                            cGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGaa                                                                                                                         <  4:338620/149‑1 (MQ=255)
                               tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                      >  2:406891/1‑149 (MQ=255)
                               tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                      >  2:230650/1‑149 (MQ=255)
                                gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                                     <  2:156996/149‑1 (MQ=255)
                                    gaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGttt                                                                                                                 >  1:386277/1‑149 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               >  1:180446/1‑149 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               >  1:206515/1‑149 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               <  1:230650/149‑1 (MQ=255)
                                      aaaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAc                                                                                                               >  1:441256/1‑149 (MQ=255)
                                       aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAca                                                                                                              <  3:202779/149‑1 (MQ=255)
                                          gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACGTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAg                                                                                                                        >  3:331665/1‑136 (MQ=255)
                                          gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAg                                                                                                                        >  1:342983/1‑136 (MQ=255)
                                          gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAg                                                                                                                        <  2:342983/136‑1 (MQ=255)
                                          gggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagca                                                                                                           <  4:449851/149‑1 (MQ=255)
                                          gggCTGATTTATGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAg                                                                                                                        <  4:331665/136‑1 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagc                                                                                                         >  1:443984/1‑149 (MQ=255)
                                            gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagc                                                                                                         >  2:330613/1‑149 (MQ=255)
                                              tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                       <  1:406891/149‑1 (MQ=255)
                                              tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                       >  3:349139/1‑149 (MQ=255)
                                                 tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTaa                                                                                                    >  4:204834/1‑149 (MQ=255)
                                                              ggATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCga                                                                                       >  1:396664/1‑149 (MQ=255)
                                                                atatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAAt                                                                                                                                            <  3:346129/94‑1 (MQ=255)
                                                                atatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAAt                                                                                                                                            >  4:346129/1‑94 (MQ=255)
                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  1:330613/149‑1 (MQ=255)
                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    >  3:462900/1‑149 (MQ=255)
                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  3:204834/149‑1 (MQ=255)
                                                                 tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  2:396664/149‑1 (MQ=255)
                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  4:342686/1‑149 (MQ=255)
                                                                           gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:139509/149‑1 (MQ=255)
                                                                            aCCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTg                                                                         >  2:141221/1‑149 (MQ=255)
                                                                               gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                      <  2:331533/149‑1 (MQ=255)
                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  3:342686/149‑1 (MQ=255)
                                                                                     ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  4:437909/149‑1 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  3:359344/1‑149 (MQ=255)
                                                                                      tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAg                                                               >  4:152881/1‑149 (MQ=255)
                                                                                       gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              <  2:206515/149‑1 (MQ=255)
                                                                                       gaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGt                                                              >  3:165602/1‑149 (MQ=255)
                                                                                          aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGAc                                                           >  2:250182/1‑149 (MQ=255)
                                                                                          aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGAc                                                           >  3:322870/1‑149 (MQ=255)
                                                                                           ttAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACa                                                          <  1:407410/149‑1 (MQ=255)
                                                                                                  cACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                        >  3:343344/1‑96 (MQ=255)
                                                                                                  cACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                        <  4:343344/96‑1 (MQ=255)
                                                                                                     cAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttct                                                >  4:89239/1‑149 (MQ=255)
                                                                                                       gCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  3:277802/1‑121 (MQ=255)
                                                                                                       gCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:277802/121‑1 (MQ=255)
                                                                                                        cccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTc                                             <  3:152881/149‑1 (MQ=255)
                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCga                                                                                       >  3:288819/1‑105 (MQ=255)
                                                                                                          cAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCga                                                                                       <  4:288819/105‑1 (MQ=255)
                                                                                                           aaCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCAc                                          >  1:51412/1‑149 (MQ=255)
                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  3:210749/1‑114 (MQ=255)
                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:210749/114‑1 (MQ=255)
                                                                                                              cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                       >  2:196292/1‑149 (MQ=255)
                                                                                                               tGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCAt                                      >  4:42069/1‑149 (MQ=255)
                                                                                                                   ttCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    >  2:5925/1‑99 (MQ=255)
                                                                                                                   ttCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                                    <  1:5925/99‑1 (MQ=255)
                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  3:14625/1‑105 (MQ=255)
                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:14625/105‑1 (MQ=255)
                                                                                                                       tGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGgcg                              >  3:343495/1‑149 (MQ=255)
                                                                                                                          gaTTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTg                           >  1:208204/1‑149 (MQ=255)
                                                                                                                            ttGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGtt                         >  2:159180/1‑149 (MQ=255)
                                                                                                                            ttGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGtt                         <  3:89239/149‑1 (MQ=255)
                                                                                                                               aaGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                      <  1:159180/149‑1 (MQ=255)
                                                                                                                                aGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                >  3:19125/1‑106 (MQ=255)
                                                                                                                                aGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                                <  4:19125/106‑1 (MQ=255)
                                                                                                                                 aataacattCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttt                    >  2:378644/8‑149 (MQ=255)
                                                                                                                                   agagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt                  >  1:440482/1‑149 (MQ=255)
                                                                                                                                   agagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt                  <  4:462900/149‑1 (MQ=255)
                                                                                                                                   agagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt                  <  4:349139/149‑1 (MQ=255)
                                                                                                                                    gagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGTAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt                 <  2:386277/149‑1 (MQ=255)
                                                                                                                                    gagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt                 <  1:432428/149‑1 (MQ=255)
                                                                                                                                       tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCt                            <  2:191396/135‑1 (MQ=255)
                                                                                                                                       tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCt                            >  1:191396/1‑135 (MQ=255)
                                                                                                                                       tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa              <  2:441256/149‑1 (MQ=255)
                                                                                                                                       tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa              <  2:180446/149‑1 (MQ=255)
                                                                                                                                        ttCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          >  3:1250/1‑88 (MQ=255)
                                                                                                                                        ttCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                          <  4:1250/88‑1 (MQ=255)
                                                                                                                                        ttCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc             >  2:109093/1‑149 (MQ=255)
                                                                                                                                        ttCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAAc             >  1:228655/1‑149 (MQ=255)
                                                                                                                                             gtgtCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                      >  2:322566/1‑87 (MQ=255)
                                                                                                                                             gtgtCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                      <  1:322566/87‑1 (MQ=255)
                                                                                                                                                    aTCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGCACCAGGAAACAGa  >  4:153800/1‑149 (MQ=255)
                                                                                                                                                    aTCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGa  >  4:34385/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGAAACAGA  >  NC_000913/3092631‑3092926

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