Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,488 T→C 5.3% Q72Q (CAA→CAG 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,4880TC 5.3% 362.8 / 6.2 113Q72Q (CAA→CAGhofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (49/57);  new base C (4/2);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 116 total aligned reads displayed.

TGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAC  >  NC_000913/3520343‑3520634
                                                                                                                                                 |                                                                                                                                                  
tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  3:7475/1‑149 (MQ=255)
 gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  2:294658/1‑148 (MQ=255)
 gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                 >  1:506438/1‑148 (MQ=255)
 gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTa                                                                                                                                                <  3:204720/149‑1 (MQ=255)
  cGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAg                                                                                                                                               <  4:209119/149‑1 (MQ=255)
     ttttCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                            <  3:196367/149‑1 (MQ=255)
          cAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCt                                                                                                                                       >  3:804336/1‑149 (MQ=255)
            gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCtt                                                                                                                                      <  1:502310/148‑1 (MQ=255)
               gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                                   <  1:61300/148‑1 (MQ=255)
               gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                                   <  2:447940/148‑1 (MQ=255)
               gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                  <  4:389234/149‑1 (MQ=255)
                cAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAgg                                                                                                                                 <  4:349447/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                                                                                                                      <  2:569689/148‑1 (MQ=255)
                            cgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      <  1:27630/148‑1 (MQ=255)
                              cccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                                                                                                                    <  2:385082/148‑1 (MQ=255)
                              cccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                                                   <  3:774618/149‑1 (MQ=255)
                               cgatctcTGCACCGCCTCCTGGCGGGCGATATTGTTATACTGCCAGGCCATGGAATGCACTTATAGAATGGTGCCTTGCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGGGCGTAGTGCCTGCTTCCAAGGAACATCTGTTAGAt                                                                                                                  <  4:823209/143‑1 (MQ=255)
                               ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                                                   >  1:815051/1‑148 (MQ=255)
                               ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                                                   >  1:747447/1‑148 (MQ=255)
                                cgcgCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGAt                                                                                                                  >  1:585535/1‑148 (MQ=255)
                                   gcCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTaa                                                                                                              <  3:491818/149‑1 (MQ=255)
                                     cTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      <  3:683571/139‑1 (MQ=255)
                                     cTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                      >  4:683571/1‑139 (MQ=255)
                                     cTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACg                                                                                                            <  3:618063/149‑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                                                                                                                       <  3:649149/128‑1 (MQ=255)
                                               tcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGCTTGTAGTGCCTGCTTCCAGGGAACATCTGt                                                                                                                       >  4:649149/1‑128 (MQ=255)
                                                cctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCTTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTAcc                                                                                                  >  2:122698/1‑148 (MQ=255)
                                                 ctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                  <  3:900143/115‑1 (MQ=255)
                                                   ggcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              >  4:432733/1‑149 (MQ=255)
                                                     cgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGa                                                                                             <  2:272418/148‑1 (MQ=255)
                                                          gATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                                                       >  4:526640/1‑149 (MQ=255)
                                                            tatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                                     <  3:211669/149‑1 (MQ=255)
                                                             atTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                            <  3:326508/93‑1 (MQ=255)
                                                                 ttATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGaca                                                                                 <  2:96148/148‑1 (MQ=255)
                                                                   aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              >  4:615262/1‑149 (MQ=255)
                                                                   aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              >  3:829349/1‑149 (MQ=255)
                                                                   aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                                              <  4:7475/149‑1 (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                                                                                   >  3:355668/1‑120 (MQ=255)
                                                                                           tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                                                   <  4:355668/120‑1 (MQ=255)
                                                                                            gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                                                      >  2:891716/1‑148 (MQ=255)
                                                                                             agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                   >  3:70034/1‑134 (MQ=255)
                                                                                             agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                   <  4:70034/134‑1 (MQ=255)
                                                                                              gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              <  3:223542/106‑1 (MQ=255)
                                                                                              gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                              >  4:223542/1‑106 (MQ=255)
                                                                                                  aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtct                                                               <  2:941714/133‑1 (MQ=255)
                                                                                                  aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtct                                                               >  1:941714/1‑133 (MQ=255)
                                                                                                   tGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGCGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGTTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                                              >  3:484889/1‑149 (MQ=255)
                                                                                                    gTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                                              <  1:826869/148‑1 (MQ=255)
                                                                                                    gTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                                              <  2:239994/148‑1 (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)
                                                                                                           tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                                        >  1:893851/1‑147 (MQ=255)
                                                                                                           tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                                        <  2:893851/147‑1 (MQ=255)
                                                                                                              tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                                    >  2:296580/1‑148 (MQ=255)
                                                                                                               gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                                   <  2:23060/148‑1 (MQ=255)
                                                                                                               gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                                  <  3:169151/149‑1 (MQ=255)
                                                                                                               gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                                  <  4:804336/149‑1 (MQ=255)
                                                                                                               gCCGCGTTATCAGCCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                                   >  2:774942/1‑148 (MQ=255)
                                                                                                                 cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg                                 >  1:144948/1‑148 (MQ=255)
                                                                                                                    gTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCAGCTTCCATGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCCGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             <  4:596045/149‑1 (MQ=255)
                                                                                                                     ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             >  2:828840/1‑148 (MQ=255)
                                                                                                                     ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                             >  2:91149/1‑148 (MQ=255)
                                                                                                                      tATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAcc                            >  2:613851/1‑148 (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                                                                           <  2:813164/88‑1 (MQ=255)
                                                                                                                                   cTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCa                                                                           >  1:813164/1‑88 (MQ=255)
                                                                                                                                      ggcACTACAAACTGTAGTGCCTGCTTCCAGGGACCATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATccacca            >  1:577014/3‑148 (MQ=255)
                                                                                                                                       gcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              <  1:793366/96‑1 (MQ=255)
                                                                                                                                       gcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              >  2:793366/2‑97 (MQ=255)
                                                                                                                                        cACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATc          >  2:286296/1‑148 (MQ=255)
                                                                                                                                          cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGc       >  4:408879/1‑149 (MQ=255)
                                                                                                                                          cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGc       >  3:780664/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)
                                                                                                                                               gTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGTCGACACGACCTGGTTCAACTTCTCCTGTGCAGCCAGCGCCTGCAACACTTGAGCTACCGGAACGTCATCCACCATCAGCGTCAc  >  3:945199/1‑149 (MQ=255)
                                                                                                                                                tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              >  3:312016/1‑88 (MQ=255)
                                                                                                                                                tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                              <  4:312016/88‑1 (MQ=255)
                                                                                                                                                tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAc  >  2:198346/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAC  >  NC_000913/3520343‑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: