Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 7.7% intergenic (+28/+75) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3150CG 7.7% 400.7 / 20.5 119intergenic (+28/+75)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (56/53);  new base G (2/7);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 120 total aligned reads displayed.

ACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACT  >  NC_000913/2286167‑2286462
                                                                                                                                                     |                                                                                                                                                   
aCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                     <  4:351742/149‑1 (MQ=255)
    gccaTCGGTGTTGAAAACGGGAATTATAGGCGGCAACAGTCGGCAACGAGCGCCACCTTGTTGCCGCACCGCAAGGGGTGGCCGGGGGCGGAAGTTCAAATCCTTTCGTGGCCACCAAAAAACCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                 <  4:383194/146‑1 (MQ=255)
    gCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                 >  3:13099/1‑149 (MQ=255)
      aaTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                                                                               >  1:562789/1‑149 (MQ=255)
       aTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGTGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                 >  3:474509/1‑146 (MQ=255)
       aTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                 <  4:474509/146‑1 (MQ=255)
        tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                                                                             <  4:771779/149‑1 (MQ=255)
        tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                                                                             <  1:469083/149‑1 (MQ=255)
        tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggtt                                                                                                                                             >  4:187806/1‑149 (MQ=255)
         cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                                                                            <  1:212630/149‑1 (MQ=255)
          ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                           <  2:848232/149‑1 (MQ=255)
          ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                           <  2:994846/149‑1 (MQ=255)
            tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         <  1:1013926/149‑1 (MQ=255)
            tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         <  3:382776/149‑1 (MQ=255)
            tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         <  3:931335/149‑1 (MQ=255)
            tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         <  4:232094/149‑1 (MQ=255)
             gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttc                                                                                                                                        <  2:192210/149‑2 (MQ=255)
               ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                                                                      <  3:629108/149‑1 (MQ=255)
                gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                     >  1:143162/1‑149 (MQ=255)
                   aaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                                                                                  >  2:175217/1‑149 (MQ=255)
                       gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                                                                              <  3:462627/149‑1 (MQ=255)
                        tagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                                                                             >  3:192621/1‑149 (MQ=255)
                           taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaa                                                                                                                          >  1:931835/1‑149 (MQ=255)
                           taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaactaa                                                                                                                          <  4:523246/149‑1 (MQ=255)
                           taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGGTTTTTTATATCTGCaattaa                                                                                                                          >  2:226141/1‑149 (MQ=255)
                            aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                                                                         <  2:669700/149‑1 (MQ=255)
                             ttAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                                                                        >  4:97698/1‑149 (MQ=255)
                              tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACTCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                       <  2:568598/149‑1 (MQ=255)
                              tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                       <  2:748835/149‑1 (MQ=255)
                              tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                       <  1:868363/149‑1 (MQ=255)
                               agcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                      >  4:834061/1‑149 (MQ=255)
                                gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTTATATCTGCAATTAATTCGa                                                                                                                     >  2:312661/1‑149 (MQ=255)
                                gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGa                                                                                                                     >  4:883513/1‑149 (MQ=255)
                                    ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                 <  2:253563/149‑1 (MQ=255)
                                             ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                                                                        >  2:937587/1‑149 (MQ=255)
                                                aCGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGaca                                                                                                     >  3:800369/1‑149 (MQ=255)
                                                   tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                                                                  <  1:1084682/149‑1 (MQ=255)
                                                   tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                                                                  <  2:734649/149‑1 (MQ=255)
                                                       gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                                                                              >  1:89960/1‑149 (MQ=255)
                                                           ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGTAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                                          <  2:305310/149‑1 (MQ=255)
                                                             tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                        <  1:312661/149‑1 (MQ=255)
                                                             tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                        <  1:1059508/149‑1 (MQ=255)
                                                              ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                                                                       <  3:466152/149‑1 (MQ=255)
                                                               gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAAGTACAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACTCTTATGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttg                                                                                      >  2:788245/1‑148 (MQ=255)
                                                                tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                                                                     >  2:731451/1‑149 (MQ=255)
                                                                        ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                                <  2:92246/147‑1 (MQ=255)
                                                                        ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                                >  1:92246/1‑147 (MQ=255)
                                                                         cGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                             <  3:623456/149‑1 (MQ=255)
                                                                             aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                         <  3:834061/149‑1 (MQ=255)
                                                                                  gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                                                                    <  3:950472/149‑1 (MQ=255)
                                                                                    gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                                                  <  4:18731/133‑1 (MQ=255)
                                                                                    gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                                                  >  3:18731/1‑133 (MQ=255)
                                                                                     tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc                                                                 <  4:451027/149‑1 (MQ=255)
                                                                                     tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc                                                                 <  3:157505/149‑1 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                 <  2:1017948/92‑1 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                 >  1:1017948/1‑92 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                                                                       >  3:867652/1‑102 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                                                                       <  4:867652/102‑1 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaca                                                        <  1:731451/149‑1 (MQ=255)
                                                                                                aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAcc                                                      <  1:607396/149‑1 (MQ=255)
                                                                                                 ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc                                                     >  4:28911/1‑149 (MQ=255)
                                                                                                  gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                    >  1:266752/1‑149 (MQ=255)
                                                                                                  gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                    >  4:552908/1‑149 (MQ=255)
                                                                                                  gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                    >  1:925242/1‑149 (MQ=255)
                                                                                                  gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGCTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                    >  2:584233/1‑149 (MQ=255)
                                                                                                    tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCATCCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                                  <  2:269366/149‑1 (MQ=255)
                                                                                                     cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt                                                 >  1:732853/1‑149 (MQ=255)
                                                                                                            tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                          >  1:240157/1‑149 (MQ=255)
                                                                                                            tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                          >  4:522431/1‑149 (MQ=255)
                                                                                                            tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCAGCCTGAGTAtt                                          >  3:292873/1‑149 (MQ=255)
                                                                                                             ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTc                                         <  1:821823/149‑1 (MQ=255)
                                                                                                              tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGCATTcc                                        >  1:388964/1‑149 (MQ=255)
                                                                                                                  gCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                    >  3:218994/1‑149 (MQ=255)
                                                                                                                     gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGACTATTCCACGGCGa                                 >  4:684447/1‑149 (MQ=255)
                                                                                                                       ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAtt                               >  4:684112/1‑149 (MQ=255)
                                                                                                                       ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAtt                               <  1:1024036/149‑1 (MQ=255)
                                                                                                                        cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTc                              <  2:89960/149‑1 (MQ=255)
                                                                                                                        cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTc                              >  4:884230/1‑149 (MQ=255)
                                                                                                                         aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                    <  3:759351/142‑1 (MQ=255)
                                                                                                                         aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                    >  4:759351/1‑142 (MQ=255)
                                                                                                                          aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTa                            >  2:403281/1‑149 (MQ=255)
                                                                                                                           aaaTCGCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCTGCCTGAGTATTCCGCGGCGATTCTAc                           >  4:670531/1‑149 (MQ=255)
                                                                                                                            aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAcc                          >  1:55819/1‑149 (MQ=255)
                                                                                                                             aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAccc                         <  3:552908/149‑1 (MQ=255)
                                                                                                                                ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                         <  4:539171/98‑1 (MQ=255)
                                                                                                                                ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                         >  3:539171/1‑98 (MQ=255)
                                                                                                                                 cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                                                                         <  1:980889/65‑1 (MQ=255)
                                                                                                                                 cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTc                                         >  1:82067/1‑129 (MQ=255)
                                                                                                                                 cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTc                                         <  2:82067/129‑1 (MQ=255)
                                                                                                                                  aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                    >  3:97518/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                 >  1:833184/1‑149 (MQ=255)
                                                                                                                                      aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                        >  4:497040/1‑77 (MQ=38)
                                                                                                                                      aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                        <  3:497040/77‑1 (MQ=37)
                                                                                                                                      aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAt                >  4:474901/1‑149 (MQ=255)
                                                                                                                                       aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTTTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt               >  3:173208/1‑149 (MQ=255)
                                                                                                                                       aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt               >  2:641598/1‑149 (MQ=255)
                                                                                                                                            caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt          <  2:554343/149‑1 (MQ=255)
                                                                                                                                                 cTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctg     >  3:924662/1‑149 (MQ=255)
                                                                                                                                                 cTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctg     >  2:74292/1‑149 (MQ=255)
                                                                                                                                                    aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctgact  <  4:127200/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACT  >  NC_000913/2286167‑2286462

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