Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,159,421 A→T 15.7% I346F (ATT→TTT)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,159,4210AT15.7% 128.1 / 2.7 69I346F (ATT→TTT) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base A (23/31);  new base T (9/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCAT  >  NC_000913/2159273‑2159566
                                                                                                                                                       |                                                                                                                                                 
ggTTGACAGCACACGGGTAAACCGACCGGCGTTGCAGCAAACCACCCCGGCGGCGATTGCGCTGCAAATTGCCCAGGCTCGCTCCC‑CCTCC‑ACATGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGa                                                                                                                                                   <  2:1684971/149‑1 (MQ=255)
          agccGGGGAAAATGACCGGAGTCGCAGGATACCATTCCGGCGGGGATTGATCCGCAAATTGCCCAGGATCGCTCCC‑CCCCC‑TTTCGCGCTTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTggtggt                                                                                                                                         <  2:74214/147‑1 (MQ=255)
           tCCGGGCAAAATTACCGGAGTTGCCGGACACCATTCCGGCGGCTATTGATCTGCAAAATGCCCAGGATCGATCGC‑CGCCC‑ACTAGCGCCTCGGTGGAAGTAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATACTGggggtg                                                                                                                                        >  1:1356120/1‑149 (MQ=255)
           tCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAGGTCGAGCAAACGCTGATTATCTCGGTGG‑CCCCGGTGTTACtgatagag                                                                                                                                        >  2:901589/1‑143 (MQ=255)
           tCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGgtggtg                                                                                                                                        >  2:1671947/1‑149 (MQ=255)
           tCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGggggtg                                                                                                                                        >  1:1344946/1‑149 (MQ=255)
             ggggCAGAA‑TACCGCAGTTGCAGGACACCATTCCGGCGGCTCTTGATCTGCAAATTGCCCAGGAGAGCGCCC‑CCCCC‑CTCCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGttt                                                                                                                                     <  2:36593/148‑1 (MQ=255)
             cGGGCAAGATTACCGGAGTGGCAGGACTCGTTTCCGGGGGCGGTTGAGCTGCCAGTGGCCCAGGGTCGCCCCC‑CCACC‑TTGTGCGCCTCGCTGGATGGAGTCGTGCAGGCGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGtt                                                                                                                                      <  2:842293/149‑1 (MQ=255)
              gggCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGGGgcgtct                                                                                                                                     >  1:728359/1‑147 (MQ=255)
               ggCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGCTTCTGGTGGTGtttt                                                                                                                                    >  1:29641/1‑149 (MQ=255)
                 cAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCGCC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTCTCTCGGTGG‑CGCGGCTGATTCGGGTGGTGgttttg                                                                                                                                  >  1:979783/1‑148 (MQ=255)
                 cAAAATTACCGGAGGTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTa                                                                                                                                  <  1:795371/149‑1 (MQ=255)
                   aaaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCCCC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTAtt                                                                                                                                <  2:100291/149‑1 (MQ=255)
                    caTTACCGGAGTCGCAGGAATCCATCCCGGCGGCGATTGCTGTGCAAATTGCCCAGGATCGCTCGT‑CCACC‑ATTCGCGGCTCGCTGGAATAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTc                                                                                                                               <  2:594901/148‑1 (MQ=255)
                    aagtACCGGAGTTGCTGGAAACCAGTCCGGCGGCGGTTGATGTGCAAATTGCCCGGGATCGCTGCC‑CCACC‑TTTCGGGCGTCGCTGGAAGGAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTc                                                                                                                               <  2:584892/146‑1 (MQ=255)
                    aaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCCGGTGTCTCTGAAGGAGATTATATTc                                                                                                                               >  2:790400/1‑149 (MQ=255)
                        aCCGGGGTTGCATGATACCATTCCGGCGGCTATCGATCCGCAAATTGCACAGGAGCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAGCGCTGGTTATCTCGGTGG‑CGCTGGTCATTCTGGTGGTGTTTTTATTCCTgc                                                                                                                           <  2:966578/149‑1 (MQ=255)
                            gAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCTCCACC‑‑TGCGCGCGTCCCTGGCAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTTTGGCGGTGGTTTTCTTCTTGCACTc                                                                                                                       >  2:655619/1‑149 (MQ=255)
                               ttGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACA‑ATTCGCGCCTCGCTGGACGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGAATCTGGTGGTGTTTGTATTCCTGCGCTCggg                                                                                                                    >  1:509570/1‑149 (MQ=255)
                                tGCACGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATGCGCGCCTCGCTGGAAGAAGGCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGt                                                                                                                   <  1:328015/149‑1 (MQ=255)
                                   aGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTcgc                                                                                                                >  2:1196212/1‑149 (MQ=255)
                                   aGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCCGGTGATCCTGGTGGAGTTTTTATTCCTGCGCTCGGGTcgc                                                                                                                >  2:1232258/1‑149 (MQ=255)
                                     gAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCGTGGTGTCTTTGGTGGTGTTTTTATTCCTGGGCTCGGGTCgcgc                                                                                                              >  2:725575/1‑149 (MQ=255)
                                     gAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAGCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTGTTCCGGCGCTCGGGTggcgg                                                                                                              >  1:285456/1‑148 (MQ=255)
                                      caaCCATCCCGGCGGCGACTGCTCTGCAAATTGCCCAGGATCGCTCCC‑CCGCC‑GTTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGcc                                                                                                             <  1:1758958/148‑1 (MQ=255)
                                               cggcggCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCAACTattatt                                                                                                    >  1:131462/1‑149 (MQ=255)
                                                       aTTGATCTGCAAATTGCCCAGGATCGC‑CCC‑CCACCGATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCCCGTGACTCTGGAAGTGAATTTATTCCTGCCCTCGGGTCGCGACACTAGTATTCCCGCCGt                                                                                            >  2:1748762/1‑149 (MQ=255)
                                                           aTCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGTTGG‑TGCTGGTGAAAATATAGCAAAGCACATCCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCg                                                                                        >  2:1737934/1‑149 (MQ=255)
                                                             cTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCCCTATTATTCCCGCGGTTTCGGt                                                                                      >  1:4215/1‑149 (MQ=255)
                                                               gCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑TGGCGGTGATATTAAAGGAGAATTTATACCTGCGCTAGAGTCGCGCCACTATTATTGCCGACGTTTAGGTGc                                                                                    >  2:1413297/1‑149 (MQ=255)
                                                               gCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGGCGGTGTCTCTGGTGATGATATTCTTCCTGCGCTCGGGTCGCGCCACCATTATTTCCGCCGTTTCGGTga                                                                                    >  2:938298/1‑148 (MQ=255)
                                                                    ttGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCCTCGTATCTCTGGAAGAGATATTATTACTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGtt                                                                               >  2:750472/1‑149 (MQ=255)
                                                                     tGCCCAGGATCGCCCCC‑CCACC‑TTTCGCGCCCCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGttt                                                                              <  1:1431671/149‑1 (MQ=255)
                                                                     tGCCCAGGATCGCCCCC‑CCACC‑AGTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGttt                                                                              <  1:1058515/149‑1 (MQ=255)
                                                                        ccgggctCGCTCCC‑CCACC‑ATGCGCGCCGCGCTGGAAGAAGTCGAGCAACCGCTGATGATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTTGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGc                                                                           <  2:754711/143‑1 (MQ=255)
                                                                          aGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑TCCCGGTGTCATTGGAGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTACTCCCGCCGTGTCGGTGCCGGTTTCGccg                                                                         >  2:1789564/1‑147 (MQ=255)
                                                                            gATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCGTGCGCTCGGGTCGCGCCACTATTATTGCCGCCGTTTCGGTGCCGGTTTCGCTGAt                                                                       >  1:376444/1‑149 (MQ=255)
                                                                             accGCCCCC‑CCACC‑ATTCGCGCCCCGCTGGAAGAAGTCCAGCAAACCCTGATTATCTCGGTGG‑CGCTCGTGATTCTGGTGGGGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGAtt                                                                      <  2:1446818/147‑1 (MQ=255)
                                                                             aTCGCTCCC‑CCAAC‑CTTCGCCCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGAtt                                                                      <  1:838549/149‑1 (MQ=255)
                                                                                  tCCC‑CCGCC‑ACTCGCGCCTCGCGGGAAGGAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTAc                                                                 <  1:680943/149‑1 (MQ=255)
                                                                                  cccc‑ccACC‑ATTCGCCCCTCGCTGGACGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTAc                                                                 <  1:1383487/148‑1 (MQ=255)
                                                                                     c‑ccACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCCGCCGTAGCTT‑TTGAAGAAAAATTATTCCTGCCCTCGGGTCTCGCCACTATTATGCACGCCGTTTCGGTGCCGGTTTCGCGGATTGGTACGtt                                                              >  2:831621/1‑149 (MQ=255)
                                                                                             aTTCGCGCCTGGCCGGAAGAAGCCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgcggc                                                        <  1:981108/149‑1 (MQ=255)
                                                                                              ttCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgcggcg                                                       >  1:1116997/1‑149 (MQ=255)
                                                                                              ttCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCGCCGTATCATTAAAAAAGAAATTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGGTTCGGTGCCGGTTTAGCTGAATGGTACCTTTgcggcg                                                       >  2:598652/1‑149 (MQ=255)
                                                                                                    ccTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTAc                                                 <  1:1671947/149‑1 (MQ=255)
                                                                                                     cTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCGCCGTATCACTGAAAAAAAAATTATTCCTGCGCTCGGGTCGCGCCACTCTTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTCGGTTTGCGGCGACGGAcc                                                >  2:670162/1‑149 (MQ=255)
                                                                                                     cTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCGCCGTAACACTAAAAAAGAAATTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGATTCGCAGATTGGTACGTTTGCGGCGATGTAcc                                                >  2:670366/1‑149 (MQ=255)
                                                                                                     cTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CCCCCGTGACTCTGGAGGAGTTTTTATTCCTGCGCTCGGGTCGAGCCACTATTATCACCGCCGTTTCGCTGCCGGGTTCGCTGATTGGTACGTTTGCGGCGATGTAct                                                >  2:374977/1‑148 (MQ=255)
                                                                                                            gaagaaGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGGGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGGACGTTTGCGGCGATGTACCTGTGCgg                                         >  1:756167/1‑149 (MQ=255)
                                                                                                             aagaagTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGa                                        >  1:1174960/1‑149 (MQ=255)
                                                                                                               gatgTCGGGCCAGCGCTGGTTACCTCGGTGG‑CGCTGGTGCTTCGGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGAtt                                      <  2:1613362/146‑1 (MQ=255)
                                                                                                                   tCGGGCAAACGCTGATTATCTGGGTGG‑CGCTGGTGATTCTGGTGGTGTTCTTATTCCTGCGCTCGGGTCGCGCCACTATTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGt                                  <  1:655619/149‑1 (MQ=255)
                                                                                                                    cgcgCAGACGCCGATTATCGCGGTGG‑CGGTGGTCATTCTGGGGGTGTTTTTATCCCGGCGGCCGGGTCCCGCCCCTATGGTGCCTGCCGTTTCGGTGCCGGTTTCGCGGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGtc                                 <  2:1069988/146‑1 (MQ=255)
                                                                                                                     gAGCAAACGCTGATTATCTCGGTGG‑CGGCGGTATCACTGGAAATGAAATTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGtct                                >  2:1640869/1‑149 (MQ=255)
                                                                                                                       ccgaccGCTGATTATCTCGGTGG‑CGCTGGTGGTTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTACTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCa                              <  2:319542/144‑1 (MQ=255)
                                                                                                                         catcGCTGATTATCTTGGTGG‑CGCGGGTGATTCTGGTGGTGTGTTTATTCCTGCGCTCGGGTGGCGCCCCTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTCCCTGTGCGGATTCAGTCTCAAt                            <  2:1778614/146‑1 (MQ=255)
                                                                                                                          agcGCTGCTTATCTCGGGGG‑CGCTGGCGATTCTGGGGGTGTGTTCATTCCTGCCGTCGGGTCGCGCCACTATTATTCCCGCGGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATa                           <  2:728359/147‑1 (MQ=255)
                                                                                                                            cgtggATTATCTCGGTGG‑CCCTGGGGATTCTGGCGGTGTTTTGACTCCTGCGCTCGGGTTGCGCCACTATTATTCCCGCCGTTTCGGTCCCCGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAAc                         <  2:1577444/145‑1 (MQ=255)
                                                                                                                            cGCTGATTATCTCGGTCG‑CGCTCGTGTTTCTGGTGGTGTCTTCATTCCTCCGCTCGGGTCGCGCCACCATGTTTCCCGCCGTTTCGGTCCCGGTTTCGCCGATCGGTACGCTCGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAAc                         <  2:885322/149‑1 (MQ=255)
                                                                                                                             gCTGATTGACCCGGGGG‑CGCTGGTGATTCTGGTGGTGTTTTTATCCCTGCGCCCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAAcc                        <  2:44879/149‑1 (MQ=255)
                                                                                                                             gCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAAcc                        <  1:1196212/149‑1 (MQ=255)
                                                                                                                               tgctttgctcGGGGG‑CGTTGGTGATTCTGGGGGTGTTTTTATTCCTGCGCTCGGGTCGGGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCtt                      <  2:736280/142‑1 (MQ=255)
                                                                                                                               tGATTAGCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGCCTACATGTACAATCAGTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGCTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCtt                      <  1:1779537/149‑1 (MQ=255)
                                                                                                                                 aTTATCTCGGGCG‑CGCTTGTGATTCTGGGGTTGTTTTTACTCCCCCCCTCGGGTCGCGCTACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGCCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTc                    <  2:826279/149‑1 (MQ=255)
                                                                                                                                   tATCTCGGTGG‑CGCTGGCGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGt                  <  1:289/149‑1 (MQ=255)
                                                                                                                                        cGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATg             <  1:1737934/149‑1 (MQ=255)
                                                                                                                                            gg‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGgcgc         >  2:1120117/1‑149 (MQ=255)
                                                                                                                                              acgcgGGTGTTTCTGGTGGTGTTTTTCTTCCTGCGCTCGGGTCGCGCCACTATTGTTCCCGCCGGTTCGGTGCCGGTTTCGCTGTTTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCt        <  2:1208128/148‑1 (MQ=255)
                                                                                                                                                    gTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACGAt  >  1:1554807/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCC‑CCACC‑ATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGG‑CGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCAT  >  NC_000913/2159273‑2159566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑