Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 66,533 3540 bp→CCCGAGACCCGCACGCGTCGACCTGCA [araD]araB [araD], araA, araB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 66533 70072 3540 226 [0] [0] 249 [araD]–araB [araD],araA,araB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 665320 (0.000)162 (0.930)
+CCCGAGACCCGCACGCGTCGACCTGCA
88/130 0.3 100% coding (19/696 nt) araD L‑ribulose‑5‑phosphate 4‑epimerase
?NC_000913 70073 = 0 (0.000)intergenic (‑25/‑314) araB/araC L‑ribulokinase/transcriptional regulator for ara operon

Only 100 of 226 total aligned reads displayed.

CCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAC  >  NC_000913/66441‑66622
                                                                                           |                                                                                          
 cgcgCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcc                                                                                          <  1:2926597‑M1/93‑3 (MQ=255)
  gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGccc                                                                                         <  1:1415285‑M1/93‑4 (MQ=255)
   cgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                        >  1:4205276‑M1/1‑89 (MQ=255)
   cgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                        <  2:1133374‑M1/93‑5 (MQ=255)
   cgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                        >  2:1648188‑M1/1‑89 (MQ=255)
     cTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgag                                                                                      >  2:879128‑M1/1‑87 (MQ=255)
       cgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagac                                                                                    <  1:4484797‑M1/93‑9 (MQ=255)
         cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagaccc                                                                                  <  1:5031492‑M1/93‑11 (MQ=255)
         cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagaccc                                                                                  <  1:3876365‑M1/93‑11 (MQ=255)
           aTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgc                                                                                >  2:2989835‑M1/1‑81 (MQ=255)
            tCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgca                                                                               >  2:3865017‑M1/1‑80 (MQ=255)
            tCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgca                                                                               <  1:2291026‑M1/93‑14 (MQ=255)
               aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                            >  2:1266597‑M1/1‑77 (MQ=255)
               aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                            >  1:4459403‑M1/1‑77 (MQ=255)
               aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                            <  1:3325126‑M1/93‑17 (MQ=255)
                cGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcg                                                                           >  2:1921739‑M1/1‑76 (MQ=255)
                cGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcg                                                                           >  1:4511987‑M1/1‑76 (MQ=255)
                 ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                          >  2:2173642‑M1/1‑75 (MQ=255)
                 ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                          >  1:3405420‑M1/1‑75 (MQ=255)
                  gcgcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtc                                                                         <  1:634280‑M1/93‑20 (MQ=255)
                      tGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacc                                                                     <  2:2376207‑M1/93‑24 (MQ=255)
                      tGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacc                                                                     <  1:5064453‑M1/93‑24 (MQ=255)
                      tGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacc                                                                     <  2:1963413‑M1/93‑24 (MQ=255)
                        aCGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctg                                                                   <  2:896534‑M1/93‑26 (MQ=255)
                             gCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggt                                                              <  2:5119137‑M1/93‑31 (MQ=255)
                             gCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggt                                                              <  1:488092‑M1/93‑31 (MQ=255)
                                ccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatg                                                           >  2:5805128‑M1/1‑60 (MQ=255)
                                ccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatg                                                           <  1:2979686‑M1/93‑34 (MQ=255)
                                ccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcacgtatg                                                           >  1:239085‑M1/1‑60 (MQ=255)
                                ccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGacctagacccgcccgcgtcgacctgcagttatg                                                           >  2:2596787‑M1/1‑60 (MQ=255)
                                 cATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgg                                                          <  1:1464878‑M1/93‑35 (MQ=255)
                                    gtgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggaga                                                       >  2:1167741‑M1/1‑56 (MQ=255)
                                      gAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaa                                                     >  2:2621924‑M1/1‑54 (MQ=255)
                                       aGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaac                                                    <  1:560290‑M1/93‑41 (MQ=255)
                                        gCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaaca                                                   <  1:1762987‑M1/93‑42 (MQ=255)
                                         gcttaccagGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacag                                                  <  2:6081195‑M1/89‑43 (MQ=255)
                                          gTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                                 >  1:2503108‑M1/1‑50 (MQ=255)
                                           tGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagta                                                <  1:3163306‑M1/93‑45 (MQ=255)
                                           tGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCCGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagta                                                <  2:1489743‑M1/93‑45 (MQ=255)
                                            gACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtag                                               <  2:4479721‑M1/93‑46 (MQ=255)
                                             aCCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtaga                                              >  2:1294063‑M1/1‑47 (MQ=255)
                                               cAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGgccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagaga                                            >  2:3987199‑M1/1‑45 (MQ=255)
                                               cAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagaga                                            <  2:4919165‑M1/93‑49 (MQ=255)
                                               cAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagaga                                            <  1:3732684‑M1/93‑49 (MQ=255)
                                                aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagag                                           <  2:4449132‑M1/93‑50 (MQ=255)
                                                aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaaaagtagagag                                           >  2:2108134‑M1/1‑44 (MQ=255)
                                                 ggTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagt                                          <  1:836948‑M1/93‑51 (MQ=255)
                                                  gTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagtt                                         <  1:5255619‑M1/93‑52 (MQ=255)
                                                  gTTGTGTTTTGGCAGCGCCAGGGTGGCTTCTAATACCCGGCGcccgagcaccgcacgcgtcgacctgcaggtatggagcaacagtagagagtt                                         >  1:2903052‑M1/1‑42 (MQ=255)
                                                   gtgtgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttg                                        <  2:652507‑M1/92‑53 (MQ=255)
                                                    tgtgTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgc                                       <  2:4057367‑M1/93‑54 (MQ=255)
                                                      tgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcga                                     <  2:4301754‑M1/93‑56 (MQ=255)
                                                      tgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcga                                     >  1:2071335‑M1/1‑38 (MQ=255)
                                                        ttttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgata                                   <  1:701506‑M1/93‑58 (MQ=255)
                                                        ttttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgata                                   >  2:3300381‑M1/1‑36 (MQ=255)
                                                          ttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaa                                 >  1:4142074‑M1/1‑34 (MQ=255)
                                                          ttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaa                                 >  2:6342733‑M1/1‑34 (MQ=255)
                                                            ggCAGCGCCAGGTTGTCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaa                               <  1:2595150‑M1/93‑62 (MQ=255)
                                                            ggCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaa                               <  2:5249910‑M1/93‑62 (MQ=255)
                                                            ggCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaa                               >  2:5905226‑M1/1‑32 (MQ=255)
                                                            ggCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggcatggagaaacagtagagagttgcggtaaaaa                               >  2:2760161‑M1/1‑32 (MQ=255)
                                                             gCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaag                              <  2:4638398‑M1/93‑63 (MQ=255)
                                                               aGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                            <  1:6076886‑M1/93‑65 (MQ=255)
                                                               aGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                            >  2:3705434‑M1/1‑29 (MQ=255)
                                                                gcgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgt                           <  1:4811451‑M1/93‑66 (MQ=255)
                                                                gcgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgt                           <  1:1261383‑M1/93‑66 (MQ=255)
                                                                  gcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtca                         <  2:2256316‑M1/93‑68 (MQ=255)
                                                                  gcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtca                         >  2:543285‑M1/1‑26 (MQ=255)
                                                                  gcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtca                         >  1:2658986‑M1/1‑26 (MQ=255)
                                                                   ccAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcag                        <  2:6355237‑M1/93‑69 (MQ=255)
                                                                     aGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggt                      <  1:2353888‑M1/93‑71 (MQ=255)
                                                                      ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                     >  2:5777557‑M1/1‑22 (MQ=255)
                                                                       gTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtag                    <  2:76399‑M1/93‑73 (MQ=255)
                                                                       gTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtag                    >  1:4334642‑M1/1‑21 (MQ=255)
                                                                        ttGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagg                   >  1:1975850‑M1/1‑20 (MQ=255)
                                                                         tGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                  <  1:2737976‑M1/93‑75 (MQ=255)
                                                                         tGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                  <  1:6209758‑M1/93‑75 (MQ=255)
                                                                         tGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                  <  2:4417979‑M1/93‑75 (MQ=255)
                                                                          ggCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggat                 <  2:5667351‑M1/93‑76 (MQ=255)
                                                                          ggCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggat                 >  1:4379794‑M1/1‑18 (MQ=255)
                                                                          ggCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggat                 <  1:922785‑M1/93‑76 (MQ=255)
                                                                           gCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatc                <  1:5266174‑M1/93‑77 (MQ=255)
                                                                           gCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatc                <  1:2552558‑M1/93‑77 (MQ=255)
                                                                           gCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatc                <  2:5295265‑M1/93‑77 (MQ=255)
                                                                           gCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatc                <  1:6377065‑M1/93‑77 (MQ=255)
                                                                             ttCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccg              >  1:5757823‑M1/1‑15 (MQ=255)
                                                                              tCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgc             >  2:3646638‑M1/1‑14 (MQ=255)
                                                                              tCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgc             >  1:1859491‑M1/1‑14 (MQ=255)
                                                                               cTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgct            >  2:280737‑M1/1‑13 (MQ=255)
                                                                               cTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgct            <  2:2715631‑M1/93‑81 (MQ=255)
                                                                               cTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgct            <  1:5668749‑M1/93‑81 (MQ=255)
                                                                                 aaTACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaa          <  2:4787062‑M1/93‑83 (MQ=255)
                                                                                  aTACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaat         >  2:1794799‑M1/1‑10 (MQ=255)
                                                                                   tACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatc        >  1:1000921‑M1/1‑9 (MQ=255)
                                                                                     cccggcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgacaaaaagcgtcaggtaggatccgctaatctt      <  2:3386271‑M1/90‑87 (MQ=255)
                                                                                      ccggcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatctta     <  2:3155977‑M1/91‑88 (MQ=255)
                                                                                      ccggcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatctta     <  2:2920340‑M1/91‑88 (MQ=255)
                                                                                       cggcGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttat    <  1:5876428‑M1/92‑89 (MQ=255)
                                                                                         gcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcatgtaggatccgctaatcttatgg  <  2:4843573‑M1/93‑91 (MQ=255)
                                                                                         gcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatgg  >  2:1349608‑M1/1‑3 (MQ=255)
                                                                                           |                                                                                          
CCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAC  >  NC_000913/66441‑66622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: