Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 934,839 T→C V798A (GTT→GCT)  ftsK → cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,8390TC98.4% 532.5 / ‑1.6 185V798A (GTT→GCT) ftsKcell division protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (89/93);  total (91/94)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 211 total aligned reads displayed.

CACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGC  >  NC_000913/934749‑934930
                                                                                          |                                                                                           
cacaCCAGCGGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTg                                                                                           <  1:354665/93‑1 (MQ=255)
cacaACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTg                                                                                           >  2:7657735/1‑93 (MQ=255)
 acaacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgc                                                                                          >  1:5375297/1‑93 (MQ=255)
 acaacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgc                                                                                          >  2:6300746/1‑93 (MQ=255)
   accaGCCGGTGCCGCCCCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgcgc                                                                                        <  2:8317163/91‑1 (MQ=18)
   aacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgcgc                                                                                        >  1:3984571/1‑93 (MQ=255)
      aGCCGGTCGCCCCTCATCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGcgccgc                                                                                     <  2:5721848/93‑1 (MQ=18)
        ccGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcag                                                                                   >  1:2003136/1‑93 (MQ=255)
           gtgccgccccAGCCACAATAGCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcag                                                                                <  2:1919275/89‑1 (MQ=18)
           gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCCACAACCGGTTGGTCCgccgcaa                                                                                >  1:6482433/1‑92 (MQ=18)
           gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcag                                                                                >  1:6482435/1‑93 (MQ=255)
           gTTCCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcag                                                                                <  2:4177813/93‑1 (MQ=255)
            ttGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcagc                                                                               <  1:5282731/93‑1 (MQ=255)
            gtGCGCCGCACCCACAATAGCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcagc                                                                               <  2:3242386/92‑1 (MQ=18)
              gcgcCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGcagcagcag                                                                             >  1:4193378/1‑93 (MQ=255)
              gcgcCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGcagcagcag                                                                             >  1:5717687/1‑93 (MQ=255)
                cccgcAGCCACAATATCAGCAGCCGCAACAACCGCTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTa                                                                           <  1:5550706/92‑1 (MQ=255)
                cccccAGCCACAATACCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTa                                                                           <  1:8840019/89‑1 (MQ=18)
                   gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATca                                                                        <  2:87843/93‑1 (MQ=255)
                      gCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagca                                                                     >  1:1702988/1‑93 (MQ=255)
                       ccACAGTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagcag                                                                    <  2:3948787/93‑1 (MQ=255)
                       ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagcag                                                                    >  2:7101034/1‑93 (MQ=255)
                        cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCagcagc                                                                   >  2:2697948/1‑93 (MQ=255)
                           aaTATCGGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                                <  1:9003068/93‑1 (MQ=255)
                           aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                                >  1:2205984/1‑93 (MQ=255)
                           aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                                >  1:1006089/1‑93 (MQ=255)
                             tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                              <  1:9458177/93‑1 (MQ=255)
                              accagcagCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCacaac                                                             <  2:3397970/91‑1 (MQ=255)
                               acaccagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaaca                                                            <  1:7529329/89‑1 (MQ=255)
                                 agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGc                                                          <  1:9370937/93‑1 (MQ=255)
                                  gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                         >  1:1774213/1‑93 (MQ=255)
                                   cagcCGCAACAACAGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCa                                                        <  2:2117433/93‑1 (MQ=255)
                                    agcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAg                                                       <  1:226979/93‑1 (MQ=255)
                                      ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                     >  1:481179/1‑93 (MQ=255)
                                      ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                     >  1:9840262/1‑93 (MQ=255)
                                      ccGCAACAACCGGTTGCGCAACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                     <  1:4625075/93‑1 (MQ=255)
                                      ccGCAAAAACCGGTTCCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                     <  1:8307544/93‑1 (MQ=255)
                                       cGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTg                                                    >  1:2925618/1‑93 (MQ=255)
                                       cGCAAAAACCGATTGCACCACAGCAGCAATATCGACAGCCGCAACAACCGGCTGGGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTg                                                    <  1:2127760/93‑1 (MQ=14)
                                        gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGttac                                                   >  1:6829740/1‑91 (MQ=255)
                                        gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgc                                                   <  2:5258861/93‑1 (MQ=255)
                                           acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgtc                                                >  2:5536847/1‑91 (MQ=255)
                                           acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGcc                                                >  2:8126839/1‑93 (MQ=255)
                                             aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGATGGGCcac                                              >  1:1640882/1‑93 (MQ=255)
                                              acgggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                             <  1:6592597/90‑1 (MQ=255)
                                              acCGGTTGCGCCACAGCAGCAATATCAGCAGGCGCAACAACCCGCTGCGGCGCAGCAGCAGCATCACCAGCCACAACTGCCATTTGGGCcaca                                             >  2:7752032/1‑93 (MQ=17)
                                               ccGGTGGCGCCACAGCAGCAATATCAGCAGCCGCAACACCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                            <  2:2308922/93‑1 (MQ=255)
                                                cGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAc                                           >  2:2197260/1‑93 (MQ=255)
                                                cGGTTGCGCCACAGCACCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAc                                           <  2:2510529/93‑1 (MQ=255)
                                                 gggtGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAcc                                          <  2:3693114/90‑1 (MQ=255)
                                                  gtggcgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                         <  2:2642482/90‑1 (MQ=255)
                                                  gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                         >  2:9199621/1‑93 (MQ=255)
                                                   ttGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGc                                        >  1:1769794/1‑93 (MQ=255)
                                                    tGCGCCACAGCAGCAATATCAGCAGCCGCAACACCCGGCTCCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                       <  1:704901/93‑1 (MQ=255)
                                                    tGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                       <  1:7782427/93‑1 (MQ=255)
                                                     gcgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAg                                      <  2:6871271/93‑1 (MQ=255)
                                                       gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                    <  1:8783187/93‑1 (MQ=255)
                                                       cccACAACAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                    <  1:1019853/92‑1 (MQ=255)
                                                          acaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAc                                 >  2:4304659/1‑93 (MQ=255)
                                                          acaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAc                                 >  1:9618514/1‑93 (MQ=255)
                                                           cagcacAAATACGAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAcc                                <  2:389746/93‑1 (MQ=21)
                                                            agcagcAATATCAGCAGCCGCAACACCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                               <  2:8874714/93‑1 (MQ=255)
                                                            agcagcAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                               <  1:5969895/93‑1 (MQ=255)
                                                            agcagcAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                               <  1:9505176/93‑1 (MQ=255)
                                                                   atatCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                        >  2:3304187/1‑93 (MQ=255)
                                                                     atCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                      >  1:380089/1‑93 (MQ=255)
                                                                     atCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                      <  2:8500279/93‑1 (MQ=255)
                                                                     atCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                      <  2:4686092/93‑1 (MQ=255)
                                                                      tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                     >  1:9205456/1‑93 (MQ=255)
                                                                      tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                     >  2:9235544/1‑93 (MQ=255)
                                                                      tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                     >  2:3218229/1‑93 (MQ=255)
                                                                      tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                     <  1:5019649/93‑1 (MQ=255)
                                                                       cagcagCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCg                    <  2:6228707/93‑1 (MQ=255)
                                                                        agcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGc                   <  2:10028478/93‑1 (MQ=255)
                                                                        agcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGc                   >  1:3201728/1‑93 (MQ=255)
                                                                        agcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTACGCCACAACCGCAGGATACCCTGCTTCATCCGc                   >  2:125298/1‑93 (MQ=255)
                                                                         gcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCt                  >  2:3476587/1‑93 (MQ=255)
                                                                          cagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTg                 <  1:8342492/93‑1 (MQ=255)
                                                                             ccGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg              <  1:5148366/93‑1 (MQ=255)
                                                                             ccGCAACAACCGGCTCCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg              <  2:2524603/93‑1 (MQ=255)
                                                                              cGCAACACCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa             <  2:1192385/93‑1 (MQ=255)
                                                                              cGCAACAACCGTCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa             <  1:3134177/93‑1 (MQ=255)
                                                                              cGCAACAACCGGCTTCCCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa             <  2:4396160/93‑1 (MQ=255)
                                                                              cGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa             <  1:7151534/93‑1 (MQ=255)
                                                                              cGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa             <  2:2994205/93‑1 (MQ=255)
                                                                               gCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt            >  2:2414354/1‑93 (MQ=255)
                                                                               gCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt            >  2:80003/1‑93 (MQ=255)
                                                                               gCAACAACCGGATGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt            <  2:9946138/93‑1 (MQ=255)
                                                                               caaacaacCGGCTGCGCCGCCGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt            <  1:4341923/91‑1 (MQ=255)
                                                                                  acaacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCg         >  1:471214/1‑93 (MQ=255)
                                                                                    aacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa       <  2:6697669/93‑1 (MQ=255)
                                                                                    aacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa       >  1:7963557/1‑93 (MQ=255)
                                                                                    aacCGGCGGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa       <  1:2828166/93‑1 (MQ=255)
                                                                                     acCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTaa      >  2:7489529/1‑93 (MQ=255)
                                                                                     acCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTaa      <  2:5989437/93‑1 (MQ=255)
                                                                                      ccGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt     <  1:5867115/93‑1 (MQ=255)
                                                                                      ccGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt     <  1:8703293/93‑1 (MQ=255)
                                                                                      acGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt     <  1:4171020/92‑1 (MQ=255)
                                                                                       cGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg    <  2:6092990/93‑1 (MQ=255)
                                                                                         gTTCCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATCCCCTGCTTCATCCGCTGTTGATGCGTAATGGc  <  2:6816639/93‑1 (MQ=255)
                                                                                          |                                                                                           
CACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGC  >  NC_000913/934749‑934930

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