Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,668 C→T 8.7% intergenic (+16/‑417) 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,6680CT 8.7% 234.8 / 13.7 93intergenic (+16/‑417)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (48/36);  new base T (6/2);  total (54/39)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAA  >  NC_000913/781519‑781816
                                                                                                                                                     |                                                                                                                                                    
tCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGccc                                                                                                                                                      >  2:88631/1‑150 (MQ=255)
tCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGccc                                                                                                                                                      >  2:1693428/1‑150 (MQ=255)
    gTCATCCTGCCCGACCCTCCATCCTGACTGCTTTAGGCAGGATAAGCCCGCCAGGGGTCGTTAGCTGTGGTGGTAGTGCTGTGGTCTTTTATTCTGGGGGTCGCAGGTTGGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                  <  2:1224849/150‑1 (MQ=12)
           tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:1086440/1‑150 (MQ=255)
           tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:1441121/1‑150 (MQ=255)
           tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:866229/1‑150 (MQ=255)
             ctcGACCCACCATCCTTAATGATTAAGTCATCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:1361570/148‑1 (MQ=255)
             cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:763534/150‑1 (MQ=255)
             cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:1361728/150‑1 (MQ=255)
              ccgcccccccATCCTGAATGATTAAGGCAGCATACGCCCGCATGGGGTCGTTAGCTCCGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCTGGTGCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                        <  1:402252/146‑1 (MQ=17)
               cGACCCACCATCCTGTATGATTAAGTCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                       <  1:2136437/150‑1 (MQ=255)
                gACCCACCATCCTGAATTATTAATGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                      <  1:1908776/150‑1 (MQ=255)
                gACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                      <  1:1666769/150‑2 (MQ=255)
                    caccaTCCTGACTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:379211/150‑1 (MQ=255)
                    caccaTCCTGAATGATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:1356015/150‑1 (MQ=255)
                    caccaTCCTGAATGACTACGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:1356001/150‑1 (MQ=255)
                            tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCATTTGGTAGAGCAGTTGACTTTTAAAAAATTGGTCGCAGGTTAGAATCCTTCACGACCCACCAAAGTAAAAAATCGCCCTAAATGCCCATTTTTGCTATCTGCGATa                                                                                                                          >  1:1374221/1‑150 (MQ=17)
                             gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAAc                                                                                                                         >  2:2269421/1‑150 (MQ=255)
                                 gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTTCGAAGCTCaa                                                                                                                     >  2:1442877/1‑150 (MQ=255)
                                 gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCaa                                                                                                                     >  2:2229477/1‑150 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCGATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCcaa                                                                                                                    >  2:1973638/1‑150 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    >  2:707560/1‑150 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    <  1:1599681/150‑1 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    >  2:2454536/1‑150 (MQ=255)
                                   ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGAGACTCAAAg                                                                                                                   >  1:1998697/1‑150 (MQ=255)
                                   ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAg                                                                                                                   >  1:1998686/1‑150 (MQ=255)
                                      aGGCAGCCTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:49970/150‑1 (MQ=255)
                                      aGGCAGCATAATCCCGCCAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:191127/150‑1 (MQ=255)
                                      aGGCACCATAATCCCGCTTGGGGTCGTTCTCTCATTTGGTAGAGCATTCGACTTTAAATCCTTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:488027/150‑1 (MQ=14)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                            >  2:75919/1‑150 (MQ=255)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                            <  2:591721/150‑1 (MQ=255)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCGaa                                                                                                            >  2:75787/1‑150 (MQ=255)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGGAATCAAGGGGTCGCCGGGTCGACTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCGGGTGTGCTAGCGGCGATACGCCAAGaggcgac                                                                                                            >  2:292828/1‑144 (MQ=14)
                                                 taccggtaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                     <  2:1951831/143‑1 (MQ=18)
                                                 tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                     <  1:1693428/150‑1 (MQ=39)
                                                   ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgaagca                                                                                                   >  2:1668521/1‑150 (MQ=38)
                                                   ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCACCGATACTAAAATATGCGAAACTgcacca                                                                                                   >  2:1107295/1‑150 (MQ=11)
                                                    cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATTGTAGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAATCGCCCGAAAGGCGCTTTTTTTCTATCTGCGATACTCAAAGATTCGAACcagcagcag                                                                                                  >  1:643601/1‑150 (MQ=11)
                                                      actcgggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTGCTAAAAGCGATACTCAAAGATTCGAACCTGCAGCAggg                                                                                                >  2:1755494/5‑149 (MQ=11)
                                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                            <  1:1755494/150‑1 (MQ=32)
                                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                            >  2:107379/1‑150 (MQ=32)
                                                           ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:787930/150‑1 (MQ=32)
                                                           ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:12895/150‑1 (MQ=32)
                                                           ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:1866242/150‑1 (MQ=32)
                                                           ggTCGTTAGCTCAGTTGGTACAGCAGTTGACTTTTAATCAATTGGTCGCAGGTACGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:1906065/150‑1 (MQ=17)
                                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCAAAGATTCGAAACTGCAGCAGGTttgagtt                                                                                         >  1:765114/1‑150 (MQ=12)
                                                               gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:1006442/150‑1 (MQ=17)
                                                               gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:264588/150‑1 (MQ=17)
                                                               gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:264585/150‑1 (MQ=17)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                  >  1:1743595/1‑150 (MQ=32)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                  >  1:2432370/1‑150 (MQ=32)
                                                                       aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                               <  1:360121/150‑1 (MQ=32)
                                                                        gTTGGTAGCGCAGTTGACTTTTCATCACTTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAACAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                              <  2:1191550/150‑1 (MQ=11)
                                                                              agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGATGAGCGCAGCGAaacacc                                                                        >  1:1263289/1‑148 (MQ=17)
                                                                               gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTCTCAGCGATACTCAAAGATTCGAACCTGCCGCAGGTGTGAGTTGAGCGCAGCGAAACAACg                                                                       >  2:2425811/1‑150 (MQ=255)
                                                                                agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                      >  2:410122/1‑150 (MQ=32)
                                                                                agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                      >  2:281553/1‑150 (MQ=17)
                                                                                    gTTGTCTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAATTTTCGAACCTGCAGCAGGTTTTAGTTGAGCGCAGCTAAACAACGGAGcc                                                                  <  2:1553320/150‑1 (MQ=11)
                                                                                    gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTAACCACCAATGTAAAAAAGCGACCTAAAGTCGCTTTTTTGCTATATGCTATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGATCTCAGCGATACAACGGAgac                                                                  >  1:972322/1‑148 (MQ=11)
                                                                                       gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                               >  1:1328716/1‑150 (MQ=32)
                                                                                           tttAATCATTTGGTCGCAGGTTCTAATCCTGCACGACCCACCAATGTAAAAATGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAATTATTCGAACCTGCAGCAGGTTTTTGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                           <  2:885541/150‑1 (MQ=11)
                                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACGACGGAGCCGCTCgcggcgc                                                       >  2:1442594/1‑149 (MQ=17)
                                                                                                  acttGGTCGCAGGTTCGACTCCTGCACGACCCCCCACTGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1743983/148‑1 (MQ=255)
                                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1619191/150‑1 (MQ=17)
                                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1619220/150‑1 (MQ=32)
                                                                                                          gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                            <  2:107018/150‑1 (MQ=12)
                                                                                                          gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                            <  2:107040/150‑1 (MQ=17)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCGAAGATTCGACCCTGCCGGAGCGGTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGATGGgcgag                                       >  1:854546/1‑150 (MQ=11)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGgcgag                                       >  2:328171/1‑150 (MQ=17)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                       >  2:505839/1‑150 (MQ=32)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                       >  1:864806/1‑150 (MQ=32)
                                                                                                                  gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  2:1521121/1‑150 (MQ=32)
                                                                                                                  gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  1:1498927/1‑150 (MQ=17)
                                                                                                                  gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  1:1049051/1‑150 (MQ=17)
                                                                                                                   taTCCTGCACGACCCACCACTGTAAACATGCGCCGTAAAGGCGCTTTTTTGCTATCTGCGTTACTCAACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGgcgagcgac                                   <  2:2534887/149‑2 (MQ=11)
                                                                                                                   aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCGTTTTTGCTATCGGCGATACTCAAAGATTAGAACCTGCAGCAGGTTTGAGTTGTGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGagcgaa                                   >  1:845174/1‑150 (MQ=11)
                                                                                                                     tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAACAGCAATATGAGAAGCTGCAGCAGGTTAGAGGGTAGCGCATCGAACCGACGGAGCCGCTCGCGCCGACGGCCCGGGGGGCGagagacgc                                 >  2:2076767/1‑150 (MQ=2)
                                                                                                                      ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                >  2:1884215/1‑150 (MQ=32)
                                                                                                                      ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATACTCAAAGATGCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACCGGAGCCGCTCGCGGCGACGGCCCGGGGGGCGagcgaagcg                                >  2:1884183/1‑150 (MQ=11)
                                                                                                                                  accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                    <  2:643601/150‑1 (MQ=32)
                                                                                                                                       tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc               >  1:638243/1‑150 (MQ=32)
                                                                                                                                         tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCATCCTGCACGACcca             >  1:1790327/1‑150 (MQ=17)
                                                                                                                                         tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca             >  1:320791/1‑150 (MQ=32)
                                                                                                                                         caaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca             >  2:1427141/2‑150 (MQ=12)
                                                                                                                                          aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:2139923/1‑150 (MQ=12)
                                                                                                                                          aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:1690145/1‑150 (MQ=12)
                                                                                                                                          aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:1689999/1‑150 (MQ=12)
                                                                                                                                          aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTAGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:303745/1‑150 (MQ=255)
                                                                                                                                          aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  1:1831401/1‑150 (MQ=32)
                                                                                                                                           aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc           >  1:2307365/1‑150 (MQ=12)
                                                                                                                                            aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca          >  1:1468114/1‑150 (MQ=32)
                                                                                                                                            aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca          >  1:847538/1‑150 (MQ=32)
                                                                                                                                             aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa         >  2:1013508/1‑150 (MQ=17)
                                                                                                                                                  gcCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa    <  1:1441121/150‑1 (MQ=17)
                                                                                                                                                  gcCCTAAAGGCGATTTTTTGATATTTGCTATACTAAAATTTTCTAACCTGCTGCTGTTTTGTTTTGCGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa    <  1:1086440/150‑1 (MQ=11)
                                                                                                                                                    ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGCCCCACCAATGTaaaa  >  1:1653709/1‑150 (MQ=17)
                                                                                                                                                    ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa  >  1:1423433/1‑150 (MQ=32)
                                                                                                                                                    ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa  >  2:264746/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                    
TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAA  >  NC_000913/781519‑781816

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