Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,693 G→T 7.5% T395K (ACG→AAG)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,6930GT 7.5% 366.2 / 14.8 105T395K (ACG→AAG) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (49/48);  new base T (4/4);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 106 total aligned reads displayed.

ACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGC  >  NC_000913/2318546‑2318832
                                                                                                                                                   |                                                                                                                                           
accaGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCg                                                                                                                                             >  2:554112/1‑148 (MQ=255)
   aGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGt                                                                                                                                            <  4:6698/146‑1 (MQ=255)
   aGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGt                                                                                                                                            >  3:6698/1‑146 (MQ=255)
    gCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                                                                                                                         <  1:554112/148‑1 (MQ=255)
                ccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                                                                                                                             >  2:674946/1‑148 (MQ=255)
                  aCGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGc                                                                                                                                        <  4:7613/135‑1 (MQ=255)
                  aCGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                >  1:95145/1‑143 (MQ=255)
                  aCGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                <  2:95145/143‑1 (MQ=255)
                     ttttCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                                                                                                                         >  4:47403/1‑131 (MQ=255)
                     ttttCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                                                                                                                         <  3:47403/131‑1 (MQ=255)
                           ttACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCat                                                                                                                                      <  2:517479/128‑1 (MQ=255)
                           ttACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCat                                                                                                                                      >  1:517479/1‑128 (MQ=255)
                            tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTa                                                                                                                                           >  2:388291/1‑122 (MQ=255)
                            tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTa                                                                                                                                           <  1:388291/122‑1 (MQ=255)
                            tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                                                                                                                             >  3:708485/1‑136 (MQ=255)
                            tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                                                                                                                             <  4:708485/136‑1 (MQ=255)
                            tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAg                                                                                                                 >  4:495448/1‑148 (MQ=255)
                                      cgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTa                                                                                                                         <  2:265001/130‑1 (MQ=255)
                                      cgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTa                                                                                                                         >  1:265001/1‑130 (MQ=255)
                                        aaTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGAt                                                                                                                              >  1:546466/1‑123 (MQ=255)
                                        aaTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGAt                                                                                                                              <  2:546466/123‑1 (MQ=255)
                                              cGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatg                                                                                                                                     >  3:1036694/1‑109 (MQ=255)
                                              cGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatg                                                                                                                                     <  4:1036694/110‑2 (MQ=255)
                                                aaGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGa                                                                                             <  3:495448/148‑1 (MQ=255)
                                                     gATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCa                                                                                                                                       <  4:937962/101‑1 (MQ=255)
                                                     gATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCa                                                                                                                                       >  3:937962/1‑101 (MQ=255)
                                                      aTTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCatat                                                                                                                                    >  2:381406/1‑103 (MQ=255)
                                                      aTTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCatat                                                                                                                                    <  1:381406/103‑1 (MQ=255)
                                                             gATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                    >  3:535211/1‑128 (MQ=255)
                                                             gATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                    <  4:535211/128‑1 (MQ=255)
                                                                tGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGTTAGGTATGCGCCAGTTCGtt                                                                                                           >  2:83740/1‑118 (MQ=255)
                                                                tGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGtt                                                                                                           <  1:83740/118‑1 (MQ=255)
                                                                   tATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                    <  3:1065556/122‑1 (MQ=255)
                                                                   tATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaa                                                                                                    >  4:1065556/1‑122 (MQ=255)
                                                                     ttgttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa                                                                                                      >  4:826463/1‑118 (MQ=255)
                                                                     ttgttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa                                                                                                      <  3:826463/118‑1 (MQ=255)
                                                                     ttgttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCaaaa                                                                        <  1:674946/148‑1 (MQ=255)
                                                                        ttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                    >  2:945490/1‑83 (MQ=255)
                                                                        ttgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                    <  1:945490/85‑3 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatga                                                                                                                                    <  3:427202/84‑3 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                <  4:509584/88‑1 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                <  1:222461/88‑1 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                >  2:222461/1‑88 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                                                                                                                >  3:509584/1‑88 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGcc                                                                    >  2:986806/1‑147 (MQ=255)
                                                                         tgCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCa                                                                    <  1:986806/148‑1 (MQ=255)
                                                                           cTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                               >  4:241396/1‑135 (MQ=255)
                                                                           cTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                               <  3:241396/135‑1 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcgca                                                                                                                   <  4:1023351/97‑2 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     <  1:227397/143‑1 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     >  2:227397/1‑143 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                >  4:1023107/1‑148 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                >  2:131473/1‑148 (MQ=255)
                                                                             ggTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                >  4:774222/1‑148 (MQ=255)
                                                                               tCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     <  4:537493/141‑1 (MQ=255)
                                                                               tCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     >  3:537493/1‑141 (MQ=255)
                                                                                cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACt                                                                                                        >  2:799324/1‑105 (MQ=255)
                                                                                cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     <  1:204291/140‑1 (MQ=255)
                                                                                cAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     >  2:204291/1‑140 (MQ=255)
                                                                                  ggACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                                                                                                                             >  4:544460/1‑82 (MQ=255)
                                                                                  ggACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                                                                                                                             <  3:544460/82‑1 (MQ=255)
                                                                                    aCATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                             <  3:293691/96‑1 (MQ=255)
                                                                                    aCATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                             >  4:293691/1‑96 (MQ=255)
                                                                                      aTCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTAATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGtt                                                                                                           >  1:234444/1‑96 (MQ=255)
                                                                                      aTCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTAATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGtt                                                                                                           <  2:234444/96‑1 (MQ=255)
                                                                                        caacaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                <  1:488982/137‑1 (MQ=255)
                                                                                          acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTg                                                                                                                                  >  3:792802/1‑69 (MQ=255)
                                                                                          acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTg                                                                                                                                  <  4:792802/69‑1 (MQ=255)
                                                                                          acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                             >  2:714198/1‑90 (MQ=255)
                                                                                          acaaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                             <  1:714198/90‑1 (MQ=255)
                                                                                           caaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATg                                                                                                                       >  1:972613/1‑79 (MQ=255)
                                                                                           caaAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATg                                                                                                                       <  2:972613/79‑1 (MQ=255)
                                                                                            aaaaTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                 <  3:1023107/148‑1 (MQ=255)
                                                                                            aaaaTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                 <  3:774222/148‑1 (MQ=255)
                                                                                               agtGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcgca                                                                                                                   <  4:819796/77‑2 (MQ=255)
                                                                                               aTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGcc                                                                                                                   >  3:819796/1‑79 (MQ=255)
                                                                                               aTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGc                                                                     >  3:9747/1‑125 (MQ=255)
                                                                                                ttGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCa                                                            <  3:252564/133‑1 (MQ=255)
                                                                                                ttGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCa                                                            >  4:252564/1‑133 (MQ=255)
                                                                                                ttGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGa                                             >  3:215806/1‑148 (MQ=255)
                                                                                                       gctgcCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCtt                                      <  1:131473/148‑1 (MQ=255)
                                                                                                             cacaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGt                                                                                                            >  4:257071/1‑72 (MQ=255)
                                                                                                             cacaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGt                                                                                                            <  3:257071/72‑1 (MQ=255)
                                                                                                              acaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTg                                     >  2:886470/1‑142 (MQ=255)
                                                                                                              acaGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTg                                     <  1:886470/142‑1 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTaa                                                                                                                                          >  4:1003310/1‑40 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTaa                                                                                                                                          <  3:1003310/40‑1 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCgg                                                                                 <  2:893483/97‑1 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCgg                                                                                 >  1:893483/1‑97 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCg                                              <  3:350196/132‑1 (MQ=255)
                                                                                                               caGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCg                                              >  4:350196/1‑132 (MQ=255)
                                                                                                                aGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                               <  3:93392/98‑1 (MQ=255)
                                                                                                                aGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTa                                                                               >  4:93392/1‑98 (MQ=255)
                                                                                                                         gCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGc                                    <  2:262035/132‑1 (MQ=255)
                                                                                                                         gCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGc                                    >  1:262035/1‑132 (MQ=255)
                                                                                                                             cAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCa                                                            <  2:708484/104‑1 (MQ=255)
                                                                                                                             cAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCa                                                            >  1:708484/1‑104 (MQ=255)
                                                                                                                               gTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                 >  2:1026817/1‑113 (MQ=255)
                                                                                                                               gTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAg                                                 <  1:1026817/113‑1 (MQ=255)
                                                                                                                                           ccTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGc  >  4:602454/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
ACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGC  >  NC_000913/2318546‑2318832

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