Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,719 A→T 8.0% intergenic (+67/‑366) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,7190AT 8.0% 535.2 / 13.0 237intergenic (+67/‑366)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (131/74);  new base T (2/24);  total (133/104)
Fisher's exact test for biased strand distribution p-value = 2.58e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  NC_000913/781570‑781868
                                                                                                                                                     |                                                                                                                                                     
 cgctatgggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                                                                                      <  1:6005516/144‑1 (MQ=17)
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                                                                     <  2:1331645/150‑1 (MQ=35)
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                                                                     <  1:2377540/150‑1 (MQ=35)
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                                                                     >  2:4781460/1‑150 (MQ=35)
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAATTTTAATAAATTGGTCGCAGGTTCGAATCCATCACCACCCACCAAATTAAAAAAGCGCCCAAAAGGCGCATTTTTGCTATCAGCGATAAACAAAGATTCGAACCAGCACCAgg                                                                                                                                                     >  1:640104/1‑150 (MQ=11)
     aGGGGTCGTTAGCACAGTTGGTAGAGCAGCTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                                                                  <  2:698629/150‑1 (MQ=18)
        ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                                                                               <  1:268214/150‑1 (MQ=32)
            gtgaGCTCAGTTGGGAGAGCAGTGGACCTTTAATCAAGAGGTCGCAGGTTCGAATCCGGCACGACCCACCCATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                                                                           <  2:3064543/147‑1 (MQ=11)
            gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                                                                           <  2:3437349/150‑1 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                          >  1:2716436/1‑150 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                          <  1:1389153/150‑1 (MQ=32)
                gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAGCCACCAATGTAATAAAGCGCCCTATAGGCGCTTTTTTGCTATGTGGGATTCTGAATGATTCGTTGCTGCTGCTGGTttgagtttggggg                                                                                                                                       >  2:5605138/1‑144 (MQ=11)
                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                                                                     >  2:3352404/1‑150 (MQ=32)
                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                                                                     <  2:2276768/150‑1 (MQ=32)
                  tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACGCAAAGATTCGAACCTGCAGCAGGTTTGAGATGATCGCAg                                                                                                                                     >  2:7191596/1‑150 (MQ=255)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                               >  1:6110850/1‑150 (MQ=32)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                               >  1:6921823/1‑150 (MQ=32)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                               >  2:3804518/1‑150 (MQ=17)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATAAGCGAAACAAAAAGATACGAACCTGAAGCAGGTTTGAGTGGAGCGCAGCGAaac                                                                                                                               >  1:7182688/1‑150 (MQ=11)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCCGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGATGAGCGCAGCGAaac                                                                                                                               >  2:5999076/1‑150 (MQ=17)
                              gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                                                                         <  2:5527074/150‑1 (MQ=32)
                                aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGTTACTCAAATATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                                                                       <  2:5133759/150‑1 (MQ=255)
                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                                                                   >  1:4448162/1‑150 (MQ=32)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATACGAACCTTCAGCAGGTTTGATTTGAGCGCAGCGAAACACCGGATCCGCTcgc                                                                                                                >  1:6135590/1‑150 (MQ=11)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                      <  2:1862613/128‑1 (MQ=32)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                                                                >  2:2874946/1‑150 (MQ=32)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                                                                >  1:1721223/1‑150 (MQ=32)
                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGCAACTAAAATATTAGAACCTTCAGCAGGTTTGAGTTGAGCGCAGCGAAACAAAGGAGCCGCTctcg                                                                                                               >  1:6772047/1‑150 (MQ=11)
                                            aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACACAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                                                                           >  1:1158256/1‑150 (MQ=17)
                                            aTCAATTGGTCGATGTTTCGAACCCTGCACGACCCACCAATGTAAAAACGCGCCCTACCGGCGATTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGCGCGCAGCTAAACAACGGAGCCGCTCGCGGCGa                                                                                                           <  2:3104588/150‑1 (MQ=11)
                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGCGCCGCTCgcggcgtc                                                                                                          >  2:4598574/1‑148 (MQ=12)
                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                                                                          >  1:3270768/1‑150 (MQ=32)
                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGAGCCCTAAAGGCGCTATTTTGCTATCTGCGAAACTCAAAGATAAGAACCTGCAGCAGGATTGAGTTGAGCGCAGCGAAACAAAGGAGCCGCTCGCGGCGAc                                                                                                          >  1:4016635/1‑150 (MQ=11)
                                                    gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa                                                                                                   >  1:2816994/1‑150 (MQ=32)
                                                    gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAAAGGAGCCGCTCGCGGCGACGGCCCGa                                                                                                   >  1:1885835/1‑150 (MQ=17)
                                                        ctggTTCGATTCCTGCACGACCCACCATTTTAAAAATGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCAAAGATTCGAACCTGAAGCAGGTTTGAGTTGAGCTCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               <  2:6825142/148‑1 (MQ=11)
                                                        cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               <  2:7024089/150‑1 (MQ=32)
                                                          ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                                                                             >  2:4152234/1‑150 (MQ=32)
                                                           gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgaga                                                                                            >  2:3209961/1‑150 (MQ=17)
                                                           gTTCGAATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAATGGCGATTTTTTGCTATCTGGTATACTCAAAGATTCGAACCTGCAGCCGGTTTGTGTTGAGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                                                                            <  2:6704675/150‑1 (MQ=11)
                                                                  tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                                                                     >  2:508619/1‑150 (MQ=32)
                                                                  tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                                                                     >  1:665316/1‑150 (MQ=32)
                                                                    cTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                                                                   <  2:1438643/150‑1 (MQ=17)
                                                                     tGCAGGCCCCACCTCTGTAACAACGCGCCCTTCAGGCGATTTTTGGCTTTCTGGGATACTCCAAGAGCCTCCCCTGCATCCGTTTTGAGTTGTGCGGAGCGAAACACCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                  <  2:6877847/150‑1 (MQ=2)
                                                                      gCACGACCCCCCCCTGTACAAACGCGCCCTCCCGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGt                                                                                 <  2:5756067/150‑1 (MQ=11)
                                                                         cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                                                                              <  2:5895106/150‑1 (MQ=32)
                                                                          gACCCACCCCGGTAAAAACGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTTAAGCTGGTTTGAGTTGAGCGCAGCGAAACCACGGAGCCGCTCTCGGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                             <  2:4700708/150‑1 (MQ=11)
                                                                          gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                             >  1:3482590/1‑150 (MQ=32)
                                                                          gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                             >  1:3181060/1‑150 (MQ=17)
                                                                             ccaccaATGTAAAAACGCGCCCTACCGGCGATTTTTTGCTATCTGCGATACTAAATTATTCCAACGTGCTGCTGGTTTTTTTTGCGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                                                                          <  1:5281128/150‑1 (MQ=11)
                                                                               accaATTAAAAAAAGCGCCCTAATGGCTCTTTTTTGATATATGCTATACTCAAAGATTCTAACTTGCTGCAGGTTTGTGTTGTGCGCAGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                        <  1:4020203/150‑1 (MQ=11)
                                                                                  atagcaaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                     >  2:1758494/6‑150 (MQ=255)
                                                                                     gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                                                                  >  1:1824005/1‑150 (MQ=32)
                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                                                                 >  2:3265432/1‑150 (MQ=32)
                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                                                                 >  2:2625590/1‑150 (MQ=32)
                                                                                       aaaaatGCGCCCTATCGGGGTTTTTTTGTTATCTGCGATACTAAAAGATTCGAACCTGCTGCTGTTTTGTGTTGAGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                <  1:6854306/150‑1 (MQ=11)
                                                                                       aaaaacGCGCCCTAATGGCGTTTTTTTGATATCTGCTATACAAAAATATTCGAACCTGCAGCTGTTTTGTGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                <  1:5241761/150‑1 (MQ=11)
                                                                                       aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                >  2:4232602/1‑150 (MQ=12)
                                                                                       aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                >  1:4548901/1‑150 (MQ=12)
                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                              >  2:449079/1‑150 (MQ=32)
                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTGCTAACTGCGATACTAAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                              >  1:5494367/1‑150 (MQ=17)
                                                                                          aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                             >  2:1038283/1‑150 (MQ=32)
                                                                                           acgcgcCATAATGGTGTTTTTTTGATATCTTCTATAATCAATTATTCCAACCTGCTGCTGGTTTGTGTTGCGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            <  1:1446335/148‑1 (MQ=11)
                                                                                           aaGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGTGCGCAGCGAACCAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            <  2:7068120/150‑1 (MQ=255)
                                                                                             gcgcCCTAACGGCGTTTTTTTGCTATCTGCTATACACAAAGATTCGAACCTTCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGt                                                          <  1:4640279/150‑1 (MQ=11)
                                                                                                 ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  1:4219503/1‑150 (MQ=32)
                                                                                                 ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  2:1557234/1‑150 (MQ=32)
                                                                                                 ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  1:4219895/1‑150 (MQ=32)
                                                                                                 ccTAAAGGCGCTTTTTTGCTATCTGCGATACACAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACCGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  1:1913240/1‑150 (MQ=17)
                                                                                                  cTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa                                                     >  1:6045381/1‑150 (MQ=17)
                                                                                                      aGGCGCTTTTTTGCTATCTGCGATACTCAACGTTTCTAACCTGCTGCAGGTTTGAGTTGAGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcg                                                 <  2:2151330/150‑1 (MQ=11)
                                                                                                       ggCGATTTTTTGCTTTCTGGGATACTCAAAGATTCGAACCTGCAGCTGTTTTGAGTTGAGCGCAGCGAAACTACGTAGCCGGTCGCTGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcgc                                                <  2:5065322/150‑1 (MQ=11)
                                                                                                        gcgcTTTTTTGCTATCTGGTATTCTCAACTATTCGAACCTGCAGCCGGTTTGAGTTGAGCGCCGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGcc                                               <  2:6710531/150‑1 (MQ=11)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                           >  2:4757811/1‑150 (MQ=32)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                           >  1:2211433/1‑150 (MQ=32)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                           >  2:2560718/1‑150 (MQ=32)
                                                                                                               tttGATATCTGCTATACTAAAAGATTCGAACCTGCTGCTGTTTTTTTTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                        <  1:5833836/150‑1 (MQ=11)
                                                                                                                ttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                       >  1:4451914/1‑150 (MQ=32)
                                                                                                                                 cAAAGATGCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACCCCGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                      <  2:61962/150‑1 (MQ=11)
                                                                                                                                  aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                     >  1:1562955/1‑150 (MQ=33)
                                                                                                                                  aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                     >  1:776781/1‑150 (MQ=33)
                                                                                                                                     gATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCAt                  >  2:6527947/1‑150 (MQ=35)
                                                                                                                                       ttCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATcc                >  1:4562309/1‑150 (MQ=38)
                                                                                                                                       ctCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACCCCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATcc                <  1:884175/149‑1 (MQ=255)
                                                                                                                                         tgAACATGGAGCATGTATCAGTCGAGCGCAGCGTACCAACGGAGCTGGTGGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              <  2:3682818/149‑1 (MQ=14)
                                                                                                                                         cGAACCTGCAGCAGGTTTGCGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              <  1:6494447/150‑1 (MQ=255)
                                                                                                                                         cGAACCTGCAGCAGGTTTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              <  1:3906266/150‑1 (MQ=255)
                                                                                                                                         cGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              <  1:4858744/150‑1 (MQ=39)
                                                                                                                                          gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa             >  1:1082082/1‑150 (MQ=40)
                                                                                                                                           atcgtGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa            <  2:3963381/146‑1 (MQ=255)
                                                                                                                                           aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa            >  1:5070391/1‑150 (MQ=40)
                                                                                                                                            aCCTGCAGCAGGTTTGAGTTGAGCGCAGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAt           <  1:5619978/150‑1 (MQ=255)
                                                                                                                                             ccttcttcagTTTTGATTTTAGCGCAGCTTATCAACGTAGACGCTATCGGCGACGGACCTTTGGTCGAGCTAAGGGTGTCATCATTCACGACCCACCTATGTAAAAATGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAtt          <  2:1607540/143‑1 (MQ=9)
                                                                                                                                             ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAtt          >  2:812532/1‑150 (MQ=255)
                                                                                                                                              cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAttt         >  2:1131781/1‑150 (MQ=255)
                                                                                                                                               tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTc        <  1:2298795/150‑1 (MQ=255)
                                                                                                                                                 cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCACTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc      >  2:1608797/1‑150 (MQ=255)
                                                                                                                                                   gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt    >  2:6409355/1‑150 (MQ=255)
                                                                                                                                                    cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCTCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt   <  2:3510692/150‑1 (MQ=255)
                                                                                                                                                     agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg  <  2:1529602/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  NC_000913/781570‑781868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: