Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,412,534 T→C E234E (GAA→GAG intR ← Rac prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,412,5340TC100.0% 217.4 / NA 61E234E (GAA→GAGintRRac prophage; integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/28);  total (33/28)

TTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGTTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTCAAAT  >  NC_000913/1412390‑1412681
                                                                                                                                                |                                                                                                                                                   
ttCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCg                                                                                                                                               >  1:206143/1‑151 (MQ=255)
 tCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGt                                                                                                                                              >  2:686067/1‑151 (MQ=255)
   aaGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTg                                                                                                                                             >  2:345053/1‑150 (MQ=255)
     gagCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACg                                                                                                                                          <  2:773175/151‑1 (MQ=255)
      agCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACg                                                                                                                                          >  1:532552/1‑150 (MQ=255)
      agCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACg                                                                                                                                          >  1:666077/1‑150 (MQ=255)
              gCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTg                                                                                                                                  >  1:52014/1‑150 (MQ=255)
                agtagtGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGt                                                                                                                               >  1:799442/1‑151 (MQ=255)
                 gtagtGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGt                                                                                                                               >  2:2824/1‑150 (MQ=255)
                    gtgtTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAAc                                                                                                                           <  1:595127/151‑1 (MQ=255)
                      gtTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACg                                                                                                                          <  2:203329/150‑1 (MQ=255)
                       ttACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGc                                                                                                                        >  2:811696/1‑151 (MQ=255)
                        tACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGc                                                                                                                        <  2:537982/150‑1 (MQ=255)
                            gTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCaaa                                                                                                                     <  2:624484/149‑1 (MQ=255)
                              tCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTc                                                                                                                 >  2:620306/1‑151 (MQ=255)
                               cGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTcc                                                                                                                <  2:181681/151‑1 (MQ=255)
                               cGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTcc                                                                                                                >  2:428652/1‑151 (MQ=255)
                                 gATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCca                                                                                                               <  1:283694/150‑1 (MQ=255)
                                    cccGGCATTTGTTTTTGGGGTACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTACcacac                                                                                                            <  2:62945/150‑2 (MQ=255)
                                       ggCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTc                                                                                                        <  1:646959/151‑1 (MQ=255)
                                         cATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTg                                                                                                      <  2:393929/151‑1 (MQ=255)
                                              gTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGctct                                                                                                 >  2:158673/1‑151 (MQ=255)
                                                 tttGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGt                                                                                               <  2:359938/150‑1 (MQ=255)
                                                      gggACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCtttt                                                                                         <  2:874884/151‑1 (MQ=255)
                                                       ggACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACTAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCttttt                                                                                        <  2:341221/151‑1 (MQ=255)
                                                       ggACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCttttt                                                                                        <  2:664810/151‑1 (MQ=255)
                                                         aCCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTc                                                                                       >  1:144766/1‑150 (MQ=255)
                                                           cGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACt                                                                                    >  1:213513/1‑151 (MQ=255)
                                                              aCATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTc                                                                                 <  2:553687/151‑1 (MQ=255)
                                                                    cAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTa                                                                            >  2:13764/1‑150 (MQ=255)
                                                                      aGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATg                                                                         <  1:162366/151‑1 (MQ=255)
                                                                       gTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGc                                                                        <  1:286341/151‑1 (MQ=255)
                                                                         aTCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCa                                                                       <  1:537947/150‑1 (MQ=255)
                                                                               aGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAAtt                                                                >  2:823044/1‑151 (MQ=255)
                                                                                  tttCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTg                                                             >  1:685472/1‑151 (MQ=255)
                                                                                   ttcttcTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTg                                                             <  2:853931/150‑1 (MQ=255)
                                                                                     cttctGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGttt                                                          >  2:75023/1‑151 (MQ=255)
                                                                                        ctGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGtttt                                                         >  1:551042/1‑149 (MQ=255)
                                                                                         tGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGtttt                                                         >  2:110311/1‑148 (MQ=255)
                                                                                          gACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGtttt                                                         >  2:606791/1‑147 (MQ=255)
                                                                                          gACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTa                                                        <  2:439847/148‑1 (MQ=255)
                                                                                              ttCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACg                                                 >  1:258804/1‑151 (MQ=255)
                                                                                               tCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACg                                                 <  1:823111/150‑1 (MQ=255)
                                                                                                    aTTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATc                                             >  1:152974/1‑149 (MQ=255)
                                                                                                     ttCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAgg                                          >  1:882261/1‑151 (MQ=255)
                                                                                                     ttCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAgg                                          >  1:8454/1‑151 (MQ=255)
                                                                                                            ttCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCCTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTAt                                    <  1:428698/150‑1 (MQ=255)
                                                                                                              cGAGATCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCt                                 <  1:708432/151‑1 (MQ=255)
                                                                                                               gagaTCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTg                                >  1:319147/1‑151 (MQ=255)
                                                                                                                 gaTCCACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGtt                              >  1:693632/1‑151 (MQ=255)
                                                                                                                    ccACATCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCg                           <  1:525318/151‑1 (MQ=255)
                                                                                                                     cacaTCCTCCCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGc                          >  1:173325/1‑151 (MQ=255)
                                                                                                                         tcctccCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGt                       >  2:378951/1‑150 (MQ=255)
                                                                                                                            tccCACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAAtt                   >  2:374162/1‑151 (MQ=255)
                                                                                                                              ccATGCCAGATCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTATGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAAtttt                 >  2:879261/1‑151 (MQ=255)
                                                                                                                              ccACGCCAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAAtttt                 >  1:161088/1‑151 (MQ=255)
                                                                                                                                   ccAGAGCTGCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTAc            >  1:763678/1‑151 (MQ=255)
                                                                                                                                           gCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTc      <  1:304968/149‑1 (MQ=255)
                                                                                                                                           gCCAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTCa     <  1:177048/150‑1 (MQ=255)
                                                                                                                                             cAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTCAAAt  <  2:667420/151‑1 (MQ=255)
                                                                                                                                             cAGCTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTCAAAt  <  1:158696/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCCACGCCAGAGCTGCCAGTTCCCCGTGACGAAGTCCTGAGTAAACGGCAAATTTCCACAAGTTCTGGCTCTGTCCTTTTTCACTTTCCATTAATGCATTGAATTCTGTTTTAGATAACGGATCAGGCTTTATTCTGTTTCGCTGTAATTTTTTTACTCCTTCAAAT  >  NC_000913/1412390‑1412681

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