Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,416,377 A→G V337A (GTC→GCC)  recE ← Rac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,416,3770AG100.0% 238.7 / NA 65V337A (GTC→GCC) recERac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/26);  total (39/26)

TTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCAA  >  NC_000913/1416240‑1416512
                                                                                                                                         |                                                                                                                                       
ttCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGcc                                                                                                                                       >  2:726764/1‑140 (MQ=255)
 tCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCt                                                                                                                                      <  2:776241/140‑1 (MQ=255)
 tCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCt                                                                                                                                      >  1:862339/1‑140 (MQ=255)
   cttcttCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCg                                                                                                                                    <  1:2757/140‑1 (MQ=255)
        tctGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGcc                                                                                                                                       >  2:1115517/1‑132 (MQ=255)
                 cGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGc                                                                                                                      >  2:169119/1‑140 (MQ=255)
                    tGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGcttc                                                                                                                   >  1:769447/1‑140 (MQ=255)
                    tGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGcttc                                                                                                                   >  2:1062734/1‑140 (MQ=255)
                    tGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGcttc                                                                                                                   >  1:182107/1‑140 (MQ=255)
                      cGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGcttctt                                                                                                                 >  2:618305/1‑140 (MQ=255)
                      cGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGcttctt                                                                                                                 >  1:922652/1‑140 (MQ=255)
                             ggggCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGa                                                                                                          <  1:479858/140‑1 (MQ=255)
                             ggggCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGa                                                                                                          <  2:539088/140‑1 (MQ=255)
                              gggCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGAc                                                                                                         >  2:990292/1‑140 (MQ=255)
                                gCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGc                                                                                                       <  1:990292/140‑1 (MQ=255)
                                              gCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGgcgc                                                                                         <  2:439553/140‑1 (MQ=255)
                                                aGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCgg                                                                                       >  2:1007306/1‑140 (MQ=255)
                                                 ggCAATATCCACGATTTCCGGATTAGGGTTGACATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCggg                                                                                      >  1:798704/1‑140 (MQ=255)
                                                   cAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGaa                                                                                    <  2:221574/140‑1 (MQ=255)
                                                      tatCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGaataa                                                                                 >  2:718457/1‑140 (MQ=255)
                                                      tatCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGa                                                                                     <  2:1127275/136‑1 (MQ=255)
                                                      tatCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGa                                                                                     >  1:1127275/1‑136 (MQ=255)
                                                          cACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCa                                                                             >  1:929170/1‑140 (MQ=255)
                                                           aCGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAg                                                                            <  1:1007306/140‑1 (MQ=255)
                                                           aCGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAg                                                                            <  1:726764/140‑1 (MQ=255)
                                                           aCGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAg                                                                            <  2:922652/140‑1 (MQ=255)
                                                           aCGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAg                                                                            <  1:742217/140‑1 (MQ=255)
                                                                ttCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGc                                                                       >  2:222338/1‑140 (MQ=255)
                                                                 tCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGcc                                                                      <  1:718457/140‑1 (MQ=255)
                                                                     gATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGc                                                                  <  2:828867/140‑1 (MQ=255)
                                                                      aTCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCa                                                                 <  1:222338/140‑1 (MQ=255)
                                                                                gCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGaa                                                       >  2:316163/1‑140 (MQ=255)
                                                                                    gATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTg                                                   <  1:169119/140‑1 (MQ=255)
                                                                                     aTCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAAt                                                      <  1:376994/136‑1 (MQ=255)
                                                                                     aTCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAAt                                                      >  2:376994/1‑136 (MQ=255)
                                                                                     aTCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGt                                                  >  1:1159413/1‑140 (MQ=255)
                                                                                     aTCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGt                                                  >  2:749611/1‑140 (MQ=255)
                                                                                         gTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTcgcg                                              <  2:798704/140‑1 (MQ=255)
                                                                                         gTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTcgcg                                              <  2:530495/140‑1 (MQ=255)
                                                                                          tttCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCgg                                             >  1:837513/1‑140 (MQ=255)
                                                                                           ttCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGa                                            >  1:552300/1‑140 (MQ=255)
                                                                                                tcagtACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACa                                       >  1:1372271/1‑140 (MQ=255)
                                                                                                    tACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGaaa                                   >  2:482171/1‑140 (MQ=255)
                                                                                                         tGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGc                              <  2:837513/140‑1 (MQ=255)
                                                                                                         tGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGc                              >  1:1109716/1‑140 (MQ=255)
                                                                                                         tGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGc                              <  2:358969/140‑1 (MQ=255)
                                                                                                         tGTTCAGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGc                              >  1:592293/1‑140 (MQ=255)
                                                                                                              aGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGaa                                    >  1:752731/1‑129 (MQ=255)
                                                                                                              aGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGttt                         >  1:299487/1‑140 (MQ=255)
                                                                                                              aGATATTCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTCACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGaa                                    <  2:752731/129‑1 (MQ=255)
                                                                                                                   ttCAGTGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTAtttt                    >  2:1131579/1‑140 (MQ=255)
                                                                                                                        tGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCg               <  2:1372271/140‑1 (MQ=255)
                                                                                                                        tGACGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCg               >  1:823902/1‑140 (MQ=255)
                                                                                                                          aCGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGAt             >  1:697187/1‑140 (MQ=255)
                                                                                                                          aCGTGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGAt             >  2:780998/1‑140 (MQ=255)
                                                                                                                             tGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAt          >  1:1139496/1‑140 (MQ=255)
                                                                                                                             tGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAt          >  1:987362/1‑140 (MQ=255)
                                                                                                                             tGCGCGGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAt          >  1:250010/1‑140 (MQ=255)
                                                                                                                               cgcgGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAttt        >  1:639863/1‑140 (MQ=255)
                                                                                                                               cgcgGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAttt        >  1:179553/1‑140 (MQ=255)
                                                                                                                               cgcgGGGATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAAttt        >  1:729004/1‑140 (MQ=255)
                                                                                                                                    ggATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCa   >  2:853627/1‑140 (MQ=255)
                                                                                                                                     gATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCaa  <  2:592293/140‑1 (MQ=255)
                                                                                                                                     gATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCaa  <  2:1109716/140‑1 (MQ=255)
                                                                                                                                     gATGGCCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCaa  <  1:1131579/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TTCCTTCTTCTGTTACTCGCTGCGGCATCGGGGCAGAGGAGCGACCGCAGGCAATATCCACGATTTCCGGATCAGGGTTGGCATGATCGGTTTCAGTCAGTACTTTGTTCAGATATTCAGTGACGTGCGCGGGGATGACCTCGATCCCAATTGGTGCTTCTTTTACGGACGCAACCACGATGGCGCGGGAATAATCCAGCCCGCCAGGCATGGTGATGAATTTGTCGCGGAAAACAGAAAAGGGCGGTTTATTTTCAGCGATAATTTCCTCAA  >  NC_000913/1416240‑1416512

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