Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AE015451 1,709,370 2640 bp→CGGCCGC [ppc] [ppc]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AE015451 1709370 1712009 2640 38 [0] [0] 38 [ppc] [ppc]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 17093690 (0.000)38 (0.550)
+CGGCCGC
38/586 0.6 100% intergenic (‑136/‑13) PP_1504/ppc conserved protein of unknown function/Phosphoenolpyruvate carboxylase
?AE015451 1712010 = 0 (0.000)intergenic (+1/‑257) ppc/adk Phosphoenolpyruvate carboxylase/Adenylate kinase

CGATCAGCGCGGCCAGGCGCGCTTCGGACTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGCCCGTAGCGAGGGCCTTCAACCAGTTCCTCAACCTGGCCAACATGGCCGAGCAGTACCAGTTGATCCGCCGCCGTGACGCCGA  >  AE015451/1709083‑1709662
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                     
cGATCAGCGCGGCCAGGCGCGCTTCGGACTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaag                                                                                                                                                                                                                                                                                         >  3:482569‑M1/1‑287 (MQ=255)
  aTCAGCGCGGCCAGGCGCGCTTCGGACTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcc                                                                                                                                                                                                                                                                                       <  1:482567‑M1/301‑17 (MQ=255)
    cAGCGCGGCCAGGCGCGCTTCGGACTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccg                                                                                                                                                                                                                                                                                     <  2:423032‑M1/301‑19 (MQ=255)
       gggggccaGGCGCGCGTCGGATTCAAGGTCAAGGGGAAGGACAGGACAGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGGGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggt                                                                                                                                                                                                                                                                                  <  2:391907‑M1/298‑22 (MQ=255)
                    gcgtgtcgggttgggtccaggggcaggacaggaccGGTACGGCTGAGGGTACCGGCGGGATTCTGTGCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgac                                                                                                                                                                                                                                                                     <  2:290651‑M1/288‑35 (MQ=255)
                            cTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGAGACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccataa                                                                                                                                                                                                                                                             >  2:504451‑M1/1‑259 (MQ=255)
                                 ggtccaggtGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGCGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacagatcggaagatc                                                                                                                                                                                                                                                        >  4:290719‑M1/1‑254 (MQ=255)
                                          tcaGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATTACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaag                                                                                                                                                                                                                                               <  4:504593‑M1/300‑57 (MQ=255)
                                              gACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacg                                                                                                                                                                                                                                           >  3:67241‑M1/1‑241 (MQ=255)
                                                   gagcggtACGGGTGAGGGTATCGGCGGCATGCTGTTCCCAGGGCTGGGCGAGTTCTAGGATGTTCATGATGACCTCTGTTGACGGCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggag                                                                                                                                                                                                                                      <  4:598398‑M1/298‑66 (MQ=255)
                                                   gACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggttacggccatacaggcccggcacaaggacggggag                                                                                                                                                                                                                                      >  3:25765‑M1/1‑236 (MQ=255)
                                                         tACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatttctggtgacggccatacaggccctgccccaggacggggacggcggg                                                                                                                                                                                                                                >  2:10457‑M1/1‑230 (MQ=255)
                                                         tACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggt                                                                                                                                                                                                                                <  3:254971‑M1/301‑72 (MQ=255)
                                                          aCGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggtt                                                                                                                                                                                                                               <  1:67239‑M1/301‑73 (MQ=255)
                                                                                                             gATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGGATTGAAGGGGAGTGCTGTGTTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgca                                                                                                                                                                            <  4:10457‑M1/301‑124 (MQ=255)
                                                                                                                  tCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaacttt                                                                                                                                                                       >  3:625064‑M1/1‑173 (MQ=255)
                                                                                                                               tGTTGACGCCTGGAAAGGCTGCACTGCATGGGTGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtct                                                                                                                                                          <  4:203739‑M1/301‑142 (MQ=255)
                                                                                                                                   gACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggct                                                                                                                                                      <  1:625062‑M1/301‑146 (MQ=255)
                                                                                                                                       ccTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccg                                                                                                                                                  >  3:244273‑M1/1‑152 (MQ=255)
                                                                                                                                                              gAGTGCTGTGCTGTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGTGCCGGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagcc                                                                                                                           <  1:25765‑M1/301‑173 (MQ=255)
                                                                                                                                                                       ctcttccgatctgtttggtttggccgatcacggcTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccg                                                                                                                  <  1:11115‑M1/289‑182 (MQ=255)
                                                                                                                                                                        ctTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccga                                                                                                                 >  4:340825‑M1/1‑119 (MQ=255)
                                                                                                                                                                            cACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttt                                                                                                             >  3:41353‑M1/1‑115 (MQ=255)
                                                                                                                                                                                   gtttggtttgGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgagatcggaagag                                                                                                      >  3:11115‑M1/1‑108 (MQ=255)
                                                                                                                                                                                       ggtttggCCGATCACGGTGGGCCTTTGGCGCTGGGGTGTGGGCACTCTTGGGGGATACCGCTCGCTCATGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcataccc                                                                                                  <  2:340746‑M1/301‑198 (MQ=255)
                                                                                                                                                                                        gtttggCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccc                                                                                                 >  3:580929‑M1/1‑103 (MQ=255)
                                                                                                                                                                                            ggCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattt                                                                                             >  3:420228‑M1/1‑99 (MQ=255)
                                                                                                                                                                                               cGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCGCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttcc                                                                                          <  1:580927‑M1/301‑206 (MQ=255)
                                                                                                                                                                                                     cGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgatttttcggtctt                                                                                    >  4:780179‑M1/1‑90 (MQ=255)
                                                                                                                                                                                                          tGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccc                                                                               >  3:345175‑M1/1‑85 (MQ=255)
                                                                                                                                                                                                                               gcgacactggggggataccgtgcGCTCGTGCACAAGGCGCGCTAGGACTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatccaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaat                                                          <  2:779899‑M1/292‑238 (MQ=255)
                                                                                                                                                                                                                                         ggggATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaa                                                >  4:783584‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                                                      ctgcggatctagggcgggctgAGGCTGATAGGGATGAAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaaaatttgaggagca                                   <  2:395977‑M1/291‑261 (MQ=255)
                                                                                                                                                                                                                                                                aGGGCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaaaatttgaggagcaagatcggaag                         >  4:396074‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                                                                  ggCGCGCTGAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaaaatttgaggagcacaagatgcgcgt                       <  1:41352‑M1/301‑273 (MQ=255)
                                                                                                                                                                                                                                                                         ggaGGCTGATAGGGAGGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaaaatttgaggagcacaagatgcgcgtaattctg                <  2:783304‑M1/300‑280 (MQ=255)
                                                                                                                                                                                                                                                                          gAGGCTGATAGGGATGCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccatgtaatttaatacaaaaatttgaggagcacaaggtgcgggtaattcagg               >  2:264934‑M1/1‑21 (MQ=255)
                                                                                                                                                                                                                                                                                        gtgCAACGcggccgcgggcaagcccgggtcatgtctggtgacggccatacaggcccggcacaaggacggggagggcggttgcgccaggcgcaaccgcctgcccggtcagccgcaagtggcgcattcctgcaactttgggacgcttgtctggctgccgggcgctgtgtatcttgagcagccttctgaccgatttatggtcatacccgattttccggacttggccctggtcgccgaatccattgaatttcatacaaaaatttgaggagcacaagatgcgcgtaattctgctgggagctcccgg  <  4:265000‑M1/300‑294 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                     
CGATCAGCGCGGCCAGGCGCGCTTCGGACTCCAGGTCCAGGTGCAGGACATGACCGGTACGGCTGAGGGTACCGGCGGCATTCTGTTCCCAGTGCTGGGCGAGTTCACGGATTTTCATGATGACCTCTGTTGACGCCTGGAAAGGCTGCACTGCATGGGAGTGCTGTGCTTTCACTTTAGTTTGGTTTGGCCGATCACGGCTTGCCATCGTCGCTCGGCTGTGGGCACTCTTGGGGGATACCGCTCGCTCCTGCACAGGGCGCGCTGAGGCTGATAGGGATGCAACGGAGAGAGCGCCAATGACCGATATCGATGTGCGTTTGCGTGAAGATGTCCATGTGTTGGGTGAACTGCTTGGCGAAACCATTCGCCAACAGCATGGCGATGCGTTCCTGCAGAAAATCGAGGACATCCGCCACAGTGCCAAGGCCGACCGTCGCGGCCCGGGTGAGCAGCTGAGCTCGACCCTGGCGGACCTGGCCGAAGAAGACCTGCTGCCCGTAGCGAGGGCCTTCAACCAGTTCCTCAACCTGGCCAACATGGCCGAGCAGTACCAGTTGATCCGCCGCCGTGACGCCGA  >  AE015451/1709083‑1709662

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