Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AP012306 1,438,342 Δ2,454 bp [pykF][sufE] [pykF], lpp, ynhG, [sufE]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AP012306 1438342 1440795 2454 365 [0] [0] 364 [pykF]–[sufE] [pykF],lpp,ynhG,[sufE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP012306 = 14383410 (0.000)352 (1.060) 199/294 0.1 100% coding (1061/1413 nt) pykF pyruvate kinase I
?AP012306 1440796 = 0 (0.000)coding (78/417 nt) sufE sulfur acceptor protein

Only 100 of 364 total aligned reads displayed.

GCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAATTCTGGCAGACGCTGGCCCAGCTCAATAATGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATG  >  AP012306/1440647‑1440945
                                                                                                                                                     |                                                                                                                                                     
ggtactgacgcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatg                                                                                                                                                       >  1:498574‑M2/150‑150 (MQ=255)
ggtactgacgcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatg                                                                                                                                                       >  1:1534740‑M2/150‑150 (MQ=255)
  tactgacgcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTa                                                                                                                                                     >  1:1024574‑M2/148‑150 (MQ=255)
   actgacgcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTag                                                                                                                                                    >  1:1829147‑M2/147‑150 (MQ=255)
        cgcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGata                                                                                                                                               <  1:2626331‑M2/9‑1 (MQ=255)
         gcagtgatgctgtctggtgaatgcgcaaaaggtaaatgcccgctggaagcgttttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGatat                                                                                                                                              <  2:743203‑M2/10‑1 (MQ=255)
         gcagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgaggtcaacaatgacaaccgtaaactgagcatGTAGAGatat                                                                                                                                              >  2:1233190‑M2/141‑150 (MQ=255)
          cagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAtt                                                                                                                                             <  1:1879726‑M2/11‑1 (MQ=255)
          cagtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAtt                                                                                                                                             <  2:256442‑M2/11‑1 (MQ=255)
           agtgatgctgtctggtgaatccgcaaaaggtaaatgcccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAttt                                                                                                                                            <  2:3562189‑M2/12‑1 (MQ=255)
           agtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAttt                                                                                                                                            <  1:1659303‑M2/12‑1 (MQ=255)
           agtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAttt                                                                                                                                            <  1:2904480‑M2/12‑1 (MQ=255)
           agtgatgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATAttt                                                                                                                                            >  1:3678614‑M2/139‑150 (MQ=255)
               atgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTctct                                                                                                                                        <  2:1341615‑M2/16‑1 (MQ=255)
               atgctgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTctct                                                                                                                                        <  2:4090949‑M2/16‑1 (MQ=255)
                   tgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTccc                                                                                                                                    >  2:657912‑M2/131‑150 (MQ=255)
                   tgtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTccc                                                                                                                                    >  2:3168201‑M2/131‑150 (MQ=255)
                    gtctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCa                                                                                                                                   >  2:51204‑M2/130‑150 (MQ=255)
                     tctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAg                                                                                                                                  <  2:3878672‑M2/22‑1 (MQ=255)
                     tctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAg                                                                                                                                  >  2:1786581‑M2/129‑150 (MQ=255)
                      ctggtgaatccgcaaaaggtgaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGt                                                                                                                                 <  2:1683269‑M2/23‑1 (MQ=255)
                      ctggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGt                                                                                                                                 >  1:1636831‑M2/128‑150 (MQ=255)
                       tggtgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGtt                                                                                                                                <  2:1560151‑M2/24‑1 (MQ=255)
                          tgaatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGgc                                                                                                                             <  2:3414335‑M2/27‑1 (MQ=255)
                            aatccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGgcgc                                                                                                                           >  2:4446129‑M2/122‑150 (MQ=255)
                             atccgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCa                                                                                                                          >  2:3594821‑M2/121‑150 (MQ=255)
                                cgcaaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCCGTTGGCGCAGCg                                                                                                                       <  2:4555784‑M2/33‑1 (MQ=255)
                                   aaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTaa                                                                                                                    >  1:3267564‑M2/115‑150 (MQ=255)
                                   aaaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTaa                                                                                                                    <  2:4221454‑M2/36‑1 (MQ=255)
                                    aaaggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTaaa                                                                                                                   >  1:2776100‑M2/114‑150 (MQ=255)
                                      aggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTaaaaa                                                                                                                 <  1:4580732‑M2/39‑1 (MQ=255)
                                       ggtaaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTaaaaaa                                                                                                                <  1:2247808‑M2/40‑1 (MQ=255)
                                         taaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAAtt                                                                                                              <  1:2077092‑M2/42‑1 (MQ=255)
                                          aaatacccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTa                                                                                                             <  2:363980‑M2/43‑1 (MQ=255)
                                               cccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACgcag                                                                                                        <  1:2712652‑M2/48‑1 (MQ=255)
                                                ccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACgcagc                                                                                                       <  1:3407231‑M2/49‑1 (MQ=255)
                                                ccgctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACgcagc                                                                                                       <  1:3166814‑M2/49‑1 (MQ=255)
                                                  gctggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCaa                                                                                                     <  1:798613‑M2/51‑1 (MQ=255)
                                                    tggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACt                                                                                                   <  2:2981029‑M2/53‑1 (MQ=255)
                                                     ggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACtt                                                                                                  >  1:4274257‑M2/97‑150 (MQ=255)
                                                     ggaagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACtt                                                                                                  >  1:688379‑M2/97‑150 (MQ=255)
                                                       aagcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACtttt                                                                                                <  2:4371931‑M2/56‑1 (MQ=255)
                                                        agcggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTc                                                                                               <  1:2453547‑M2/57‑1 (MQ=255)
                                                          cggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCtt                                                                                             >  2:2534042‑M2/92‑150 (MQ=255)
                                                          cggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCtt                                                                                             <  2:3222920‑M2/59‑1 (MQ=255)
                                                          cggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCtt                                                                                             <  2:1263951‑M2/59‑1 (MQ=255)
                                                           ggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCttt                                                                                            >  1:105452‑M2/91‑150 (MQ=255)
                                                           ggtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCttt                                                                                            <  1:2134259‑M2/60‑1 (MQ=255)
                                                            gtttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTa                                                                                           >  1:3321606‑M2/90‑150 (MQ=255)
                                                             tttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTAt                                                                                          <  1:1768520‑M2/62‑1 (MQ=255)
                                                              ttctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATc                                                                                         <  2:1289180‑M2/63‑1 (MQ=255)
                                                                ctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCg                                                                                       <  1:3541755‑M2/65‑1 (MQ=255)
                                                                ctatcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCg                                                                                       >  1:4174413‑M2/86‑150 (MQ=255)
                                                                    catggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCaa                                                                                   <  2:3061221‑M2/69‑1 (MQ=255)
                                                                    catggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCaa                                                                                   >  1:2405921‑M2/82‑150 (MQ=255)
                                                                      tggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATa                                                                                 <  1:896579‑M2/71‑1 (MQ=255)
                                                                       ggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATaa                                                                                <  2:2486625‑M2/72‑1 (MQ=255)
                                                                       ggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATaa                                                                                >  2:1877208‑M2/79‑150 (MQ=255)
                                                                         cgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAg                                                                              <  2:1066023‑M2/74‑1 (MQ=255)
                                                                         cgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAg                                                                              <  1:2185826‑M2/74‑1 (MQ=255)
                                                                           accatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGcc                                                                            >  2:3388825‑M2/75‑150 (MQ=255)
                                                                             catctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCAt                                                                          <  1:1902155‑M2/78‑1 (MQ=255)
                                                                               tctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAg                                                                        <  2:884304‑M2/80‑1 (MQ=255)
                                                                                  gcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGc                                                                     <  2:4161916‑M2/83‑1 (MQ=255)
                                                                                  gcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGc                                                                     <  1:3442895‑M2/83‑1 (MQ=255)
                                                                                     aacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTc                                                                  <  2:1984206‑M2/86‑1 (MQ=255)
                                                                                             gaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTAt                                                          >  2:3672721‑M2/57‑150 (MQ=255)
                                                                                                   gtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGc                                                    >  2:2676413‑M2/51‑150 (MQ=255)
                                                                                                    tgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCa                                                   <  2:1140354‑M2/101‑1 (MQ=255)
                                                                                                          acagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGt                                             >  2:581224‑M2/44‑150 (MQ=255)
                                                                                                           cagccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTg                                            >  2:2916602‑M2/43‑150 (MQ=255)
                                                                                                              ccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAAt                                         <  1:116769‑M2/111‑1 (MQ=255)
                                                                                                              ccgtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAAt                                         >  1:4521250‑M2/40‑150 (MQ=255)
                                                                                                                gtctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATAc                                       <  2:2124768‑M2/113‑1 (MQ=255)
                                                                                                                 tctcgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACg                                      <  1:3765092‑M2/114‑1 (MQ=255)
                                                                                                                    cgagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTg                                   >  2:2060667‑M2/34‑150 (MQ=255)
                                                                                                                     gagttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGc                                  <  2:2183809‑M2/118‑1 (MQ=255)
                                                                                                                      agttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCa                                 >  1:1571533‑M2/32‑150 (MQ=255)
                                                                                                                        ttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAgg                               <  1:982868‑M2/121‑1 (MQ=255)
                                                                                                                        ttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAgg                               <  1:1245371‑M2/121‑1 (MQ=255)
                                                                                                                        ttcaacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAgg                               >  2:229384‑M2/30‑150 (MQ=255)
                                                                                                                           aacaatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCg                            <  1:1993308‑M2/124‑1 (MQ=255)
                                                                                                                             caatgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGt                          <  2:442907‑M2/126‑1 (MQ=255)
                                                                                                                               atgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCa                        <  2:2355189‑M2/128‑1 (MQ=255)
                                                                                                                                tgacaaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCAc                       >  2:764724‑M2/22‑150 (MQ=255)
                                                                                                                                   caaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAg                    <  1:297362‑M2/132‑1 (MQ=255)
                                                                                                                                    aaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGa                   >  2:1488202‑M2/18‑150 (MQ=255)
                                                                                                                                    aaccgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGa                   >  1:3626687‑M2/18‑150 (MQ=255)
                                                                                                                                       cgtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGa                >  2:2550470‑M2/15‑150 (MQ=255)
                                                                                                                                        gtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGAt               >  2:3840900‑M2/14‑150 (MQ=255)
                                                                                                                                        gtaaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGAt               <  2:2763663‑M2/137‑1 (MQ=255)
                                                                                                                                         taaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATc              <  1:4056305‑M2/138‑1 (MQ=255)
                                                                                                                                          aaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATcc             >  1:1505275‑M2/12‑150 (MQ=255)
                                                                                                                                          aaactgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCATCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATcc             >  1:3905975‑M2/12‑150 (MQ=255)
                                                                                                                                           aactgagcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCa            <  2:3311140‑M2/140‑1 (MQ=255)
                                                                                                                                              tgcgcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCActt         >  2:399444‑M2/8‑150 (MQ=255)
                                                                                                                                                 gcatGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCActtctt      <  1:144243‑M2/146‑1 (MQ=255)
                                                                                                                                                  catGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACttcttc     >  1:4397233‑M2/4‑150 (MQ=255)
                                                                                                                                                   atGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCa    >  2:2581256‑M2/3‑150 (MQ=255)
                                                                                                                                                     gTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATg  <  2:1408129/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAATTCTGGCAGACGCTGGCCCAGCTCAATAATGTAGAGATATTTCTCTTCCCAGTTGGCGCAGCGTAAAAAATTACGCAGCAACTTTTCTTTATCCGGCAATAAAGCCATAGTGCCTCCCTGTTATCCCAGCAAACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATG  >  AP012306/1440647‑1440945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: