Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 6.6% intergenic (+20/‑413) 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,6720AG 6.6% 600.0 / 27.3 233intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (122/86);  new base G (4/11);  total (126/107)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 235 total aligned reads displayed.

GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781523‑781817
                                                                                                                                                     |                                                                                                                                                 
gTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                   >  1:5451685/1‑150 (MQ=255)
 taaTCCTGCACGACCCAACATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                  <  2:5072270/148‑1 (MQ=255)
 tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                  >  1:1319748/1‑150 (MQ=255)
       tGCTCGACACACCAGCCTGAATTATTAAGTCAGCCTAAACCCGGTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                            <  2:4496286/150‑1 (MQ=14)
       tGCACGACCCACCATCCTGAATGATTAAGTCAGCATAAGCCCGCTAGGGGTCTTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCtt                                                                                                                                            <  1:614475/150‑1 (MQ=18)
       tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCtt                                                                                                                                            >  2:3389496/1‑150 (MQ=255)
        gCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                           <  1:1197788/150‑1 (MQ=255)
        gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                           <  2:3566775/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:4341397/1‑150 (MQ=255)
           cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                                        <  2:1786100/150‑1 (MQ=255)
              cccACCATCCTGAATGATTAAGGCAGCATACCCCCGCACGGGGTCGTTATCTCAGTGGTTAGAGCAGTTCACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTTTTTTCt                                                                                                                                     <  1:1147844/150‑1 (MQ=17)
               ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTTTTTacat                                                                                                                                    >  2:5242796‑M2/1‑146 (MQ=255)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  2:4069015/1‑150 (MQ=255)
                 accaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTacattg                                                                                                                                  >  2:2842051‑M2/1‑144 (MQ=255)
                  ccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                 >  2:3306028/1‑150 (MQ=255)
                    cgcccgAAAGATTAAGGCAGCATACTCCCGCAAGGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:418786/145‑1 (MQ=17)
                    aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:1266848/150‑1 (MQ=255)
                      ccggcaagATTTAGGCCGCCTACGCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                             <  2:3804495/143‑1 (MQ=11)
                      ccTTAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                             <  1:929359/150‑1 (MQ=255)
                       cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATAAGCGAt                                                                                                                            >  2:4834408/1‑150 (MQ=255)
                       cTGAAGGATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTcattggtgggtc                                                                                                                             >  2:1910387‑M2/1‑138 (MQ=255)
                        tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATa                                                                                                                           >  2:2073183/1‑150 (MQ=255)
                         gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATAc                                                                                                                          >  2:4570804/1‑150 (MQ=255)
                          aaTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTATTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGGTTTTTTGGTATGTGCGATACt                                                                                                                         >  2:1035747/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)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGCTACTCAAAGATTCGaa                                                                                                             >  1:5089308/1‑150 (MQ=255)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACAAAAAGATTAGaa                                                                                                             >  2:1360882/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)
                                                    aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  2:4018495/150‑1 (MQ=33)
                                                    aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:4678950/150‑1 (MQ=33)
                                                    aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:662897/150‑1 (MQ=33)
                                                      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)
                                                              aGCTCAGTTGGTAGAGCATTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                                                                     <  2:473240/150‑1 (MQ=17)
                                                                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)
                                                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                              <  2:434794/150‑1 (MQ=32)
                                                                      tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                             <  1:4743045/150‑1 (MQ=32)
                                                                        gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGGaaca                                                                           >  1:2794146/1‑150 (MQ=12)
                                                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                          <  2:4542379/150‑1 (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)
                                                                                       tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                            >  1:4888794/1‑150 (MQ=17)
                                                                                           aTCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGAAGTTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:3233045/150‑1 (MQ=11)
                                                                                           aTCAATTGGCCGCAGGTTCGAACCCTGCCCGACCCACCACTGTAAAAACGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:2360506/150‑1 (MQ=11)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCAGGTTTGAGATTAGCGCAGCGAAACAACGGACCCGCTCGCGCCGAc                                                       >  2:2961926/1‑150 (MQ=11)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  2:2296447/1‑150 (MQ=32)
                                                                                            tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  1:404248/1‑150 (MQ=17)
                                                                                             cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAACAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg                                                      <  2:3136326/150‑1 (MQ=17)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                     <  2:3989127/150‑1 (MQ=32)
                                                                                              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)
                                                                                                       cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGGTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                            >  2:1323409/1‑150 (MQ=17)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                        >  1:3682750/1‑150 (MQ=32)
                                                                                                                aTCCTGCACGACCCACCAATTTAAAAATGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                   <  2:754114/150‑1 (MQ=255)
                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:3032598/150‑1 (MQ=32)
                                                                                                                   cTGCACGACCAACCAATGTAAAAAAGCGCCCTAATGGCGCTTTTTTGATATATGCGATACTAAAAGATTCTAACCTGCAGCTGGTTTTTGTTGAGCGCAGCTAAACAACGGATCCGATCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:4137630/150‑1 (MQ=11)
                                                                                                                    tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                               >  2:5110497/1‑150 (MQ=32)
                                                                                                                       aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgcaggattcgaacctgcgaccaagtgagtacaagtcaactgctctaccaactgggctaacgacccgacgtggtgggtgatgacgggatcgaac                            <  2:4857347‑M2/150‑109 (MQ=255)
                                                                                                                         gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTCTGAAAGATTAGAACCTGCAGCAGGTTTGAGTTGAGCGGAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCGGAAGGGCGAGCGAAGCGAGTCATc                          >  1:1843078/1‑150 (MQ=11)
                                                                                                                          aCCCACCACTGAAAACAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGACACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATcc                         <  2:4445137/150‑1 (MQ=11)
                                                                                                                             caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCATAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGTGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                      >  2:1641537/1‑150 (MQ=17)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAATCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:2005787/1‑150 (MQ=17)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:3410600/1‑150 (MQ=32)
                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2094423/1‑150 (MQ=32)
                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2535318/1‑150 (MQ=32)
                                                                                                                                 atagcaaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                  >  1:2693030/6‑150 (MQ=255)
                                                                                                                                   tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                >  1:461630/1‑150 (MQ=17)
                                                                                                                                     tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca              >  2:2871559/1‑150 (MQ=32)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAACCTGCACGACccac             >  1:4071702/1‑150 (MQ=17)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac             >  2:2682603/1‑150 (MQ=32)
                                                                                                                                       aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:155835/1‑150 (MQ=32)
                                                                                                                                       aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  2:665264/1‑150 (MQ=32)
                                                                                                                                        aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:4992355/150‑1 (MQ=17)
                                                                                                                                         aaaGCGCCCTAATGGCTATTTTTTGCTATCTGCGATACTCAAATATTCTAACCTGATGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa          <  1:4116458/150‑1 (MQ=11)
                                                                                                                                          atgcgcCCTTAACTCGATTTTTTGCTTTCTGCGATCCTCAACGACACGAAACTGCAGCTGGTTTGAGTTGTTGAGAGCGTAACTACGTAGGCTCTCGCTGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         <  2:5456368/148‑1 (MQ=2)
                                                                                                                                              gcCCAAAGGGCGTTTTTTTGCTATCTGCGATAGTCAAAGATTCGAACCTGAAGCTGTTTTTTGTTGAGGGCTGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa     <  2:1575397/150‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                                                 
GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781523‑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: