Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,206,609 C→T S315S (AGC→AGT ygjE → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,6090CT84.4% 800.0 / 102.3 302S315S (AGC→AGTygjEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (106/149);  total (130/172)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

Only 100 of 307 total aligned reads displayed.

TGAAAAACGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTCTGGCTGGCCT  >  NC_000913/3206519‑3206700
                                                                                          |                                                                                           
 gAAAAACGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGtgtg                                                                                          <  2:532588/93‑1 (MQ=255)
 gAAAAACGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGtgtg                                                                                          >  2:6169431/1‑93 (MQ=255)
  aaaaaCGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTgtgg                                                                                         <  2:823055/93‑1 (MQ=255)
  aaaaaCGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTgtgg                                                                                         <  2:7240353/93‑1 (MQ=255)
   aaaaCGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTgtggt                                                                                        <  1:5474740/93‑1 (MQ=255)
   aaaaCGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTgtggt                                                                                        <  2:2555816/93‑1 (MQ=255)
    aaaCGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTgtggtg                                                                                       >  1:6152969/1‑93 (MQ=255)
        ggATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCAc                                                                                   <  1:162022/93‑1 (MQ=255)
        ggATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCAc                                                                                   <  1:6069365/93‑1 (MQ=255)
            gCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGAt                                                                               >  1:6936403/1‑93 (MQ=255)
             cTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGAtg                                                                              >  2:591800/1‑93 (MQ=255)
             cTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGAtg                                                                              <  2:1232543/93‑1 (MQ=255)
                gggCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGAtgctg                                                                           <  1:4365491/93‑1 (MQ=255)
                 ggCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGAtgctgc                                                                          >  2:660142/1‑93 (MQ=255)
                   cTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCtt                                                                        <  2:3152380/93‑1 (MQ=255)
                   cTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCtt                                                                        >  2:1539024/1‑93 (MQ=255)
                     gATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCt                                                                      <  1:2427732/93‑1 (MQ=255)
                     gATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCt                                                                      <  2:1125533/93‑1 (MQ=255)
                       tGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTgc                                                                    >  2:3273977/1‑93 (MQ=255)
                       tGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTgc                                                                    <  2:244543/93‑1 (MQ=255)
                          tAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCa                                                                 <  1:6442437/93‑1 (MQ=255)
                          tAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCa                                                                 >  2:278977/1‑93 (MQ=255)
                             gcgcgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTa                                                              <  1:3158006/93‑1 (MQ=255)
                             gcgcgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTa                                                              <  2:6731239/93‑1 (MQ=255)
                             gcgcgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTa                                                              <  2:7095924/93‑1 (MQ=255)
                              gccgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTAt                                                             <  1:1209698/91‑1 (MQ=255)
                              cgcgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTAt                                                             <  2:6223218/93‑1 (MQ=255)
                               gcgcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATc                                                            <  1:3557615/93‑1 (MQ=255)
                                 gcTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAg                                                          >  1:5780874/1‑93 (MQ=255)
                                  cTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGt                                                         <  1:7060591/93‑1 (MQ=255)
                                   tGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGtt                                                        <  1:4679724/93‑1 (MQ=255)
                                    ggTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTg                                                       >  2:1641292/1‑93 (MQ=255)
                                     gTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTgg                                                      <  2:4228192/93‑1 (MQ=255)
                                     gTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTgg                                                      >  1:7323654/1‑93 (MQ=255)
                                     gTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTgg                                                      <  2:4875037/93‑1 (MQ=255)
                                     gTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTgg                                                      <  1:208984/93‑1 (MQ=255)
                                     gTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTgg                                                      <  2:1373960/93‑1 (MQ=255)
                                       gCTGTGGATTTTCGGCGGGGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGga                                                    >  2:44455/1‑93 (MQ=255)
                                        cTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgac                                                   >  2:3389458/1‑93 (MQ=255)
                                        cTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgac                                                   >  2:3483138/1‑93 (MQ=255)
                                        cTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgac                                                   <  1:2621448/93‑1 (MQ=255)
                                          gtgGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgacga                                                 >  2:5027845/1‑93 (MQ=255)
                                          gtgGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgacga                                                 <  1:6856778/93‑1 (MQ=255)
                                          gtgGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCATTGAGGCTGCTTCTGCGCATTATCAGTTGGgacgc                                                 >  1:1852164/1‑92 (MQ=255)
                                           tgGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGgacgac                                                >  2:4013333/1‑93 (MQ=255)
                                              aTTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACAtt                                             <  1:413572/93‑1 (MQ=255)
                                                  tCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCa                                         <  2:896323/93‑1 (MQ=255)
                                                  tCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCa                                         <  1:7482769/93‑1 (MQ=255)
                                                  tCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCa                                         <  1:339247/93‑1 (MQ=255)
                                                   cggcggTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAg                                        <  1:5409671/93‑1 (MQ=255)
                                                   cggcggTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAg                                        <  2:4540739/93‑1 (MQ=255)
                                                    ggcggTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGt                                       >  2:4300809/1‑93 (MQ=255)
                                                    ggcggTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGt                                       >  1:4215662/1‑93 (MQ=255)
                                                      cggTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGtaa                                     <  1:6105423/93‑1 (MQ=255)
                                                       ggTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGtaat                                    >  2:4640554/1‑93 (MQ=255)
                                                        gTGATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGtaata                                   >  2:6304904/1‑93 (MQ=255)
                                                          gATTATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATaaa                                 <  1:6583613/93‑1 (MQ=255)
                                                          gATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATaaa                                 >  1:2486015/1‑93 (MQ=255)
                                                           aTTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAg                                >  1:4903967/1‑93 (MQ=255)
                                                            ttATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgc                               <  2:5148832/93‑1 (MQ=255)
                                                            ttATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgc                               >  2:1066639/1‑93 (MQ=255)
                                                            ttATATCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgc                               >  1:3972759/1‑93 (MQ=255)
                                                             tataTCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcg                              <  1:4464785/93‑1 (MQ=255)
                                                             tataTCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcg                              <  2:7512864/93‑1 (MQ=255)
                                                             tataTCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcg                              >  1:5065496/1‑93 (MQ=255)
                                                             tataTCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcg                              <  1:3710837/93‑1 (MQ=255)
                                                               tatCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcggc                            >  1:2229722/1‑93 (MQ=255)
                                                               tatCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcggc                            >  1:936423/1‑93 (MQ=255)
                                                                atCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAgcggcg                           >  1:5185619/1‑93 (MQ=255)
                                                                 tCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGt                          <  2:6763418/93‑1 (MQ=255)
                                                                 tCGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGt                          <  2:635950/93‑1 (MQ=255)
                                                                  cGATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTg                         <  1:312662/93‑1 (MQ=255)
                                                                   gATGCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTgg                        >  2:3021079/1‑93 (MQ=255)
                                                                      gCCGCGATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAAc                     <  2:6972587/93‑1 (MQ=255)
                                                                       ccGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACg                    <  1:835599/93‑1 (MQ=255)
                                                                         gcgATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGGCATTGTCAGTAATAAAGCGGCGTGGAACGtt                  <  2:2086643/93‑1 (MQ=255)
                                                                         gcgATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGtt                  <  2:4812057/93‑1 (MQ=255)
                                                                         gcgATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGtt                  >  1:3788395/1‑93 (MQ=255)
                                                                          cgATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGttt                 <  1:326554/93‑1 (MQ=255)
                                                                          cgATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGCCAGTAATAAAGCGGCGTGGAACGttt                 <  2:5235963/93‑1 (MQ=255)
                                                                          cgATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGttt                 >  1:191567/1‑93 (MQ=255)
                                                                           gATGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAGCGtttt                >  1:1519246/1‑93 (MQ=255)
                                                                            aTGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttc               <  2:4973640/93‑1 (MQ=255)
                                                                            aTGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttc               >  1:716260/1‑93 (MQ=255)
                                                                            aTGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttc               >  2:1331400/1‑93 (MQ=255)
                                                                             tGGTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttct              <  2:1862389/93‑1 (MQ=255)
                                                                              ggTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttctt             <  1:1029003/93‑1 (MQ=255)
                                                                              ggTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttctt             >  2:221791/1‑93 (MQ=255)
                                                                               gTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGGGTGGAACGTTttcttc            <  1:2802402/93‑1 (MQ=255)
                                                                               gTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttcttc            <  2:3892356/93‑1 (MQ=255)
                                                                               gTCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTttcttc            >  2:6017304/1‑93 (MQ=255)
                                                                                tCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTtcttct           <  1:1394593/93‑1 (MQ=255)
                                                                                tCGGTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTtcttct           <  1:112623/93‑1 (MQ=255)
                                                                                   gTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggc        >  2:5509428/1‑93 (MQ=255)
                                                                                   gTTACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggc        >  2:2644273/1‑93 (MQ=255)
                                                                                    ttACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggct       <  1:4263866/93‑1 (MQ=255)
                                                                                    ttACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggct       <  2:3679746/93‑1 (MQ=255)
                                                                                     tACAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggctg      <  1:1753332/93‑1 (MQ=255)
                                                                                      aCAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTctggctgg     <  1:6751813/93‑1 (MQ=255)
                                                                                       cAGTGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTCtggctggc    >  1:2678473/1‑93 (MQ=255)
                                                                                          |                                                                                           
TGAAAAACGGATGCTGGGGCTGATGGTAGGCGCGCTGGTGCTGTGGATTTTCGGCGGTGATTATATCGATGCCGCGATGGTCGGTTACAGCGTGGTGGCACTGATGCTGCTTCTGCGCATTATCAGTTGGGACGACATTGTCAGTAATAAAGCGGCGTGGAACGTTTTCTTCTGGCTGGCCT  >  NC_000913/3206519‑3206700

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