Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,487 T→A 5.4% T73S (ACT→TCT) ‡ hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,4870TA 5.4% 370.5 / 9.2 113T73S (ACT→TCT) ‡hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (48/56);  new base A (4/2);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 117 total aligned reads displayed.

TACTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAC  >  NC_000913/3520340‑3520634
                                                                                                                                                   |                                                                                                                                                   
   tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  3:7475/1‑149 (MQ=255)
   tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  3:864485/1‑149 (MQ=255)
    gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  2:826869/1‑148 (MQ=255)
    gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  1:506438/1‑148 (MQ=255)
    gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  2:294658/1‑148 (MQ=255)
    gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  2:163116/1‑148 (MQ=255)
    gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTa                                                                                                                                                <  3:204720/149‑1 (MQ=255)
     cGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAg                                                                                                                                               <  4:209119/149‑1 (MQ=255)
       aTTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGt                                                                                                                                              >  2:276676/1‑148 (MQ=255)
        ttttCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                            <  3:196367/149‑1 (MQ=255)
             cAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCt                                                                                                                                       >  3:804336/1‑149 (MQ=255)
               gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCtt                                                                                                                                      <  1:502310/148‑1 (MQ=255)
                  gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                                   <  2:447940/148‑1 (MQ=255)
                  gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                  <  4:389234/149‑1 (MQ=255)
                  gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                  <  3:499051/149‑1 (MQ=255)
                      aTTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTAGCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGa                                                                                                                               <  2:865239/148‑1 (MQ=255)
                             tGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                                                                                                                        <  2:252070/148‑1 (MQ=255)
                             tGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                                                                                                                        >  2:704259/1‑148 (MQ=255)
                               cgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      <  1:27630/148‑1 (MQ=255)
                               cgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      <  2:569689/148‑1 (MQ=255)
                                  cgatctcTGCACCGCCTCCTGGCGGGCGATATTGTTATACTGCCAGGCCATGGAATGCACTTATAGAATGGTGCCTTGCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGGGCGTAGTGCCTGCTTCCAAGGAACATCTGTTAGAt                                                                                                                  <  4:823209/143‑1 (MQ=255)
                                  ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                                                   >  1:815051/1‑148 (MQ=255)
                                  ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                                                   >  1:747447/1‑148 (MQ=255)
                                  ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGAt                                                                                                                  <  4:672048/149‑1 (MQ=255)
                                   cgcgCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGAt                                                                                                                  >  1:585535/1‑148 (MQ=255)
                                      gcCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTaa                                                                                                              <  3:491818/149‑1 (MQ=255)
                                        cTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      <  3:683571/139‑1 (MQ=255)
                                         tGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATgagatcg                                                                                                           >  3:823209/1‑143 (MQ=255)
                                             ccGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac                                                                                                        >  2:210237/1‑148 (MQ=255)
                                               gCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGt                                                                                                     >  4:225212/1‑149 (MQ=255)
                                                cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGt                                                                                                     <  2:231377/148‑1 (MQ=255)
                                                cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGt                                                                                                     <  1:276676/148‑1 (MQ=255)
                                                cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTa                                                                                                    <  4:896837/149‑1 (MQ=255)
                                                  tcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTAcc                                                                                                  >  4:485180/1‑149 (MQ=255)
                                                  tcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGCTTGTAGTGCCTGCTTCCAGGGAACATCTGt                                                                                                                       >  4:649149/1‑128 (MQ=255)
                                                  tcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGCTTGTAGTGCCTGCTTCCAGGGAACATCTGt                                                                                                                       <  3:649149/128‑1 (MQ=255)
                                                    ctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                  <  3:900143/115‑1 (MQ=255)
                                                      ggcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              >  4:432733/1‑149 (MQ=255)
                                                        cgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGa                                                                                             <  2:272418/148‑1 (MQ=255)
                                                            ctatatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTc                                                                                         <  1:477564/146‑1 (MQ=255)
                                                               tatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                                     <  3:211669/149‑1 (MQ=255)
                                                                atTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                            >  4:326508/1‑93 (MQ=255)
                                                                atTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                            <  3:326508/93‑1 (MQ=255)
                                                                    ttATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGaca                                                                                 <  2:96148/148‑1 (MQ=255)
                                                                      aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              <  4:7475/149‑1 (MQ=255)
                                                                      aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              >  3:829349/1‑149 (MQ=255)
                                                                      aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              >  4:615262/1‑149 (MQ=255)
                                                                       ttCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                                             >  3:567502/1‑149 (MQ=255)
                                                                         cTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACGCCGTACCGCTGACGTCTGGCGACACGAcc                                                                            <  1:163116/148‑1 (MQ=255)
                                                                         cTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAcc                                                                            <  1:915127/148‑1 (MQ=255)
                                                                                        aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                                             <  1:474433/148‑1 (MQ=255)
                                                                                        aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCt                                                            <  3:303944/149‑1 (MQ=255)
                                                                                              tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                                                   <  4:355668/120‑1 (MQ=255)
                                                                                              tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                                                   >  3:355668/1‑120 (MQ=255)
                                                                                                agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                   <  4:70034/134‑1 (MQ=255)
                                                                                                agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                   >  3:70034/1‑134 (MQ=255)
                                                                                                 gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              <  3:223542/106‑1 (MQ=255)
                                                                                                 gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              >  4:223542/1‑106 (MQ=255)
                                                                                                   gaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagccagc                                                 >  4:488586/1‑149 (MQ=255)
                                                                                                     aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtct                                                               >  1:941714/1‑133 (MQ=255)
                                                                                                     aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtct                                                               <  2:941714/133‑1 (MQ=255)
                                                                                                      tGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGCGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGTTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                                              >  3:484889/1‑149 (MQ=255)
                                                                                                        ttGCCTTCCTGCCGCGTTATCAGTCCGCCGCTTTTCACTACAAACTGTAGt                                                                                                                                              <  4:62410/51‑1 (MQ=255)
                                                                                                        ttGCCTTCCTACCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGt                                                                                                                                              >  3:62410/1‑51 (MQ=255)
                                                                                                         tGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                                           <  4:567502/149‑1 (MQ=255)
                                                                                                              tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                                                       <  2:852204/116‑1 (MQ=255)
                                                                                                              tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                                                       >  1:852204/1‑116 (MQ=255)
                                                                                                              tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                                        >  1:893851/1‑147 (MQ=255)
                                                                                                              tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                                        <  2:893851/147‑1 (MQ=255)
                                                                                                                cTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTGCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAcc                                     <  2:585535/148‑1 (MQ=255)
                                                                                                                 tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                                    >  1:736813/1‑148 (MQ=255)
                                                                                                                 tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                                    >  2:296580/1‑148 (MQ=255)
                                                                                                                  gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                                   <  2:930846/148‑1 (MQ=255)
                                                                                                                  gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                                   <  2:23060/148‑1 (MQ=255)
                                                                                                                  gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                                  <  4:804336/149‑1 (MQ=255)
                                                                                                                    cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg                                 >  1:144948/1‑148 (MQ=255)
                                                                                                                       gTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCAGCTTCCATGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCCGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             <  4:596045/149‑1 (MQ=255)
                                                                                                                        ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAcc                            <  3:285102/149‑1 (MQ=255)
                                                                                                                        ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             >  2:91149/1‑148 (MQ=255)
                                                                                                                        ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             >  2:828840/1‑148 (MQ=255)
                                                                                                                         tATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAcc                            >  2:613851/1‑148 (MQ=255)
                                                                                                                             aGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACGTGAGCTACCGGAAc                       <  3:74557/149‑1 (MQ=255)
                                                                                                                             aGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGaa                        <  2:641590/148‑1 (MQ=255)
                                                                                                                                ccGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGt                     <  1:828840/148‑1 (MQ=255)
                                                                                                                                  ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCAt                  >  4:274844/1‑149 (MQ=255)
                                                                                                                                   gcgcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCAt                  <  1:296580/148‑1 (MQ=255)
                                                                                                                                      cTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCa                                                                           >  1:813164/1‑88 (MQ=255)
                                                                                                                                      cTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCa                                                                           <  2:813164/88‑1 (MQ=255)
                                                                                                                                         ggcACTACAAACTGTAGTGCCTGCTTCCAGGGACCATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATccacca            >  1:577014/3‑148 (MQ=255)
                                                                                                                                          gcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              >  2:793366/2‑97 (MQ=255)
                                                                                                                                          gcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              <  1:793366/96‑1 (MQ=255)
                                                                                                                                           cACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATc          >  2:286296/1‑148 (MQ=255)
                                                                                                                                             cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGc       >  3:780664/1‑149 (MQ=255)
                                                                                                                                             cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGc       >  4:408879/1‑149 (MQ=255)
                                                                                                                                               aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCg      >  2:104748/1‑148 (MQ=255)
                                                                                                                                               aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCg      <  1:706182/148‑1 (MQ=255)
                                                                                                                                                cAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGaa                        >  2:159336/1‑129 (MQ=255)
                                                                                                                                                cAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGaa                        <  1:159336/129‑1 (MQ=255)
                                                                                                                                                   tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              <  4:312016/88‑1 (MQ=255)
                                                                                                                                                   tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAc  >  2:198346/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAC  >  NC_000913/3520340‑3520634

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