Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,719 A→T 6.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 6.0% 840.3 / 5.3 361intergenic (+67/‑366)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (195/123);  new base T (5/26);  total (202/159)
Fisher's exact test for biased strand distribution p-value = 1.38e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTT  >  NC_000913/781570‑781867
                                                                                                                                                     |                                                                                                                                                    
 cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAAGCGCCCTAAAGGCGCTTTTTTTCTATATGCGATACTAAAATATTAGAACcagcagcag                                                                                                                                                     >  1:8264046/1‑150 (MQ=11)
    caGGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAAGTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                                                                                  <  1:5680287/149‑1 (MQ=18)
    aaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                                                                                  <  1:558456/150‑1 (MQ=34)
     aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                                                                 <  1:7323938/150‑1 (MQ=33)
       gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                                                                               <  1:2566796/150‑1 (MQ=32)
         gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTCCTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAGACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                                                                             <  2:8776624/150‑1 (MQ=12)
         gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                                                                             <  2:717327/150‑1 (MQ=32)
          tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                                                                            <  1:2998609/150‑1 (MQ=32)
          tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                                                                            <  2:4539920/150‑1 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                         <  1:5094542/150‑1 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                         >  2:6914316/1‑150 (MQ=32)
               aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                                                                                                                       >  2:6085124/1‑150 (MQ=32)
                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                     <  1:2771717/150‑1 (MQ=32)
                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                     >  1:9279951/1‑150 (MQ=32)
                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                     <  1:7794658/150‑1 (MQ=32)
                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                     <  2:3216956/150‑1 (MQ=32)
                   cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                                                                   <  1:4130841/150‑1 (MQ=32)
                      ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                                                                                <  1:5333976/150‑1 (MQ=32)
                       ttgTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                                                               <  2:6811329/148‑1 (MQ=17)
                       tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                                                               <  2:3035323/150‑1 (MQ=32)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGAGCGCCGCGAaac                                                                                                                              >  2:6389307/1‑150 (MQ=17)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                              >  2:3598941/1‑150 (MQ=32)
                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                                                                             >  2:8526686/1‑150 (MQ=32)
                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                                                                             >  2:1912963/1‑150 (MQ=17)
                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTTAGTTGAGCGCAGCGAaaca                                                                                                                             >  1:6676131/1‑150 (MQ=17)
                              gctgTTGACTTTTAATAAATTTGTCGCAGGTTCGACTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                                                                        <  2:6580225/147‑1 (MQ=11)
                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTGTAAAAAAGCGCCCTCAAGGCGCTTTGTGGCTCTCTGCGATCCTCGGAGATTCGAACCGGCAGCAGGTTGGAGTTGTGCGCAGCGAAACACCGGAGCCg                                                                                                                    >  2:7193618/1‑150 (MQ=11)
                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                                                                    >  2:1875839/1‑150 (MQ=32)
                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAACCTGCAGCAGGTGTGATTTGAGCGCAGCGAAACAACGGAGACg                                                                                                                    >  1:9125305/1‑150 (MQ=255)
                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACACACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGCGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                                                                  >  2:5779339/1‑150 (MQ=17)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                                                               >  1:4636962/1‑150 (MQ=32)
                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                                                                               >  1:9434276/1‑150 (MQ=32)
                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTTCGATACTCGTAGTTTCGAACCTGCAGCTGGTTTGAGTTGAGCGGAGCGAAGCAACGGAGCCGCTcgcg                                                                                                              >  2:276013/1‑150 (MQ=11)
                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                                                                              <  2:5948784/150‑1 (MQ=32)
                                         ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGTTTCGAACCTGAAGCTGGTTTGAGTTTCGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                             >  1:4955084/1‑150 (MQ=11)
                                            ctCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                                                                          <  2:4705875/149‑1 (MQ=17)
                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                                                                         >  1:6301358/1‑150 (MQ=32)
                                             tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGAAACACAAAGATACGAAACTGCACCAGGTTTGAGGTTAACGCAGCGAACCAACCGAGACTCTCTCGCCGAc                                                                                                         >  2:1609665/1‑150 (MQ=11)
                                                aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                                                                      >  2:9190091/1‑150 (MQ=32)
                                                     tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                                                                 >  2:5318438/1‑150 (MQ=32)
                                                            ttCTAATCCTGCCCGACACACCAATGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCTAACCTGCAGAAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                                                                          <  2:9306612/150‑1 (MQ=11)
                                                               gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGCTTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                       >  1:3292813/1‑150 (MQ=17)
                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                      >  1:1514058/1‑150 (MQ=32)
                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                      >  1:2307750/1‑150 (MQ=32)
                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTGGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                      >  1:6838339/1‑150 (MQ=17)
                                                                 aTCCCGCACGACCCACCTATTTAAAAAAGCGCCCTAACGGCGCTTTTTTGATATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAATCAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                                                                     <  2:8951507/150‑1 (MQ=11)
                                                                    cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAAGCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAAAGGAGCCGCTCGCGGCGACGGCCCGTAGGGCGagagaagcga                                                                                  >  2:3288637/1‑150 (MQ=11)
                                                                     tGCACGTCCCACCAATGTAAAAAAGCGCCCTATTGGCGCTTTTTTGATATCTGCGTTACTCAAATATTCTAACCTGCATCAGGTTTGAGTTGAGCTCAGCGAAACAACGGAGCCTCTCGCGGCGTCTGACCGAAGGGCGAGCGAAGCGAg                                                                                 <  2:1561793/150‑1 (MQ=11)
                                                                     tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                 >  2:3785606/1‑150 (MQ=32)
                                                                     tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                 <  2:5595766/150‑1 (MQ=32)
                                                                          gACCCACAAATTTAAAAACGCGCCCTAACGGCTTTTTTTTGCTTTCTGCTATACTCAAATTTTCGACCCGGCAGCTGGTTTGCTTTGAGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                            <  1:8749703/150‑1 (MQ=11)
                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                       >  1:5039778/1‑150 (MQ=32)
                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGGCATCCTGCa                                                                       >  2:1020967/1‑150 (MQ=17)
                                                                                ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                      >  1:1032628/1‑150 (MQ=32)
                                                                                   aTGAAAAAACCCCCCCAAACGGCGCTTTTTTGCTAACTGCGATACTAAAAGATCCCACCCTGCTGCTGGTTTGCGTGGCGCGCCGCCAACCAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc                                                                   <  1:3885787/150‑1 (MQ=2)
                                                                                    tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                  >  1:1241212/1‑150 (MQ=32)
                                                                                       aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                               >  1:606911/1‑150 (MQ=12)
                                                                                       aaaaaaGCGCCCTAACGGCGATTTTTTGCTATCTGCGATCCTAAAAGATTCTAACGTGCAGCTGTTTTGTGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                               <  1:8668550/150‑1 (MQ=11)
                                                                                       aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                               >  1:6447751/1‑150 (MQ=32)
                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                             >  2:2224821/1‑150 (MQ=32)
                                                                                           aaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                           >  2:9146392/1‑150 (MQ=12)
                                                                                           aaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                           >  2:346029/1‑150 (MQ=12)
                                                                                                   tAAAGGCGTTTTTTTGCTATATGCGATACTCAAAGATCCGAACCTGCCGCAGGTTTGCGTTTAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa                                                   <  2:2483314/150‑1 (MQ=11)
                                                                                                   tAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa                                                   >  2:4389366/1‑150 (MQ=32)
                                                                                                       ggCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcgc                                               <  2:9683019/150‑1 (MQ=32)
                                                                                                       ggCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcgc                                               <  2:4894977/150‑1 (MQ=32)
                                                                                                        gcgcTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGcc                                              >  2:5061312/1‑150 (MQ=32)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                            >  1:1978990/1‑132 (MQ=32)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                            <  2:1978990/132‑1 (MQ=32)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCtca                                          >  1:7572449/1‑148 (MQ=17)
                                                                                                            ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                          >  1:3480333/1‑150 (MQ=32)
                                                                                                            ttttttGCTATCTGATATACTCAAAGTTTCGAACCTGCAGGAGGTTTGATTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                          <  2:2184938/150‑1 (MQ=11)
                                                                                                                ttGCTATCTTCGATACTAAAATATTCGAACCTTCAGGTGGTTTGTGTTGTGCGCAGGGAAACAACGGTGCCGCTCGCGGCGACGGCCCGTAGGGCGAGCGATGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                      <  1:7933915/150‑1 (MQ=11)
                                                                                                                ttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                      <  2:4680949/150‑1 (MQ=32)
                                                                                                                  gCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                    >  2:7132273/1‑150 (MQ=11)
                                                                                                                              aCTCACAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAtt                        <  1:9644553/150‑1 (MQ=11)
                                                                                                                              aCTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAtt                        >  1:3678215/1‑150 (MQ=11)
                                                                                                                              aCTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCACTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAtt                        >  2:4182095/1‑150 (MQ=11)
                                                                                                                                tCAAAGATTCTAACCTGCAGCTGGTTTGAGGTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      <  1:8554914/150‑1 (MQ=12)
                                                                                                                                tCAAAGATTCTAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCTCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      <  2:3855428/150‑1 (MQ=17)
                                                                                                                                tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      >  2:8230363/1‑150 (MQ=32)
                                                                                                                                tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      >  1:2148266/1‑150 (MQ=32)
                                                                                                                                tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      >  2:8472452/1‑150 (MQ=32)
                                                                                                                                tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                      >  2:9496367/1‑150 (MQ=32)
                                                                                                                                 cAAAGATTCGAACCTGCAGCAGGTTTGTGTTGAGCGCCGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGCAGTGCGAGCGCAGCGTGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                     <  2:7974963/150‑1 (MQ=11)
                                                                                                                                 cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                     <  1:595106/150‑1 (MQ=17)
                                                                                                                                  aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATTCAgg                    >  1:6344260/1‑150 (MQ=18)
                                                                                                                                  aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                    >  1:810248/1‑150 (MQ=33)
                                                                                                                                  aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                    >  2:1310821/1‑150 (MQ=33)
                                                                                                                                   cagtttCGAACCTGCACCCGTTTTGAGTGGAGGGCAGGGAACCAACGGAGCCTCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCTAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGc                   <  2:1715834/146‑1 (MQ=11)
                                                                                                                                          gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa            >  2:2863358/1‑150 (MQ=40)
                                                                                                                                          gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa            >  2:6140324/1‑150 (MQ=40)
                                                                                                                                          gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa            >  2:7407680/1‑150 (MQ=255)
                                                                                                                                           aacaatcagcagTTTTTAGTTGACGTCAGCTAAACACCGTAGCCTCTGGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa           <  1:4803978/144‑1 (MQ=14)
                                                                                                                                             ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAtt         >  1:4529668/1‑150 (MQ=255)
                                                                                                                                                gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa      <  1:716947/150‑1 (MQ=255)
                                                                                                                                                   gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt   >  2:7869004/1‑150 (MQ=255)
                                                                                                                                                   gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt   >  2:7602936/1‑150 (MQ=255)
                                                                                                                                                   gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt   >  2:8981008/1‑150 (MQ=255)
                                                                                                                                                   gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt   >  1:6923490/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTT  >  NC_000913/781570‑781867

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: