Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,995,552 C→A A264E (GCG→GAG)  dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,995,5520CA100.0% 232.0 / NA 66A264E (GCG→GAG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (44/22);  total (44/22)

GCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCG  >  NC_000913/3995417‑3995687
                                                                                                                                       |                                                                                                                                       
gcggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTa                                                                                                                                     <  2:1291555/140‑1 (MQ=255)
 cggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTAc                                                                                                                                    <  1:1395291/140‑1 (MQ=255)
 cggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTAc                                                                                                                                    <  1:332183/140‑1 (MQ=255)
           gcggcgGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTacga                                                                                                                          <  2:296818/140‑1 (MQ=255)
                               tCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGaa                                                                                                      >  2:1406951/1‑140 (MQ=255)
                               tCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGaa                                                                                                      >  2:1284754/1‑140 (MQ=255)
                                cAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAg                                                                                                     <  1:249738/140‑1 (MQ=255)
                                  gCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCa                                                                                                   >  1:722477/1‑140 (MQ=255)
                                     aGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAcc                                                                                                >  2:160620/1‑140 (MQ=255)
                                             tGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAgggg                                                                                           <  1:609760/137‑1 (MQ=255)
                                             tGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAgggg                                                                                           >  2:609760/1‑137 (MQ=255)
                                                       agTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGa                                                                              <  2:564224/140‑1 (MQ=255)
                                                           cgcgTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAaca                                                                          >  1:341192/1‑140 (MQ=255)
                                                             cgTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACt                                                                        >  2:41351/1‑140 (MQ=255)
                                                               tGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCa                                                                      <  1:1284754/140‑1 (MQ=255)
                                                               tGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCa                                                                      <  2:722477/140‑1 (MQ=255)
                                                               tGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCa                                                                      <  2:216532/140‑1 (MQ=255)
                                                                  aaCTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACgg                                                                   <  1:249563/140‑1 (MQ=255)
                                                                  aaCTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACgg                                                                   >  2:249563/1‑140 (MQ=255)
                                                                   aCTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGa                                                                  >  1:1302107/1‑140 (MQ=255)
                                                                       cccGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCtt                                                              >  1:975477/1‑140 (MQ=255)
                                                                                 gtcTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCggg                                                    >  2:443862/1‑140 (MQ=255)
                                                                                 gtcTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCggg                                                    <  2:341192/140‑1 (MQ=255)
                                                                                     tGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGt                                                <  1:1406951/140‑1 (MQ=255)
                                                                                                 gAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACTggg                                                    >  2:1141299/1‑124 (MQ=255)
                                                                                                 gAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACTggg                                                    <  1:1141299/124‑1 (MQ=255)
                                                                                                 gAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCt                                    >  2:112654/1‑140 (MQ=255)
                                                                                                           ggTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTg                                   <  1:578402/131‑1 (MQ=255)
                                                                                                           ggTCACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTg                                   >  2:578402/1‑131 (MQ=255)
                                                                                                              cACCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCt                       <  1:443862/140‑1 (MQ=255)
                                                                                                               aCCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTg                      <  1:1324319/140‑1 (MQ=255)
                                                                                                               aCCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTg                      <  1:160620/140‑1 (MQ=255)
                                                                                                               aCCCGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTg                      <  1:1184160/140‑1 (MQ=255)
                                                                                                                cccGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGa                     >  2:1237139/1‑140 (MQ=255)
                                                                                                                 ccGTTATATATGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAg                    >  1:269339/1‑140 (MQ=255)
                                                                                                                           tGACTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATcc          >  2:154972/1‑140 (MQ=255)
                                                                                                                             aCTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGt        >  1:201766/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  1:1044984/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  2:1404474/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  2:448439/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  1:447699/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  2:775232/1‑140 (MQ=255)
                                                                                                                              cTGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTa       >  1:234427/1‑140 (MQ=255)
                                                                                                                               tGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAGCACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTaa      >  2:1251060/1‑140 (MQ=255)
                                                                                                                               tGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTaa      >  2:1165024/1‑140 (MQ=255)
                                                                                                                               tGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTaa      >  1:1059154/1‑140 (MQ=255)
                                                                                                                               tGGCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTaa      >  2:287157/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGttttt               >  2:309751/1‑129 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGttttt               <  1:309751/129‑1 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:885683/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:628211/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:1007446/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:303614/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:1337365/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:1317949/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:1292031/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  2:1124138/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:978409/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    <  1:775232/140‑1 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:621752/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:52940/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:410991/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:272545/1‑140 (MQ=255)
                                                                                                                                 gCCCGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATg    >  1:241804/1‑140 (MQ=255)
                                                                                                                                  cccGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATgc   >  2:621787/1‑140 (MQ=255)
                                                                                                                                   ccGGAGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATgcg  <  2:1007446/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCG  >  NC_000913/3995417‑3995687

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