Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,215,179 Δ4 bp coding (900‑903/930 nt) metA → homoserine O‑transsuccinylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4215179 4215182 4 101 [0] [0] 103 metA homoserine O‑transsuccinylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42151780 (0.000)100 (1.050) 73/196 0.1 100% coding (899/930 nt) metA homoserine O‑transsuccinylase
?NC_000913 4215183 = 0 (0.000)coding (904/930 nt) metA homoserine O‑transsuccinylase

Only 100 of 101 total aligned reads displayed.

GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTA  >  NC_000913/4215081‑4215277
                                                                                                 |                                                                                                   
gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcg                                                                                                   >  2:1357742‑M1/1‑98 (MQ=255)
 aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcgg                                                                                                  >  2:1926888‑M1/1‑97 (MQ=255)
  tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcgg                                                                                                  <  2:70050‑M1/99‑4 (MQ=255)
   ccGCAAAATCCACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggca                                                                                                >  2:1777752‑M1/1‑95 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggca                                                                                                <  1:1085074‑M1/100‑6 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggca                                                                                                >  1:483761‑M1/1‑95 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggca                                                                                                >  2:2345886‑M1/1‑95 (MQ=255)
   ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggca                                                                                                >  2:2338946‑M1/1‑95 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcaca                                                                                              >  2:1419454‑M1/1‑93 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcaca                                                                                              >  2:445990‑M1/1‑93 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcaca                                                                                              >  1:1729363‑M1/1‑93 (MQ=255)
     gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcaca                                                                                              >  2:2288052‑M1/1‑93 (MQ=255)
       aaaaTCAACCGAGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatg                                                                                            <  1:445990‑M1/100‑10 (MQ=255)
       aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatg                                                                                            <  2:2034715‑M1/100‑10 (MQ=255)
       aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatg                                                                                            <  1:1841300‑M1/100‑10 (MQ=255)
        aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatga                                                                                           >  1:554810‑M1/1‑90 (MQ=255)
          aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaat                                                                                         >  1:2098210‑M1/1‑88 (MQ=255)
          aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaat                                                                                         <  2:372269‑M1/100‑13 (MQ=255)
           tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatc                                                                                        <  1:121233‑M1/100‑14 (MQ=255)
            acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatcc                                                                                       <  1:243905‑M1/100‑15 (MQ=255)
             cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatcca                                                                                      >  2:1000373‑M1/1‑85 (MQ=255)
               ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaac                                                                                    >  1:1027674‑M1/1‑83 (MQ=255)
                  cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgct                                                                                 >  1:1115206‑M1/1‑80 (MQ=255)
                    agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctgg                                                                               <  2:309858‑M1/100‑23 (MQ=255)
                    agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctgg                                                                               <  2:578782‑M1/100‑23 (MQ=255)
                      cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggat                                                                             <  1:906725‑M1/100‑25 (MQ=255)
                      cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggat                                                                             >  1:1325070‑M1/1‑76 (MQ=255)
                       gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggatt                                                                            <  2:1297595‑M1/100‑26 (MQ=255)
                        agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggatta                                                                           >  2:1977412‑M1/1‑74 (MQ=255)
                          cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaat                                                                         >  2:135829‑M1/1‑72 (MQ=255)
                          cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaat                                                                         >  1:2310194‑M1/1‑72 (MQ=255)
                           tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatc                                                                        >  1:1914714‑M1/1‑71 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttc                                                                     <  1:1700939‑M1/100‑33 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttc                                                                     >  2:254752‑M1/1‑68 (MQ=255)
                              cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttc                                                                     <  2:790677‑M1/100‑33 (MQ=255)
                               gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttct                                                                    <  2:82244‑M1/100‑34 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgt                                                                  <  2:2312732‑M1/100‑36 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgt                                                                  <  2:1861874‑M1/100‑36 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgt                                                                  <  2:1729364‑M1/100‑36 (MQ=255)
                                 aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgt                                                                  >  1:1560075‑M1/1‑65 (MQ=255)
                                  gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtg                                                                 >  1:1085986‑M1/1‑64 (MQ=255)
                                    cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgat                                                               >  2:26517‑M1/1‑62 (MQ=255)
                                     aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgat                                                               <  2:1560076‑M1/99‑39 (MQ=255)
                                        gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtc                                                           <  2:2310192‑M1/100‑43 (MQ=255)
                                         tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcg                                                          <  2:1920535‑M1/100‑44 (MQ=255)
                                           aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgat                                                        >  2:2017574‑M1/1‑55 (MQ=255)
                                            tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatc                                                       >  2:2248315‑M1/1‑54 (MQ=255)
                                             ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcg                                                      >  1:225966‑M1/1‑53 (MQ=255)
                                               aCTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgtt                                                    >  1:1127647‑M1/1‑51 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgtta                                                   <  1:114910‑M1/100‑51 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgtta                                                   <  1:688676‑M1/100‑51 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgtta                                                   >  1:924951‑M1/1‑50 (MQ=255)
                                                cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgtta                                                   >  2:688677‑M1/1‑50 (MQ=255)
                                                  gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaag                                                 >  2:1553307‑M1/1‑48 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagc                                                <  1:1290553‑M1/100‑54 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagc                                                <  1:2248317‑M1/100‑54 (MQ=255)
                                                   tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagc                                                <  1:604397‑M1/100‑54 (MQ=255)
                                                     tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcga                                              >  1:1303637‑M1/1‑45 (MQ=255)
                                                       ccAACTGGCTCAACTATTACGTTTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgatt                                            >  2:1529511‑M1/1‑43 (MQ=255)
                                                       ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcgacacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgatt                                            >  1:358783‑M1/1‑43 (MQ=255)
                                                        cAACTGGCTCAACTCTTACGTCTACCAGATCACGCCATACGAcggcgcatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattc                                           >  1:385166‑M1/1‑42 (MQ=255)
                                                        cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtaatagtcgatcgttaagcgattc                                           <  1:1777751‑M1/100‑59 (MQ=255)
                                                         aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattca                                          >  1:1027623‑M1/1‑41 (MQ=255)
                                                         aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcgacacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattca                                          <  2:358783‑M1/100‑60 (MQ=255)
                                                          aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcag                                         >  1:2070693‑M1/1‑40 (MQ=255)
                                                          aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcag                                         <  1:2288054‑M1/100‑61 (MQ=255)
                                                           cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagc                                        <  1:254752‑M1/100‑62 (MQ=255)
                                                           cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagc                                        <  2:1720939‑M1/100‑62 (MQ=255)
                                                             ggCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcac                                      >  2:449706‑M1/1‑37 (MQ=255)
                                                             ggCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcac                                      >  1:1376340‑M1/1‑37 (MQ=255)
                                                              gCTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcacc                                     <  1:2338949‑M1/100‑65 (MQ=255)
                                                               cTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcacct                                    >  1:1858726‑M1/1‑35 (MQ=255)
                                                               cTCAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcacct                                    >  1:2280680‑M1/1‑35 (MQ=255)
                                                                 cAACTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcacctta                                  >  2:458313‑M1/1‑33 (MQ=255)
                                                                  aaCTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttac                                 <  1:1357741‑M1/100‑69 (MQ=255)
                                                                    cTATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacct                               >  2:1654978‑M1/1‑30 (MQ=255)
                                                                     tATTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctc                              <  1:1926889‑M1/100‑72 (MQ=255)
                                                                      aTTACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctca                             <  2:1027676‑M1/100‑73 (MQ=255)
                                                                       ttACGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcag                            <  2:1376341‑M1/100‑74 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggc                          <  2:2114834‑M1/100‑76 (MQ=255)
                                                                         aCGTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggc                          >  1:383659‑M1/1‑25 (MQ=255)
                                                                           gTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcac                        >  1:1934913‑M1/1‑23 (MQ=255)
                                                                           gTCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcac                        >  1:847374‑M1/1‑23 (MQ=255)
                                                                            tCTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcacc                       >  2:172408‑M1/1‑22 (MQ=255)
                                                                             cTACCAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcacct                      >  2:518768‑M1/1‑21 (MQ=255)
                                                                                 cAGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgg                  >  2:1651651‑M1/1‑17 (MQ=255)
                                                                                  aGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcggg                 >  2:1917139‑M1/1‑16 (MQ=255)
                                                                                  aGATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcggg                 >  2:1917371‑M1/1‑16 (MQ=255)
                                                                                   gATCACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggt                <  1:2017576‑M1/100‑86 (MQ=255)
                                                                                      cACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgcc             >  1:140661‑M1/1‑12 (MQ=255)
                                                                                      cACGCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcacccgaaggtgcc             <  1:2345889‑M1/100‑89 (MQ=255)
                                                                                         gCCATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgccttt          >  1:1606712‑M1/1‑9 (MQ=255)
                                                                                           cATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgccttttt        <  1:2019357‑M1/100‑94 (MQ=255)
                                                                                           cATACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgccttttt        <  1:2019195‑M1/100‑94 (MQ=255)
                                                                                            aTACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgcctttttt       <  1:1419453‑M1/100‑95 (MQ=255)
                                                                                             tACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgcctttttta      <  1:1529510‑M1/100‑96 (MQ=255)
                                                                                             tACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgcctttttta      <  2:1303638‑M1/100‑96 (MQ=255)
                                                                                             tACGAcggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgcctttttta      <  2:554811‑M1/100‑96 (MQ=255)
                                                                                               cgacggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgccttttttatt    <  2:1115209‑M1/100‑98 (MQ=255)
                                                                                                 acggcacatgaatccaacgctggattaatcttctgtgatagtcgatcgttaagcgattcagcaccttacctcaggcaccttcgggtgccttttttatttc  <  1:1553306‑M1/100‑100 (MQ=255)
                                                                                                 |                                                                                                   
GATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTA  >  NC_000913/4215081‑4215277

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