Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 268,921 G→T 6.7% L14I (CTT→ATT)  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,9210GT 6.7% 522.9 / 19.7 147L14I (CTT→ATT) ykfACP4‑6 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (69/68);  new base T (4/6);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 149 total aligned reads displayed.

GATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCA  >  NC_000913/268776‑269060
                                                                                                                                                 |                                                                                                                                           
gATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGa                                                                                                                                          <  3:490509/149‑1 (MQ=255)
 atatatCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGAt                                                                                                                                         >  2:948680/1‑149 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                   <  3:1176726/149‑1 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                   <  3:435934/149‑1 (MQ=255)
       ctcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAATCGCTTGCTGa                                                                                                                                   >  4:435934/1‑149 (MQ=255)
        tcCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaa                                                                                                                                  >  3:652207/1‑149 (MQ=255)
         ccGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                 >  1:932687/1‑149 (MQ=255)
          cGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTg                                                                                                                                       <  1:231193/142‑1 (MQ=255)
          cGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTg                                                                                                                                       >  2:231193/1‑142 (MQ=255)
            gCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGa                                                                                                                              <  2:632313/149‑1 (MQ=255)
             cAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                             <  2:932687/149‑1 (MQ=255)
                aCAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTc                                                                                                                          <  1:948680/149‑1 (MQ=255)
                  aGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                 <  4:1450664/140‑1 (MQ=255)
                  aGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                 >  3:1450664/1‑140 (MQ=255)
                   gggCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCaaa                                                                                                                       >  3:1238309/1‑149 (MQ=255)
                     gCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAAcc                                                                                                                     >  1:1188217/1‑149 (MQ=255)
                       aTTGCACAGACTGCTCTTCCCTACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                                                                                   >  1:1050862/1‑149 (MQ=255)
                       aTTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCt                                                                                                                                     >  3:978444/1‑131 (MQ=255)
                       aTTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCtt                                                                                                                   >  1:841180/1‑149 (MQ=255)
                        ttGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTc                                                                                                                  >  2:1325310/1‑149 (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)
                                   gCTCTTCCCAACACCGGGTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCTATTGCTGAAACGACTTCAAACCTTCAGAATTATtca                                                                                                       <  4:1238309/149‑1 (MQ=255)
                                   cctcaTCCCAACACCGGTTATACCCATAATGCCGATCAATGGCTCGTAGTTAGTTAACTGCTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtca                                                                                                       <  2:1050862/148‑1 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                                                                      <  2:790695/117‑1 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGc                                                                                                                                      >  1:790695/1‑117 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                <  3:223597/123‑1 (MQ=255)
                                    ctctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                >  4:223597/1‑123 (MQ=255)
                                      ctTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAACCGACTTCAAACCTTCAGAATTATtcatca                                                                                                    <  3:877208/149‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                <  3:1403906/120‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                             <  3:675685/123‑1 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                             >  1:586220/1‑123 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGAc                                                                                                                             >  4:675685/1‑123 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAAc                                                                                                                      >  4:486585/1‑130 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCac                                                                                                   >  3:435879/1‑149 (MQ=255)
                                       ttCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAAAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAAc                                                                                                                                >  4:1403906/1‑120 (MQ=255)
                                          ccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                 >  1:1206737/1‑116 (MQ=255)
                                          ccAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGaaa                                                                                                                                 <  2:1206737/116‑1 (MQ=255)
                                           cAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                              <  3:123194/134‑1 (MQ=255)
                                           cAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGaa                                                                                                              >  4:123194/1‑134 (MQ=255)
                                                    ttttACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtc                                                                                                        <  1:635281/131‑1 (MQ=255)
                                                          ccATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGc                                                                                                                                            <  4:452120/89‑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                                                                                  >  1:712803/1‑136 (MQ=255)
                                                                     aTCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAaata                                                                                  <  2:712803/136‑1 (MQ=255)
                                                                          tGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACggcggaatc                                                                >  2:842073/1‑142 (MQ=255)
                                                                             cTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTg                                                             <  4:652207/149‑1 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcatca                                                                                                    <  2:494791/109‑1 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATtcatca                                                                                                    >  1:494791/1‑109 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            >  2:308352/1‑149 (MQ=255)
                                                                              tCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            >  1:616625/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                          >  2:546387/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                          >  4:381164/1‑149 (MQ=255)
                                                                                gTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgaga                                                          >  2:1034697/1‑149 (MQ=255)
                                                                                   gttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCa                                                                                                                 >  4:203265/1‑91 (MQ=255)
                                                                                   gttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCa                                                                                                                 <  3:203265/91‑1 (MQ=255)
                                                                                   gttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTa                                                                     >  4:423311/1‑135 (MQ=255)
                                                                                    ttagttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTcc                                                      >  3:453462/1‑149 (MQ=255)
                                                                                      agttaACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTccgc                                                    <  1:308352/149‑1 (MQ=255)
                                                                                          aaCTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgag                                                           <  4:488170/138‑1 (MQ=255)
                                                                                          aaCTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTgag                                                           >  3:488170/1‑138 (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)
                                                                                                  aaTTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGc                                          >  4:163415/1‑147 (MQ=255)
                                                                                                   aTTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTa                                       <  2:616625/149‑1 (MQ=255)
                                                                                                     ttGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            >  4:844193/1‑126 (MQ=255)
                                                                                                     ttGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            <  3:844193/126‑1 (MQ=255)
                                                                                                       gctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGtt                                   <  3:381164/149‑1 (MQ=255)
                                                                                                        ctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            >  4:132027/1‑123 (MQ=255)
                                                                                                        ctgcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTga                                                            <  3:132027/123‑1 (MQ=255)
                                                                                                          gcAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCag                                >  2:1357867/1‑149 (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)
                                                                                                                           aCCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCt               <  1:1357867/149‑1 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                                                                         <  3:1114738/87‑1 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAAc                                                                         >  4:1114738/1‑87 (MQ=255)
                                                                                                                               aTTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAg           <  4:1131252/149‑1 (MQ=255)
                                                                                                                                    ggCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaa                                                                                       <  3:1342721/68‑1 (MQ=255)
                                                                                                                                    ggCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGaa                                                                                       >  4:1342721/1‑68 (MQ=255)
                                                                                                                                      cAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGAtttttt    >  2:1158261/1‑149 (MQ=255)
                                                                                                                                        gATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCa  >  3:312011/1‑149 (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)
                                                                                                                                                 |                                                                                                                                           
GATATATCTCCGGCAAACAGGGCATTGCACAGACTGCTCTTCCCAACACCGGTTTTACCCATAATGCCGATCACTGGCTCGTAGTTAGTTAACTGGTTAATTTGCTGCAGAAGCCGCTCTGATACCCATTGCGGCAGATCAGCAAGCGATTGCTGAAACGACTTCAAACCTTCAGAATTATTCATCACTACTCCTCTGAAAAATAATAAAAAACGGTAGAATCGTGAGATTCCGCCGTTAATTGCGTATGTTCAGAGTGATGATATATATCTGAAGATTTTTTCA  >  NC_000913/268776‑269060

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: