Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 870,890 A→T 7.1% T494S (ACC→TCC) ‡ gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913870,8900AT 7.1% 381.6 / 16.1 113T494S (ACC→TCC) ‡gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (53/52);  new base T (3/5);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 113 total aligned reads displayed.

ATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGT  >  NC_000913/870744‑871036
                                                                                                                                                  |                                                                                                                                                  
aTGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGAc                                                                                                                                                   >  1:905215/1‑148 (MQ=255)
aTGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGAc                                                                                                                                                   >  1:514494/1‑148 (MQ=255)
        cTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                           <  2:514494/148‑1 (MQ=255)
        cTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTg                                                                                                                                          <  4:4142/149‑1 (MQ=255)
        cTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGttttta                                                                                                                                          <  4:588491/149‑2 (MQ=255)
        cTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGttttta                                                                                                                                          <  3:728458/149‑2 (MQ=255)
        cTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCACAAGATATCATCTGGCAAGAATCGCCGTGGATACCGCTGGTGGTAGAAAAATTGGTTTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                           <  1:131445/148‑1 (MQ=255)
           gCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                        <  1:180538/148‑1 (MQ=255)
           gCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                        <  2:657070/148‑1 (MQ=255)
           gCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                        <  2:905215/148‑1 (MQ=255)
            cTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGtttttactg                                                                                                                                       <  1:101621/148‑5 (MQ=255)
            cTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGtttttactg                                                                                                                                       <  1:480174/148‑5 (MQ=255)
             tCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCa                                                                                                                                     >  3:772999/1‑149 (MQ=255)
             tCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCa                                                                                                                                     >  4:309014/1‑149 (MQ=255)
              cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCa                                                                                                                                     >  2:796782/1‑148 (MQ=255)
              cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCa                                                                                                                                     >  2:258775/1‑148 (MQ=255)
                ggCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                   >  2:558076/1‑148 (MQ=255)
                ggCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGc                                                                                                                                  <  3:256738/149‑1 (MQ=255)
                    cTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAg                                                                                                                               >  2:258445/1‑148 (MQ=255)
                     tGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGTTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGac                                                                                                                             >  4:502500/1‑149 (MQ=255)
                      gAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGac                                                                                                                             >  1:533982/1‑148 (MQ=255)
                            tAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtttt                                                                                                                                            >  4:170777/1‑127 (MQ=255)
                            tAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCt                                                                                                                      >  3:731021/1‑149 (MQ=255)
                             aaTGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCt                                                                                                                      >  1:312517/1‑148 (MQ=255)
                              aTGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTc                                                                                                                    <  3:749715/149‑1 (MQ=255)
                              aTGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTc                                                                                                                    >  4:329636/1‑149 (MQ=255)
                                 aTCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGc                                                                                                                 >  3:28858/1‑149 (MQ=255)
                                    cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                               <  1:602756/148‑1 (MQ=255)
                                    cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                               <  2:435066/148‑1 (MQ=255)
                                    cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                               <  1:199462/148‑1 (MQ=255)
                                    cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCAGCTTCAGCttt                                                                                                              <  3:470668/149‑1 (MQ=255)
                                       cggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGa                                                                                                            >  2:402669/1‑148 (MQ=255)
                                          aaaaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAc                                                                                                        >  4:590870/1‑149 (MQ=255)
                                           aaaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAc                                                                                                        >  2:704237/1‑148 (MQ=255)
                                                 ccGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAtt                                                                                                 >  4:821250/1‑149 (MQ=255)
                                                  cGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAtt                                                                                                 >  1:878995/1‑148 (MQ=255)
                                                    cTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGc                                                                                              <  4:219519/149‑1 (MQ=255)
                                                    cTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGc                                                                                              <  4:625677/149‑1 (MQ=255)
                                                         ataAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                   <  2:777892/107‑1 (MQ=255)
                                                         ataAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                   >  1:777892/1‑107 (MQ=255)
                                                           aaGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGc                                                                                       >  3:679589/1‑149 (MQ=255)
                                                                   gcAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAg                                                                                <  1:402669/148‑1 (MQ=255)
                                                                     aGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGa                                                                             <  3:502500/149‑1 (MQ=255)
                                                                       gATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAg                                                                            >  1:793416/1‑148 (MQ=255)
                                                                        atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                           >  2:481511/1‑148 (MQ=255)
                                                                        atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                           <  1:704237/148‑1 (MQ=255)
                                                                         tatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                         >  3:346386/1‑149 (MQ=255)
                                                                          atcatcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                         >  1:680507/1‑148 (MQ=255)
                                                                          atcatcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                        <  3:237055/149‑1 (MQ=255)
                                                                             atcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                     >  4:585872/1‑149 (MQ=255)
                                                                             atcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                     >  3:803523/1‑149 (MQ=255)
                                                                              tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                     <  1:244834/148‑1 (MQ=255)
                                                                              tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                     >  1:658817/1‑148 (MQ=255)
                                                                               cTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGc                                                                    >  2:757455/1‑148 (MQ=255)
                                                                                 ggCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTa                                                                 >  4:30480/1‑149 (MQ=255)
                                                                                  gCAAGAATCGCCGTAGATCCCGCTGTTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGATTTTTGGATCATGCCAGACACCGGCTTAAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTa                                                                 <  2:374782/148‑1 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAtt                                                              >  2:65385/1‑148 (MQ=255)
                                                                                        aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                           <  1:602931/148‑1 (MQ=255)
                                                                                        aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                           <  1:464968/148‑1 (MQ=255)
                                                                                        aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGt                                                          <  4:346386/149‑1 (MQ=255)
                                                                                             cGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATc                                                      <  2:680507/148‑1 (MQ=255)
                                                                                                  aTCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACg                                                 >  1:45909/1‑148 (MQ=255)
                                                                                                     ccGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCtt                                              >  1:754563/1‑148 (MQ=255)
                                                                                                     ccGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTa                                             <  3:595093/149‑1 (MQ=255)
                                                                                                      cGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTa                                             >  1:65419/1‑148 (MQ=255)
                                                                                                        cTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTg                                          <  3:585872/149‑1 (MQ=255)
                                                                                                             ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCaa                                                                                     <  4:843089/101‑1 (MQ=255)
                                                                                                             ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGtt                                     <  4:679589/149‑1 (MQ=255)
                                                                                                              gtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                           >  2:739525/1‑110 (MQ=255)
                                                                                                              gtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTg                                    >  3:496247/1‑149 (MQ=255)
                                                                                                               tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                   <  3:329636/149‑1 (MQ=255)
                                                                                                               tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                   <  4:803523/149‑1 (MQ=255)
                                                                                                               tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                   <  4:772999/149‑1 (MQ=255)
                                                                                                                aGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                   >  2:514940/1‑148 (MQ=255)
                                                                                                                aGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAt                                  <  3:590870/149‑1 (MQ=255)
                                                                                                                aGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAt                                  >  3:393229/1‑149 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTggg                                        <  3:531828/142‑1 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTggg                                        >  4:531828/1‑142 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAtt                                 >  4:601988/1‑149 (MQ=255)
                                                                                                                    aaaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTcc                               >  1:512935/1‑148 (MQ=255)
                                                                                                                     aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCg                              >  1:149188/1‑148 (MQ=255)
                                                                                                                      aCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGa                             <  2:65419/148‑1 (MQ=255)
                                                                                                                      aCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGAc                            >  3:431881/1‑149 (MQ=255)
                                                                                                                         ggTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGc                          >  2:584137/1‑148 (MQ=255)
                                                                                                                           tgtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTg                        <  2:878995/148‑1 (MQ=255)
                                                                                                                                gCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTAt                   <  2:352702/148‑1 (MQ=255)
                                                                                                                                 cacacaGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATc                  >  2:777037/1‑148 (MQ=255)
                                                                                                                                     caGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGtctc             <  4:731021/149‑1 (MQ=255)
                                                                                                                                      aGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGtctc             <  2:45909/148‑1 (MQ=255)
                                                                                                                                       ccaaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATGACGTTATCAAACGCTTACTGGGTTTGATTCCGACGCTGTTTATCGTCTCgg           >  3:177817/3‑149 (MQ=255)
                                                                                                                                         aaaaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGt          <  1:258775/148‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTg         >  3:781070/1‑149 (MQ=255)
                                                                                                                                          aaaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCGATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTg         <  1:584137/148‑1 (MQ=255)
                                                                                                                                          aaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTg         <  1:622456/148‑1 (MQ=255)
                                                                                                                                               cTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtg   <  3:821250/149‑1 (MQ=255)
                                                                                                                                                tGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtgt  <  4:709937/149‑1 (MQ=255)
                                                                                                                                                tGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtgt  >  4:111668/1‑149 (MQ=255)
                                                                                                                                                tGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtgt  <  4:810596/149‑1 (MQ=255)
                                                                                                                                                tGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtg   <  2:658817/148‑1 (MQ=255)
                                                                                                                                                 gACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGgtgt  <  1:514940/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGT  >  NC_000913/870744‑871036

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