Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,456,536 A→G S83S (TCA→TCG)  b1394 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,456,5360AG87.4% 510.4 / 62.0 207S83S (TCA→TCG) b1394hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/16);  new base G (96/85);  total (106/101)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 210 total aligned reads displayed.

CTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGG  >  NC_000913/1456444‑1456628
                                                                                            |                                                                                            
      ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGaa                                                                                        >  2:3748422/1‑93 (MQ=255)
      ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGaa                                                                                        <  2:3241167/93‑1 (MQ=255)
       gTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAAc                                                                                       <  1:1799702/93‑1 (MQ=255)
        tGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACg                                                                                      >  2:7420616/1‑93 (MQ=255)
          ccGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGtt                                                                                    <  1:1729894/93‑1 (MQ=255)
           cGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGttt                                                                                   <  1:2757691/93‑1 (MQ=255)
           cGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGttt                                                                                   <  2:2345816/93‑1 (MQ=255)
            gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTc                                                                                  <  1:5136575/93‑1 (MQ=255)
            gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTc                                                                                  <  2:2809717/93‑1 (MQ=255)
              gcgcGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCta                                                                                <  2:6283365/93‑1 (MQ=255)
                gcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCtcta                                                                              <  2:4736259/93‑1 (MQ=255)
                  ggATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAc                                                                            >  1:3881000/1‑93 (MQ=255)
                  ggATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAc                                                                            >  1:6892432/1‑93 (MQ=255)
                   gATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAcc                                                                           >  2:5624516/1‑93 (MQ=255)
                   gATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAcc                                                                           <  2:4080632/93‑1 (MQ=255)
                     ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCa                                                                         >  2:5761826/1‑93 (MQ=255)
                      tttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCAc                                                                        <  2:1617423/93‑1 (MQ=255)
                      tttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCAc                                                                        <  2:3136051/93‑1 (MQ=255)
                        tgtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTg                                                                      <  1:638796/93‑1 (MQ=255)
                        tgtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTg                                                                      <  1:5466629/93‑1 (MQ=255)
                         gtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTgg                                                                     <  2:894539/93‑1 (MQ=255)
                            cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTAc                                                                  >  1:2624546/1‑93 (MQ=255)
                             tGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACg                                                                 <  1:5566227/93‑1 (MQ=255)
                               gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTc                                                               >  1:7315216/1‑93 (MQ=255)
                                   ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGCACGTTTCTATAACCCACTGGTACGTCGCCt                                                           <  2:4432322/93‑1 (MQ=255)
                                    gATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTg                                                          >  2:3817516/1‑93 (MQ=255)
                                      tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGc                                                        <  2:5705570/93‑1 (MQ=255)
                                       cTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCa                                                       >  1:2482282/1‑93 (MQ=255)
                                        ttAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaa                                                      >  2:3972172/1‑93 (MQ=255)
                                        ttAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaa                                                      >  2:2381531/1‑93 (MQ=255)
                                         tAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaa                                                     >  1:4266226/1‑93 (MQ=255)
                                          aaCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaaa                                                    >  2:1503232/1‑93 (MQ=255)
                                           aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTATGTCGCCTGGCAAAAc                                                   >  2:1230650/1‑93 (MQ=255)
                                           aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                                                   <  1:1312162/93‑1 (MQ=255)
                                           aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                                                   <  2:4671935/93‑1 (MQ=255)
                                            cGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACt                                                  >  1:5094457/1‑93 (MQ=255)
                                              aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                                <  2:1666124/93‑1 (MQ=255)
                                               tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                               <  2:5449336/93‑1 (MQ=255)
                                                  tAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaaa                                            <  1:6080352/93‑1 (MQ=255)
                                                  tAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaaa                                            >  1:2797938/1‑93 (MQ=255)
                                                    aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                          <  1:3530833/93‑1 (MQ=255)
                                                      gTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCg                                        <  1:87213/93‑1 (MQ=255)
                                                       ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                       <  2:4204240/93‑1 (MQ=255)
                                                       ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                       >  2:3547153/1‑93 (MQ=255)
                                                       ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                       >  2:5088545/1‑93 (MQ=255)
                                                       ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                       <  1:4514577/93‑1 (MQ=255)
                                                       ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                       <  1:5428863/93‑1 (MQ=255)
                                                          aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                    <  1:3830394/93‑1 (MQ=255)
                                                          aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                    <  2:3675264/93‑1 (MQ=255)
                                                           tCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGAt                                   <  1:3008143/93‑1 (MQ=255)
                                                            cccACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc                                  <  2:2994359/93‑1 (MQ=255)
                                                             ccACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCt                                 >  1:3339103/1‑93 (MQ=255)
                                                               acaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                               >  2:840750/1‑93 (MQ=255)
                                                               acaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                               <  1:6505397/93‑1 (MQ=255)
                                                               acaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                               <  2:1995101/93‑1 (MQ=255)
                                                               acaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                               <  1:7010891/93‑1 (MQ=255)
                                                               acaGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                               <  1:5755947/93‑1 (MQ=255)
                                                                caGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgtg                              >  2:6102436/1‑93 (MQ=255)
                                                                 aGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                             <  2:647977/93‑1 (MQ=255)
                                                                 aGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                             <  1:1059585/93‑1 (MQ=255)
                                                                 aGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                             >  2:6935467/1‑93 (MQ=255)
                                                                  ggCCCCGCACCGGATTTAGGGAAGGCGGGTGAACGGTTCTATAACCCATTGGTACGTTGCCCGGCAAAACTGGCAAAACCGGGGATCTGTGCa                            >  1:6182291/1‑93 (MQ=255)
                                                                   gCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAg                           >  2:7167420/1‑93 (MQ=255)
                                                                     cccGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTc                         <  2:3840947/93‑1 (MQ=255)
                                                                      ccGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCa                        >  1:6685631/1‑93 (MQ=255)
                                                                      ccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCa                        >  2:7282310/1‑93 (MQ=255)
                                                                       cGCACCTGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCaa                       <  2:6645144/93‑1 (MQ=255)
                                                                        gCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt                      >  2:6087920/1‑93 (MQ=255)
                                                                        gCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt                      <  2:5404927/93‑1 (MQ=255)
                                                                         cACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg                     <  2:4491556/93‑1 (MQ=255)
                                                                         cACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg                     <  1:4683361/93‑1 (MQ=255)
                                                                          aCCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtgg                    <  2:5129381/93‑1 (MQ=255)
                                                                          aCCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtgg                    >  2:3370602/1‑93 (MQ=255)
                                                                           ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                   >  1:1390126/1‑93 (MQ=255)
                                                                           ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                   >  1:1183750/1‑93 (MQ=255)
                                                                           ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                   <  1:6405346/93‑1 (MQ=255)
                                                                           ccGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                   >  1:4372195/1‑93 (MQ=255)
                                                                            cGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTGCGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcg                  <  1:6348631/93‑1 (MQ=255)
                                                                             ggATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt                 <  2:5998714/93‑1 (MQ=255)
                                                                             ggATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt                 <  2:2927731/93‑1 (MQ=255)
                                                                              gATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtg                <  2:4820580/93‑1 (MQ=255)
                                                                                tttAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc              >  1:2208208/1‑93 (MQ=255)
                                                                                    gggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca          >  2:5919217/1‑93 (MQ=255)
                                                                                    gggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca          >  2:2094100/1‑93 (MQ=255)
                                                                                     ggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcag         >  1:4167749/1‑93 (MQ=255)
                                                                                      gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg        <  1:4022057/93‑1 (MQ=255)
                                                                                      gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg        <  1:4487890/93‑1 (MQ=255)
                                                                                      gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg        >  2:4387524/1‑93 (MQ=255)
                                                                                      gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg        <  2:6855650/93‑1 (MQ=255)
                                                                                      gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg        <  2:1824787/93‑1 (MQ=255)
                                                                                       aTGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggg       <  1:5067404/93‑1 (MQ=255)
                                                                                        tGTCGGTTGAACGTTTCTATAACCCACTGGTGCGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggaagggg      <  1:4961803/93‑1 (MQ=255)
                                                                                        tGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagggg      >  2:4334622/1‑93 (MQ=255)
                                                                                        tGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagggg      >  1:657545/1‑93 (MQ=255)
                                                                                        tGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagggg      >  2:2163806/1‑93 (MQ=255)
                                                                                         gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc     >  1:3225401/1‑93 (MQ=255)
                                                                                         gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc     >  1:980245/1‑93 (MQ=255)
                                                                                         gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc     >  1:1611691/1‑93 (MQ=255)
                                                                                         gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc     >  1:3685155/1‑93 (MQ=255)
                                                                                            aGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggcagg  >  2:5262657/1‑93 (MQ=255)
                                                                                            |                                                                                            
CTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGG  >  NC_000913/1456444‑1456628

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