Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,668 C→T 6.0% 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 6.0% 644.4 / 23.6 233intergenic (+16/‑417)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (122/97);  new base T (4/10);  total (126/107)
Fisher's exact test for biased strand distribution p-value = 5.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781519‑781817
                                                                                                                                                      |                                                                                                                                                     
tCGAGTCATCCTGCACGACCCACCATCCTGAATGATAAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGccc                                                                                                                                                       >  2:275921/1‑150 (MQ=255)
  gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                                     >  1:5470101/1‑150 (MQ=255)
     taaTCCTGCACGACCCAACATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                  <  2:5072270/148‑1 (MQ=255)
     tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                  >  1:1319748/1‑150 (MQ=255)
        tCCTGCACGACCCACCATCCTGAATGAGTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                               >  2:3016209/1‑150 (MQ=255)
        tCCTGCAAGACCCACCATCCTGATTTATTAAGTCATCATAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                               <  2:3037766/150‑1 (MQ=17)
           tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCtt                                                                                                                                            >  2:3389496/1‑150 (MQ=255)
           tGCAAGACCCACCATCCTTATTGATTTAGGCAGCATAAGCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                            <  2:4496273/150‑1 (MQ=17)
            gCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                           <  1:1197788/150‑1 (MQ=255)
            gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                           >  1:4495136/1‑150 (MQ=255)
             cccGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                          <  1:3117922/148‑1 (MQ=255)
             cACGACCGACCATCCTTTTTTATTTATGCAGCCTAATCCCGGTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                          <  2:5344934/150‑1 (MQ=14)
             cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                                          <  2:1936042/150‑1 (MQ=255)
              aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                         >  1:4341404/1‑150 (MQ=255)
               cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                                        <  2:1786100/150‑1 (MQ=255)
                  cccACCATCCTGAATGATTAAGGCAGCATACCCCCGCACGGGGTCGTTATCTCAGTGGTTAGAGCAGTTCACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTTTTTTCt                                                                                                                                     <  1:1147844/150‑1 (MQ=17)
                  cacaccaTCCTGTTTGATTAATGCATCTTAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCt                                                                                                                                     <  2:2062591/148‑1 (MQ=17)
                   ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTTTTTacat                                                                                                                                    >  2:5242796‑M2/1‑146 (MQ=255)
                    caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  1:4633906/1‑150 (MQ=255)
                    caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  2:4069015/1‑150 (MQ=255)
                      ccaTCCTGAATGATTCAGGCAGCATACTCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                 <  2:3285419/150‑1 (MQ=255)
                        cgcccgAAAGATTAAGGCAGCATACTCCCGCAAGGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:418786/145‑1 (MQ=17)
                        aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:1266848/150‑1 (MQ=255)
                          ccggcaagATTTAGGCCGCCTACGCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                             <  2:3804495/143‑1 (MQ=11)
                           cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATAAGCGAt                                                                                                                            >  2:4834408/1‑150 (MQ=255)
                            tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                           >  2:2633753/1‑150 (MQ=255)
                             gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATAc                                                                                                                          >  2:4570804/1‑150 (MQ=255)
                                tGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCa                                                                                                                       <  2:2265849/150‑1 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:4426116/150‑1 (MQ=255)
                                  aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:782611/150‑1 (MQ=255)
                                        gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                               >  2:1520730/1‑150 (MQ=255)
                                          agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                             >  2:1341104/1‑150 (MQ=255)
                                            caTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACCAAcc                                                                                                           >  2:4113771/1‑150 (MQ=255)
                                             aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                        <  2:3465197/120‑1 (MQ=255)
                                              tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                         <  1:1198174/150‑1 (MQ=255)
                                                 tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      >  2:3960605/1‑150 (MQ=39)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  2:4018495/150‑1 (MQ=33)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:662897/150‑1 (MQ=33)
                                                        aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:4678950/150‑1 (MQ=33)
                                                          gggTCGTTTGTCAAGTTGGTAGAGCTGTTGAATTTTAAACAATTGGCCGACGGTTCGAAGCCTGCACGACCCACCAATTTAAAAAAGCGCCCAAACGCCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  2:4632020/150‑1 (MQ=11)
                                                          gggTCGGTAGCTGAGTTTGTTGTTCTGTTGTCTTTTAATCTATTGGTCGCATGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  1:1413928/150‑1 (MQ=11)
                                                            gTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                           <  1:5209880/150‑1 (MQ=17)
                                                               gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                        >  1:2000925/1‑150 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            <  1:5292019/129‑1 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            >  2:5292019/1‑129 (MQ=32)
                                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                       >  1:4995909/1‑150 (MQ=32)
                                                                 tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                                      <  1:4289101/150‑1 (MQ=32)
                                                                   gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                    <  2:1319748/150‑1 (MQ=32)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   <  1:818955/150‑1 (MQ=32)
                                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   >  1:5063542/1‑150 (MQ=32)
                                                                     tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                  <  1:2633628/150‑1 (MQ=32)
                                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCAAGGATTCGAACCTGCAGCAGGTGTGTGTTGAGCGCAGCGa                                                                               >  2:599863/1‑150 (MQ=255)
                                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                               >  2:3279174/1‑150 (MQ=17)
                                                                          tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                             <  1:4743045/150‑1 (MQ=32)
                                                                            gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                           >  1:3438665/1‑150 (MQ=17)
                                                                             tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                          <  2:4542379/150‑1 (MQ=32)
                                                                                 gaaGTTTAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCTAAACAACGGa                                                                      <  2:1869540/148‑1 (MQ=255)
                                                                                   aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    >  1:5101210/1‑150 (MQ=32)
                                                                                    gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                   >  1:3053346/1‑150 (MQ=32)
                                                                                     ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCAGAAGCAGTTAGGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                  >  1:2219301/1‑150 (MQ=11)
                                                                                      tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCCAACCTGCAGCAGGTTTGAGAGGAGCGCAGAGAAACAAAGGAGCCGc                                                                 >  1:1696659/1‑150 (MQ=11)
                                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                             >  2:1988189/1‑150 (MQ=32)
                                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACGCAAAGATACGATCCTGCACCAGTTTTCAGATGATCGCAGCGAAGCGCACGAGCCGCTcgc                                                             >  1:3612563/1‑150 (MQ=11)
                                                                                               aTCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGAAGTTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:3233045/150‑1 (MQ=11)
                                                                                               aTCAATTGGCCGCAGGTTCGAACCCTGCCCGACCCACCACTGTAAAAACGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:2360506/150‑1 (MQ=11)
                                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGATCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  2:2296445/1‑150 (MQ=17)
                                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  1:404248/1‑150 (MQ=17)
                                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                     >  2:1472973/1‑150 (MQ=32)
                                                                                                    ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   >  2:381594/1‑150 (MQ=17)
                                                                                                    ttGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   <  2:4120873/150‑1 (MQ=17)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTTAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               <  2:1110520/150‑1 (MQ=12)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:419337/1‑150 (MQ=32)
                                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGCAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:4223706/1‑150 (MQ=17)
                                                                                                         cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATGTGCGATACTCAAATATTCGTACCTGCACCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg                                              >  1:4378854/1‑150 (MQ=11)
                                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                        >  1:3682750/1‑150 (MQ=32)
                                                                                                                tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAGCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                       >  2:2846183/1‑150 (MQ=17)
                                                                                                                  gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                     >  1:3017862/1‑150 (MQ=32)
                                                                                                                    aTCCTGCACGACCCACCAATTTAAAAATGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                   <  2:754114/150‑1 (MQ=255)
                                                                                                                       cTGCACGACCAACCAATGTAAAAAAGCGCCCTAATGGCGCTTTTTTGATATATGCGATACTAAAAGATTCTAACCTGCAGCTGGTTTTTGTTGAGCGCAGCTAAACAACGGATCCGATCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:4137630/150‑1 (MQ=11)
                                                                                                                        tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                               >  2:5110497/1‑150 (MQ=32)
                                                                                                                           aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgcaggattcgaacctgcgaccaagtgagtacaagtcaactgctctaccaactgggctaacgacccgacgtggtgggtgatgacgggatcgaac                            <  2:4857347‑M2/150‑109 (MQ=255)
                                                                                                                           aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  2:4287498/150‑1 (MQ=32)
                                                                                                                           aCGACCCACAATTTAAAAAACGCCCCCTACGGGGGTTTTTTTGCTTTATGCTATACTCAAAGATTCCACCCTGCTGCTGTTTTGTGTTGCGGGCTCCCAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  1:4954983/150‑1 (MQ=2)
                                                                                                                            cgccccACCACTGTAAAAAACCGCCCTAAAGGCGCTTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCTGTTTTTAGTTGTGCGCAGCGAAACAACGGAGCCGCTCGCTGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                           <  2:1022029/147‑1 (MQ=11)
                                                                                                                             gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTCTGAAAGATTAGAACCTGCAGCAGGTTTGAGTTGAGCGGAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCGGAAGGGCGAGCGAAGCGAGTCATc                          >  1:1843078/1‑150 (MQ=11)
                                                                                                                              aCCCACCACTGAAAACAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGACACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATcc                         <  2:4445137/150‑1 (MQ=11)
                                                                                                                                  accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAATCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:2005787/1‑150 (MQ=17)
                                                                                                                                  accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:5650/1‑150 (MQ=17)
                                                                                                                                    caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2094423/1‑150 (MQ=32)
                                                                                                                                    caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2535318/1‑150 (MQ=32)
                                                                                                                                     atagcaaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                  >  1:2693030/6‑150 (MQ=255)
                                                                                                                                     aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGAAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTTCACGa                  >  2:3527676/1‑150 (MQ=12)
                                                                                                                                         tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca              >  2:2871559/1‑150 (MQ=32)
                                                                                                                                          aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAACCTGCACGACccac             >  1:4071702/1‑150 (MQ=17)
                                                                                                                                           aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  2:665264/1‑150 (MQ=32)
                                                                                                                                           aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:155835/1‑150 (MQ=32)
                                                                                                                                            aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:4992355/150‑1 (MQ=17)
                                                                                                                                            aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:1634060/150‑1 (MQ=32)
                                                                                                                                             aaaGCGCCCTAATGGCTATTTTTTGCTATCTGCGATACTCAAATATTCTAACCTGATGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa          <  1:4116458/150‑1 (MQ=11)
                                                                                                                                               atgcgcCCTTAACTCGATTTTTTGCTTTCTGCGATCCTCAACGACACGAAACTGCAGCTGGTTTGAGTTGTTGAGAGCGTAACTACGTAGGCTCTCGCTGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         <  2:5456368/148‑1 (MQ=2)
                                                                                                                                                      |                                                                                                                                                     
TCGAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781519‑781817

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: