Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 268,924 A→C 6.8% S13A (TCG→GCG)  ykfA ← CP4‑6 prophage; putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913268,9240AC 6.8% 527.5 / 24.8 147S13A (TCG→GCG) ykfACP4‑6 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (70/67);  new base C (4/6);  total (74/73)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 148 total aligned reads displayed.

GATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCAATCCTCTGCGC  >  NC_000913/268776‑269071
                                                                                                                                                    |                                                                                                                                                   
gATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGa                                                                                                                                                     <  3:490509/149‑1 (MQ=255)
 atatatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAt                                                                                                                                                    >  2:948680/1‑149 (MQ=255)
      tctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAt                                                                                                                                                    >  3:862588/1‑144 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                              <  3:435934/149‑1 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                              >  4:435934/1‑149 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                              >  4:1466117/1‑149 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                              <  3:1176726/149‑1 (MQ=255)
         ccGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGatc                                                                                                                                            <  4:1170979/149‑3 (MQ=255)
         ccGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                            >  1:932687/1‑149 (MQ=255)
          cGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTg                                                                                                                                                  <  1:231193/142‑1 (MQ=255)
            gCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGa                                                                                                                                         <  2:632313/149‑1 (MQ=255)
             cAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                                        <  2:932687/149‑1 (MQ=255)
                aCAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTc                                                                                                                                     <  1:948680/149‑1 (MQ=255)
                  aGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                            >  3:1450664/1‑140 (MQ=255)
                   gggCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCaaa                                                                                                                                  >  3:1238309/1‑149 (MQ=255)
                   gggCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCaaa                                                                                                                                  >  2:1350293/1‑149 (MQ=255)
                     gCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAAcc                                                                                                                                >  1:1188217/1‑149 (MQ=255)
                       aTTGCACAGACTGCTCTTCCCTACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                                                                                              >  1:1050862/1‑149 (MQ=255)
                       aTTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                                                                                              >  1:841180/1‑149 (MQ=255)
                        ttGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                                                                                             >  2:1325310/1‑149 (MQ=255)
                        ttGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                                                                                             <  2:1188217/149‑1 (MQ=255)
                            acaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                                         >  3:1093546/1‑149 (MQ=255)
                            acaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                                         <  3:901803/149‑1 (MQ=255)
                            acaGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                                         <  1:1350293/149‑1 (MQ=255)
                                aCTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGTTCTGATACACATTGCGGGAGATCAGCAAGCGATTGCTGAAACGACTTGAAACCTTCAGAattat                                                                                                                     >  4:877208/1‑149 (MQ=255)
                                   cctcaTCCCAACACCGGTTATACCCATAATGCCGATCAATGGCTCGTAGTTAGTTAACTGCTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtca                                                                                                                  <  2:1050862/148‑1 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                                                                                 >  1:790695/1‑117 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                                                                                 <  2:790695/117‑1 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                           <  3:223597/123‑1 (MQ=255)
                                     tctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtc                                                                                                                   <  3:821485/146‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                           <  3:1403906/120‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                                        >  1:586220/1‑123 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                                        <  3:675685/123‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                                        >  4:675685/1‑123 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAAc                                                                                                                                 <  3:486585/130‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAAc                                                                                                                                 >  4:486585/1‑130 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                                                                                              >  3:435879/1‑149 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAAAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                           >  4:1403906/1‑120 (MQ=255)
                                          ccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                              >  4:1286418/1‑114 (MQ=255)
                                          ccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                            >  1:1206737/1‑116 (MQ=255)
                                           cAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                                         <  3:123194/134‑1 (MQ=255)
                                           cAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                                         >  4:123194/1‑134 (MQ=255)
                                            aaCACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                                                                                 <  4:1393112/109‑1 (MQ=255)
                                               acCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                            >  4:279741/1‑111 (MQ=255)
                                                    ttttACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtc                                                                                                                   <  1:635281/131‑1 (MQ=255)
                                                       tACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaat                                                                                              <  4:435879/149‑1 (MQ=255)
                                                           cATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaataa                                                                                            >  2:1100427/1‑147 (MQ=255)
                                                                    gATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcatca                                                                                                               <  4:1017310/119‑1 (MQ=255)
                                                                     aTCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGatc                                                                                                                                            >  1:769333/1‑87 (MQ=255)
                                                                     aTCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaata                                                                                             <  2:712803/136‑1 (MQ=255)
                                                                     aTCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaata                                                                                             >  1:712803/1‑136 (MQ=255)
                                                                          tGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACggcggaatc                                                                           >  2:842073/1‑142 (MQ=255)
                                                                             cTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTg                                                                        <  4:652207/149‑1 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcatca                                                                                                               >  1:494791/1‑109 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       >  1:616625/1‑149 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       >  2:308352/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                                     >  2:546387/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                                     >  2:1034697/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                                     >  4:381164/1‑149 (MQ=255)
                                                                                   gttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCa                                                                                                                            <  3:203265/91‑1 (MQ=255)
                                                                                   gttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCa                                                                                                                            >  4:203265/1‑91 (MQ=255)
                                                                                      agttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTccgc                                                               <  1:308352/149‑1 (MQ=255)
                                                                                        ttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaataa                                                                                            >  2:366171/1‑118 (MQ=255)
                                                                                          aaCTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgag                                                                      <  4:488170/138‑1 (MQ=255)
                                                                                             tGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAt                                                        <  1:1034697/149‑1 (MQ=255)
                                                                                             tGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAAt                                                        <  1:546387/149‑1 (MQ=255)
                                                                                               gTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                           >  1:875472/1‑64 (MQ=255)
                                                                                               gTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                           <  2:875472/64‑1 (MQ=255)
                                                                                               gTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                                        <  2:838883/67‑1 (MQ=255)
                                                                                               gTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAattatt                                                                                                                    >  2:569867/1‑87 (MQ=255)
                                                                                                ttAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGc                                                     <  2:1425336/149‑1 (MQ=255)
                                                                                                  aaTTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGc                                                     <  3:163415/147‑1 (MQ=255)
                                                                                                   aTTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                                  <  2:616625/149‑1 (MQ=255)
                                                                                                   aTTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                                  <  4:1093546/149‑1 (MQ=255)
                                                                                                    tttGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                                                                                              >  4:20/1‑72 (MQ=255)
                                                                                                     ttGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       >  4:844193/1‑126 (MQ=255)
                                                                                                       gctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGtt                                              <  3:381164/149‑1 (MQ=255)
                                                                                                        ctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       >  4:132027/1‑123 (MQ=255)
                                                                                                        ctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       <  3:132027/123‑1 (MQ=255)
                                                                                                            aGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAatta                                                                                                                      >  3:1467779/1‑72 (MQ=255)
                                                                                                            aGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAatta                                                                                                                      <  4:1467779/72‑1 (MQ=255)
                                                                                                                gCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaat                                                                                              <  1:1016012/92‑1 (MQ=255)
                                                                                                                           aCCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       >  4:941577/1‑104 (MQ=255)
                                                                                                                           aCCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                                       <  3:941577/104‑1 (MQ=255)
                                                                                                                           aCCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCt                          <  1:1357867/149‑1 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                                                                                    <  3:1114738/87‑1 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                                                                                    >  4:1114738/1‑87 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAg                      >  2:1446462/1‑149 (MQ=255)
                                                                                                                                    ggCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaa                                                                                                  <  3:1342721/68‑1 (MQ=255)
                                                                                                                                      cAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtttttt               >  2:1158261/1‑149 (MQ=255)
                                                                                                                                        gATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCa             >  3:312011/1‑149 (MQ=255)
                                                                                                                                          tCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGc                                                     <  3:574107/107‑1 (MQ=255)
                                                                                                                                            aGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACgg                                                                                  <  3:84872/76‑1 (MQ=255)
                                                                                                                                            aGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACgg                                                                                  >  4:84872/1‑76 (MQ=255)
                                                                                                                                            aGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCaga                                          <  1:543262/116‑1 (MQ=255)
                                                                                                                                              cAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCag                                           >  2:724445/1‑113 (MQ=255)
                                                                                                                                              cAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGtgat                                     >  3:478118/1‑119 (MQ=255)
                                                                                                                                                 gCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCaga                                          >  2:504047/1‑111 (MQ=255)
                                                                                                                                                 gCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCaga                                          <  1:504047/111‑1 (MQ=255)
                                                                                                                                                   gATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCAATCCTCTgcgc  >  2:537275/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCAATCCTCTGCGC  >  NC_000913/268776‑269071

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: ‑