Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,670 A→T 6.0% intergenic (+18/‑415) 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,6700AT 6.0% 603.3 / 29.5 233intergenic (+18/‑415)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (117/96);  new base T (4/10);  total (126/107)
Fisher's exact test for biased strand distribution p-value = 9.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

Only 100 of 236 total aligned reads displayed.

GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781521‑781817
                                                                                                                                                      |                                                                                                                                                   
gAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTa                                                                                                                                                     >  1:5470101/1‑150 (MQ=255)
   taaTCCTGCACGACCCAACATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAg                                                                                                                                                  <  2:5072270/148‑1 (MQ=255)
   tCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAg                                                                                                                                                  >  1:1319748/1‑150 (MQ=255)
         tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTCAAGGCGCtt                                                                                                                                            >  2:3389496/1‑150 (MQ=255)
         tGCAAGACCCACCATCCTTATTGATTTAGGCAGCATAAGCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtt                                                                                                                                            <  2:4496273/150‑1 (MQ=17)
          gCACGACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttt                                                                                                                                           <  1:1197788/150‑1 (MQ=255)
          gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttt                                                                                                                                           >  1:4495136/1‑150 (MQ=255)
           cccGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttt                                                                                                                                          <  1:3117922/148‑1 (MQ=255)
           cACGACCGACCATCCTTTTTTATTTATGCAGCCTAATCCCGGTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtttt                                                                                                                                          <  2:5344934/150‑1 (MQ=14)
           cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttt                                                                                                                                          <  2:1936042/150‑1 (MQ=255)
            aCGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCttttt                                                                                                                                         >  1:4341404/1‑150 (MQ=255)
             cGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCTTTAGGGCGCtttttt                                                                                                                                        <  2:1786100/150‑1 (MQ=255)
                cccACCATCCTGAATGATTAAGGCAGCATACCCCCGCACGGGGTCGTTATCTCAGTGGTTAGAGCAGTTCACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAA‑GCGCCTTTAGGGCGCTTTTTTTCt                                                                                                                                     <  1:1147844/150‑1 (MQ=17)
                cacaccaTCCTGTTTGATTAATGCATCTTAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCt                                                                                                                                     <  2:2062591/148‑1 (MQ=17)
                 ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAA‑GCGCCTTTAGGGCGCTTTTTTacat                                                                                                                                    >  2:5242796‑M2/1‑146 (MQ=255)
                  caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  1:4633906/1‑150 (MQ=255)
                  caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                   >  2:4069015/1‑150 (MQ=255)
                   accaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTacattg                                                                                                                                  >  2:2842051‑M2/1‑144 (MQ=255)
                    ccaTCCTGAATGATTCAGGCAGCATACTCCCGCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCt                                                                                                                                 <  2:3285419/150‑1 (MQ=255)
                      cgcccgAAAGATTAAGGCAGCATACTCCCGCAAGGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:418786/145‑1 (MQ=17)
                      aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  1:1266848/150‑1 (MQ=255)
                        ccggcaagATTTAGGCCGCCTACGCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                                             <  2:3804495/143‑1 (MQ=11)
                         cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATAAGCGAt                                                                                                                            >  2:4834408/1‑150 (MQ=255)
                          tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                           >  2:2633753/1‑150 (MQ=255)
                           gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATAc                                                                                                                          >  2:4570804/1‑150 (MQ=255)
                             aTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCGAAAGGCGCTTTTTTGCTATCTGCGATACTc                                                                                                                        <  2:3087409/150‑1 (MQ=255)
                              tGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCa                                                                                                                       <  2:2265849/150‑1 (MQ=255)
                                aTTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:4426116/150‑1 (MQ=255)
                                aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                     <  2:782611/150‑1 (MQ=255)
                                      gcagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                               >  2:1520730/1‑150 (MQ=255)
                                        agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                             >  2:1341104/1‑150 (MQ=255)
                                          caTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACCAAcc                                                                                                           >  2:4113771/1‑150 (MQ=255)
                                           aTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCtttttt                                                                                                                                        <  2:3465197/120‑1 (MQ=255)
                                            tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                         <  1:1198174/150‑1 (MQ=255)
                                               tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      >  2:3960605/1‑150 (MQ=39)
                                                      aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:4678950/150‑1 (MQ=33)
                                                      aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:662897/150‑1 (MQ=33)
                                                      aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  2:4018495/150‑1 (MQ=33)
                                                        gggTCGTTTGTCAAGTTGGTAGAGCTGTTGAATTTTAAACAATTGGCCGACGGTTCGAAGCCTGCACGACCCACCAATTTAAAAAA‑GCGCCCAAACGCCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  2:4632020/150‑1 (MQ=11)
                                                        gggTCGGTAGCTGAGTTTGTTGTTCTGTTGTCTTTTAATCTATTGGTCGCATGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  1:1413928/150‑1 (MQ=11)
                                                          gTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                           <  1:5209880/150‑1 (MQ=17)
                                                             gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                        >  1:2000925/1‑150 (MQ=32)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            <  1:5292019/129‑1 (MQ=32)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAc                                                                                                            >  2:5292019/1‑129 (MQ=32)
                                                              ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                       >  1:4995909/1‑150 (MQ=32)
                                                               tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                                      <  1:4289101/150‑1 (MQ=32)
                                                                 gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                    <  2:1319748/150‑1 (MQ=32)
                                                                  cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   <  1:818955/150‑1 (MQ=32)
                                                                  cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   >  1:5063542/1‑150 (MQ=32)
                                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                  <  1:2633628/150‑1 (MQ=32)
                                                                      gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCAAGGATTCGAACCTGCAGCAGGTGTGTGTTGAGCGCAGCGa                                                                               >  2:599863/1‑150 (MQ=255)
                                                                      gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                               >  2:3279174/1‑150 (MQ=17)
                                                                        tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                             <  1:4743045/150‑1 (MQ=32)
                                                                          gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                           >  1:3438665/1‑150 (MQ=17)
                                                                           tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                          <  2:4542379/150‑1 (MQ=32)
                                                                               gaaGTTTAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCTAAACAACGGa                                                                      <  2:1869540/148‑1 (MQ=255)
                                                                                 aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    >  1:5101210/1‑150 (MQ=32)
                                                                                  gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                   >  1:3053346/1‑150 (MQ=32)
                                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCAGAAGCAGTTAGGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                                  >  1:2219301/1‑150 (MQ=11)
                                                                                    tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCCAACCTGCAGCAGGTTTGAGAGGAGCGCAGAGAAACAAAGGAGCCGc                                                                 >  1:1696659/1‑150 (MQ=11)
                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                             >  2:1988189/1‑150 (MQ=32)
                                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACGCAAAGATACGATCCTGCACCAGTTTTCAGATGATCGCAGCGAAGCGCACGAGCCGCTcgc                                                             >  1:3612563/1‑150 (MQ=11)
                                                                                             aTCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGAAGTTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:3233045/150‑1 (MQ=11)
                                                                                             aTCAATTGGCCGCAGGTTCGAACCCTGCCCGACCCACCACTGTAAAAAC‑GCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                        <  2:2360506/150‑1 (MQ=11)
                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGATCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  2:2296445/1‑150 (MQ=17)
                                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       >  1:404248/1‑150 (MQ=17)
                                                                                                aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                     >  2:1472973/1‑150 (MQ=32)
                                                                                                  ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCACAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   >  2:381594/1‑150 (MQ=17)
                                                                                                  ttGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   <  2:4120873/150‑1 (MQ=17)
                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTTAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               <  2:1110520/150‑1 (MQ=12)
                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:419337/1‑150 (MQ=32)
                                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGCAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                               >  2:4223706/1‑150 (MQ=17)
                                                                                                       cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATGTGCGATACTCAAATATTCGTACCTGCACCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg                                              >  1:4378854/1‑150 (MQ=11)
                                                                                                            gTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                         >  2:2923622/1‑150 (MQ=32)
                                                                                                             ttCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                        >  1:3682750/1‑150 (MQ=32)
                                                                                                                gAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                     >  1:3017862/1‑150 (MQ=32)
                                                                                                                  aTCCTGCACGACCCACCAATTTAAAAAT‑GCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                   <  2:754114/150‑1 (MQ=255)
                                                                                                                     cTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:3032598/150‑1 (MQ=32)
                                                                                                                     cTGCACGACCAACCAATGTAAAAAA‑GCGCCCTAATGGCGCTTTTTTGATATATGCGATACTAAAAGATTCTAACCTGCAGCTGGTTTTTGTTGAGCGCAGCTAAACAACGGATCCGATCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                <  2:4137630/150‑1 (MQ=11)
                                                                                                                      tGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                               >  2:5110497/1‑150 (MQ=32)
                                                                                                                         aCGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgcaggattcgaacctgcgaccaagtgagtacaagtcaactgctctaccaactgggctaacgacccgacgtggtgggtgatgacgggatcgaac                            <  2:4857347‑M2/150‑109 (MQ=255)
                                                                                                                         aCGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  2:4287498/150‑1 (MQ=32)
                                                                                                                         aCGACCCA‑CAATTTAAAAAACGCCCCCTACGGGGGTTTTTTTGCTTTATGCTATACTCAAAGATTCCACCCTGCTGCTGTTTTGTGTTGCGGGCTCCCAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                            <  1:4954983/150‑1 (MQ=2)
                                                                                                                          cgccccACCACTGTAAAAAA‑CCGCCCTAAAGGCGCTTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCTGTTTTTAGTTGTGCGCAGCGAAACAACGGAGCCGCTCGCTGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                           <  2:1022029/147‑1 (MQ=11)
                                                                                                                           gACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTCTGAAAGATTAGAACCTGCAGCAGGTTTGAGTTGAGCGGAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCGGAAGGGCGAGCGAAGCGAGTCATc                          >  1:1843078/1‑150 (MQ=11)
                                                                                                                            aCCCACCACTGAAAACAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGACACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGAGCGAAGCGAGTCATcc                         <  2:4445137/150‑1 (MQ=11)
                                                                                                                                accaATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAATCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:2005787/1‑150 (MQ=17)
                                                                                                                                accaATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     >  1:5650/1‑150 (MQ=17)
                                                                                                                                  caATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2094423/1‑150 (MQ=32)
                                                                                                                                  caATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   >  2:2535318/1‑150 (MQ=32)
                                                                                                                                   atagcaaaaaa‑GCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                  >  1:2693030/6‑150 (MQ=255)
                                                                                                                                   aaTGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGAAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTTCACGa                  >  2:3527676/1‑150 (MQ=12)
                                                                                                                                       tAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca              >  2:2871430/1‑150 (MQ=32)
                                                                                                                                        aaaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAACCTGCACGACccac             >  1:4071702/1‑150 (MQ=17)
                                                                                                                                         aaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:155835/1‑150 (MQ=32)
                                                                                                                                         aaaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  2:665264/1‑150 (MQ=32)
                                                                                                                                          aaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:4992355/150‑1 (MQ=17)
                                                                                                                                          aaaa‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:1634060/150‑1 (MQ=32)
                                                                                                                                           aaa‑GCGCCCTAATGGCTATTTTTTGCTATCTGCGATACTCAAATATTCTAACCTGATGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa          <  1:4116458/150‑1 (MQ=11)
                                                                                                                                             atgcgcCCTTAACTCGATTTTTTGCTTTCTGCGATCCTCAACGACACGAAACTGCAGCTGGTTTGAGTTGTTGAGAGCGTAACTACGTAGGCTCTCGCTGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         <  2:5456368/148‑1 (MQ=2)
                                                                                                                                                      |                                                                                                                                                   
GAGTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAA‑GCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781521‑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: ‑