Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,330,060 A→G D330G (GAC→GGC)  topA → DNA topoisomerase type I omega protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,330,0600AG91.7% 748.5 / 47.0 266D330G (GAC→GGC) topADNA topoisomerase type I omega protein
Reads supporting (aligned to +/- strand):  ref base A (11/11);  new base G (125/119);  total (136/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 268 total aligned reads displayed.

TGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329968‑1330149
                                                                                            |                                                                                         
   cGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGgcgc                                                                                        <  2:1486868/93‑1 (MQ=255)
    gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       >  2:2177991/1‑93 (MQ=255)
    gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       <  2:1904619/93‑1 (MQ=255)
    gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       >  1:7158793/1‑93 (MQ=255)
         aaaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTaaa                                                                                  <  1:2979237/93‑1 (MQ=255)
          aaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAat                                                                                 <  2:2122909/93‑1 (MQ=255)
          aaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAat                                                                                 <  2:5025488/93‑1 (MQ=255)
            aaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAatat                                                                               <  1:2831110/93‑1 (MQ=255)
             aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              <  2:7529150/93‑1 (MQ=255)
             aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              >  1:1763821/1‑93 (MQ=255)
              ccATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATgg                                                                             >  1:590643/1‑93 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           >  1:2674676/1‑93 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  2:2627549/93‑1 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  2:1602861/93‑1 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  2:6102871/93‑1 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           >  2:623434/1‑93 (MQ=255)
                atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt                                                                           <  2:4308632/93‑1 (MQ=255)
                 tgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTc                                                                          >  2:4611329/1‑93 (MQ=255)
                 tgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTc                                                                          >  2:1722413/1‑93 (MQ=255)
                 tgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTc                                                                          >  2:3831760/1‑93 (MQ=255)
                  gatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcg                                                                         >  2:4380066/1‑93 (MQ=255)
                   atgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcgc                                                                        >  1:5251462/1‑93 (MQ=255)
                    tgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcgcg                                                                       <  1:3529614/93‑1 (MQ=255)
                     gatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCgg                                                                      >  2:7093978/1‑93 (MQ=255)
                     gatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCgg                                                                      >  1:4821001/1‑93 (MQ=255)
                      atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGt                                                                     <  2:4007172/93‑1 (MQ=255)
                      atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGt                                                                     >  1:7437727/1‑93 (MQ=255)
                      atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGCACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGt                                                                     <  1:2260386/93‑1 (MQ=255)
                          cgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTtata                                                                 <  2:6475279/93‑1 (MQ=255)
                           gcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTatat                                                                <  1:7372021/93‑1 (MQ=255)
                            cAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATc                                                               <  1:5106440/93‑1 (MQ=255)
                            cAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATc                                                               >  1:2109265/1‑93 (MQ=255)
                             aGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCa                                                              <  2:124613/93‑1 (MQ=255)
                              gCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAg                                                             >  2:4205612/1‑93 (MQ=255)
                               cGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGc                                                            <  2:5999825/93‑1 (MQ=255)
                                gTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCg                                                           >  2:7012139/1‑93 (MQ=255)
                                 tttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGa                                                          <  2:7453194/93‑1 (MQ=255)
                                 tttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGa                                                          >  1:3924799/1‑93 (MQ=255)
                                 tttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGa                                                          >  1:4386016/1‑93 (MQ=255)
                                   tGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATa                                                        <  2:6793725/93‑1 (MQ=255)
                                    gTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATaa                                                       >  2:4042884/1‑93 (MQ=255)
                                     tATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAt                                                      <  1:3595181/93‑1 (MQ=255)
                                      aTGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtt                                                     <  1:6111646/93‑1 (MQ=255)
                                      aTGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtt                                                     >  2:2226623/1‑93 (MQ=255)
                                       tGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAttt                                                    <  2:6409182/93‑1 (MQ=255)
                                         aaGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTg                                                  <  1:3143026/93‑1 (MQ=255)
                                          aGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTgg                                                 >  1:5614897/1‑93 (MQ=255)
                                           gCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGt                                                >  1:970553/1‑93 (MQ=255)
                                           gCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGt                                                >  2:4193577/1‑93 (MQ=255)
                                           gCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGt                                                >  1:3421355/1‑93 (MQ=255)
                                             nggCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              <  2:1334948/92‑1 (MQ=255)
                                               gCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGa                                            <  1:6435213/93‑1 (MQ=255)
                                                cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGgaagaa                                           >  1:6173880/1‑93 (MQ=255)
                                                cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGag                                           <  1:5847530/93‑2 (MQ=255)
                                                cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaa                                           >  2:5525234/1‑93 (MQ=255)
                                                 tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          <  2:1809889/93‑1 (MQ=255)
                                                 tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          >  1:4707660/1‑93 (MQ=255)
                                                 tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          <  2:3739084/93‑1 (MQ=255)
                                                  atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         >  1:5525296/1‑93 (MQ=255)
                                                  atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         >  1:2855315/1‑93 (MQ=255)
                                                  atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         >  1:1379406/1‑93 (MQ=255)
                                                   tatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        <  1:5329658/93‑1 (MQ=255)
                                                    atCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAatat                                       >  2:2446577/1‑93 (MQ=255)
                                                     tCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                      <  2:1469836/93‑1 (MQ=255)
                                                     tCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                      <  2:6273601/93‑1 (MQ=255)
                                                     tCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                      <  1:3759841/93‑1 (MQ=255)
                                                     tCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                      >  1:6253779/1‑93 (MQ=255)
                                                      cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCt                                     <  1:3291884/93‑1 (MQ=255)
                                                      cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCt                                     <  1:5867470/93‑1 (MQ=255)
                                                      cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCt                                     <  1:6856640/93‑1 (MQ=255)
                                                      cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAACGATAATTTTGGTAAGAAATATCt                                     >  2:376996/1‑93 (MQ=255)
                                                       aCTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTg                                    <  2:2757773/93‑1 (MQ=255)
                                                        cTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGc                                   >  2:2640701/1‑93 (MQ=255)
                                                          tACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCg                                 >  2:5260389/1‑93 (MQ=255)
                                                             aTGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaa                              >  2:5105180/1‑93 (MQ=255)
                                                             aTGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaa                              <  1:1270777/93‑1 (MQ=255)
                                                             aTGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaa                              <  1:4792673/93‑1 (MQ=255)
                                                              tGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             <  2:5447288/93‑1 (MQ=255)
                                                               gCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAg                            >  1:1023071/1‑93 (MQ=255)
                                                                  tACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTcc                         >  1:3901838/1‑93 (MQ=255)
                                                                  tACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTcc                         >  2:4990975/1‑93 (MQ=255)
                                                                     cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      <  2:2489775/93‑1 (MQ=255)
                                                                       aCTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATc                    <  2:4932610/93‑1 (MQ=255)
                                                                       aCTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATc                    <  1:4570990/93‑1 (MQ=255)
                                                                       aCTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATc                    >  2:5931799/1‑93 (MQ=255)
                                                                       aCTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGCAAGAAATATCTGCCGGAAAGGCCGAATc                    >  1:2573588/1‑93 (MQ=255)
                                                                         tCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAg                  <  1:1569167/93‑1 (MQ=255)
                                                                            aCTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTAc               >  2:5769778/1‑93 (MQ=255)
                                                                             cTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACg              >  1:6610157/1‑93 (MQ=255)
                                                                              tAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGc             >  1:868575/1‑93 (MQ=255)
                                                                              tAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGc             <  2:3952372/93‑1 (MQ=255)
                                                                              tAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGc             <  2:337189/93‑1 (MQ=255)
                                                                               aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            >  2:2544963/1‑93 (MQ=255)
                                                                                aCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCa           >  1:594673/1‑93 (MQ=255)
                                                                                 ccTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAg          <  2:3159953/93‑1 (MQ=255)
                                                                                  cTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGc         <  1:4320138/93‑1 (MQ=255)
                                                                                   tGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCa        <  2:6660760/93‑1 (MQ=255)
                                                                                      gTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAg     >  1:4049910/1‑93 (MQ=255)
                                                                                       tCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGa    >  1:5925118/1‑93 (MQ=255)
                                                                                         aGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGaaa  >  2:4950792/1‑93 (MQ=255)
                                                                                            |                                                                                         
TGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329968‑1330149

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