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,0600AG99.7% 1013.9 / ‑4.7 293D330G (GAC→GGC) topADNA topoisomerase type I omega protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (141/151);  total (142/151)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 299 total aligned reads displayed.

GCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329970‑1330149
                                                                                          |                                                                                         
 cGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGgcgc                                                                                        <  1:3738264/93‑1 (MQ=255)
  gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       >  1:1535337/1‑93 (MQ=255)
  gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       >  1:4397544/1‑93 (MQ=255)
   tGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCgg                                                                                      <  2:71533/93‑1 (MQ=255)
   tGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCgg                                                                                      <  1:6012215/93‑1 (MQ=255)
    gAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGt                                                                                     >  2:149535/1‑93 (MQ=255)
    gAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGt                                                                                     >  2:4257819/1‑93 (MQ=255)
       aaaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTaaa                                                                                  <  2:1378694/93‑1 (MQ=255)
        aaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAat                                                                                 <  2:5392754/93‑1 (MQ=255)
        aaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAat                                                                                 >  2:4650883/1‑93 (MQ=255)
        aaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAat                                                                                 <  1:3517125/93‑1 (MQ=255)
          aaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAatat                                                                               >  2:2305415/1‑93 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              >  1:6177627/1‑93 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              >  2:1889722/1‑93 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              >  1:4016776/1‑93 (MQ=255)
            ccATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATgg                                                                             >  1:4204894/1‑93 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAGTATGGtt                                                                           >  2:6279489/1‑93 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  1:891217/93‑1 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           >  1:4754283/1‑93 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           >  1:695895/1‑93 (MQ=255)
                gatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcg                                                                         >  2:5490138/1‑93 (MQ=255)
                   gatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCgg                                                                      >  1:525207/1‑93 (MQ=255)
                    atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGt                                                                     >  2:2644976/1‑93 (MQ=255)
                    atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGt                                                                     <  1:1387497/93‑1 (MQ=255)
                    atgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGt                                                                     <  2:4009814/93‑1 (MQ=255)
                     tgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGtt                                                                    <  1:5839811/93‑1 (MQ=255)
                     tgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGtt                                                                    <  2:5945909/93‑1 (MQ=255)
                       gcgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTtat                                                                  <  1:2006167/93‑1 (MQ=255)
                        cgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTtata                                                                 <  1:3034680/93‑1 (MQ=255)
                         gcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTatat                                                                >  2:3576106/1‑93 (MQ=255)
                           aGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCa                                                              >  2:75906/1‑93 (MQ=255)
                            gCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAg                                                             >  2:4292914/1‑93 (MQ=255)
                             cGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGc                                                            <  2:1254917/93‑1 (MQ=255)
                              gTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCg                                                           <  1:5622543/93‑1 (MQ=255)
                                ttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGAt                                                         >  1:2221428/1‑93 (MQ=255)
                                  gTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATCTCAGCGATaa                                                       >  2:960675/1‑93 (MQ=255)
                                  gTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATaa                                                       >  1:2404692/1‑93 (MQ=255)
                                  gTATAAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATaa                                                       >  1:1142210/1‑93 (MQ=255)
                                   tATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAt                                                      <  2:2481380/93‑1 (MQ=255)
                                   tATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAt                                                      <  2:5258733/93‑1 (MQ=255)
                                    aTGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtt                                                     >  2:226888/1‑93 (MQ=255)
                                     tGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAttt                                                    <  1:5491275/93‑1 (MQ=255)
                                      gAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtttt                                                   >  2:4070753/1‑93 (MQ=255)
                                      gAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtttt                                                   <  2:2975579/93‑1 (MQ=255)
                                       aaGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTg                                                  >  1:4519543/1‑93 (MQ=255)
                                        aGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTgg                                                 >  1:4059305/1‑93 (MQ=255)
                                         gCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTGAATATGGTTCGCGGTTATATCAGCGATAATTTTGGt                                                >  2:1693502/1‑93 (MQ=255)
                                          cAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTa                                               <  1:2104289/93‑1 (MQ=255)
                                          cAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTa                                               <  1:2993062/93‑1 (MQ=255)
                                          cAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTa                                               <  2:2064284/93‑1 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              <  2:3007070/93‑1 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              <  1:3216075/93‑1 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              <  1:3626540/93‑1 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              >  1:1882838/1‑93 (MQ=255)
                                            ggCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAg                                             <  1:4810732/93‑1 (MQ=255)
                                            ggCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAg                                             <  2:1068977/93‑1 (MQ=255)
                                              cTATATCACTTACATGCGTACCGACTCCACTAACCTGCGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaa                                           >  2:2761882/1‑93 (MQ=255)
                                              cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaa                                           >  2:3453855/1‑93 (MQ=255)
                                               tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          >  1:5377282/1‑93 (MQ=255)
                                               tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          >  2:4571185/1‑93 (MQ=255)
                                                atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         >  1:2736150/1‑93 (MQ=255)
                                                atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         <  2:6201611/93‑1 (MQ=255)
                                                 tatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        >  1:2271932/1‑93 (MQ=255)
                                                  atCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAatat                                       >  2:517767/1‑93 (MQ=255)
                                                  atCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAatat                                       >  2:3528848/1‑93 (MQ=255)
                                                    cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCt                                     <  2:5713681/93‑1 (MQ=255)
                                                     aCTTACATGCGTCCCGACTCCACTAACCCGAGTCAGGGCGCGGTAAATATGGTTCGCGGCTATATCAGCGATAATTTTGGTAAGAAATATCTg                                    >  1:934351/1‑93 (MQ=255)
                                                      cTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGc                                   >  1:1241277/1‑93 (MQ=255)
                                                      cTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGc                                   <  2:4805160/93‑1 (MQ=255)
                                                        tACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCg                                 <  2:4152896/93‑1 (MQ=255)
                                                           aTGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaa                              >  2:4982566/1‑93 (MQ=255)
                                                            tGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             <  1:2548787/93‑1 (MQ=255)
                                                            tGCGTACCGACTCCACCAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             <  2:5971874/93‑1 (MQ=255)
                                                              cGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGt                           >  2:1999311/1‑93 (MQ=255)
                                                               gTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          >  1:1033238/1‑93 (MQ=255)
                                                               gTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          <  1:2755755/93‑1 (MQ=255)
                                                                tACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTcc                         >  2:4607713/1‑93 (MQ=255)
                                                                  ccGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGa                       >  1:631354/1‑93 (MQ=255)
                                                                  ccGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGa                       <  2:4962751/93‑1 (MQ=255)
                                                                   cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      <  2:823193/93‑1 (MQ=255)
                                                                   cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      <  1:430508/93‑1 (MQ=255)
                                                                   cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      <  1:4933251/93‑1 (MQ=255)
                                                                    gACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAAt                     >  1:554091/1‑93 (MQ=255)
                                                                    gACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAAt                     >  2:2591027/1‑93 (MQ=255)
                                                                    gACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAAt                     >  2:6069784/1‑93 (MQ=255)
                                                                     aCTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATc                    <  2:1208177/93‑1 (MQ=255)
                                                                      cTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCa                   <  2:3272317/93‑1 (MQ=255)
                                                                        ccACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGt                 <  2:3694320/93‑1 (MQ=255)
                                                                         cACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTa                <  2:5839157/93‑1 (MQ=255)
                                                                             aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            >  2:3913843/1‑93 (MQ=255)
                                                                             aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            <  1:257200/93‑1 (MQ=255)
                                                                             aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            <  2:1139826/93‑1 (MQ=255)
                                                                             aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            >  2:2719507/1‑93 (MQ=255)
                                                                                 tGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCa        >  1:3368166/1‑93 (MQ=255)
                                                                                 tGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCa        >  2:5099329/1‑93 (MQ=255)
                                                                                    gTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAg     >  1:818017/1‑93 (MQ=255)
                                                                                     tCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGa    >  2:2906269/1‑93 (MQ=255)
                                                                                     tCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGa    <  1:4641886/93‑1 (MQ=255)
                                                                                      cAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGaa   <  1:64193/93‑1 (MQ=255)
                                                                                       aGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGaaa  >  1:5148902/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329970‑1330149

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