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.1% 1095.2 / ‑2.2 321D330G (GAC→GGC) topADNA topoisomerase type I omega protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (150/168);  minor base T (3/0);  total (153/168)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 329 total aligned reads displayed.

GCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329970‑1330149
                                                                                          |                                                                                         
  gTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGcgcg                                                                                       <  1:4223860/93‑1 (MQ=255)
   tGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCgg                                                                                      >  1:2033834/1‑93 (MQ=255)
   tGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCgg                                                                                      <  1:3567424/93‑1 (MQ=255)
   ggAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCgg                                                                                      <  1:2126779/92‑1 (MQ=255)
     aaaaaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTa                                                                                    <  2:7633296/93‑1 (MQ=255)
     aaaaaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTa                                                                                    >  2:415139/1‑93 (MQ=255)
       aaaaaCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTaaa                                                                                  <  2:4983284/93‑1 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              <  2:8099117/93‑1 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              <  2:1597246/93‑1 (MQ=255)
           aCCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATg                                                                              <  2:3891352/93‑1 (MQ=255)
            ccATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATgg                                                                             >  2:7394256/1‑93 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  1:2634319/93‑1 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  2:1172102/93‑1 (MQ=255)
              atgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGtt                                                                           <  1:4083121/93‑1 (MQ=255)
               tgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTc                                                                          >  2:4421059/1‑93 (MQ=255)
               tgatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTc                                                                          <  2:5581137/93‑1 (MQ=255)
                gatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcg                                                                         >  2:2338282/1‑93 (MQ=255)
                gatgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTcg                                                                         <  1:6129719/93‑1 (MQ=255)
                 atgatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTTCACTAACCTGAGTCAGGGCGCGGTAAATATAGTTcgc                                                                        >  1:5984299/1‑93 (MQ=255)
                   gatgGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTAAGTCAGGGCGCGGGAAATATGGTTCGCgg                                                                      >  1:4950026/1‑93 (MQ=255)
                      ggCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAGCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTta                                                                   <  1:1044998/93‑1 (MQ=255)
                       gcgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTtat                                                                  <  1:8550074/93‑1 (MQ=255)
                        cgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCTGTAAATATGGTTCGCGGTtata                                                                 <  2:6177994/93‑1 (MQ=255)
                        cgcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTtata                                                                 <  1:322210/93‑1 (MQ=255)
                         gcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTatat                                                                >  2:606092/1‑93 (MQ=255)
                         gcAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTatat                                                                >  1:1672612/1‑93 (MQ=255)
                          cAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATc                                                               <  2:5477401/93‑1 (MQ=255)
                           aGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCa                                                              <  1:7836759/93‑1 (MQ=255)
                            gCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAg                                                             <  1:4319957/93‑1 (MQ=255)
                            gCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAg                                                             >  1:7552269/1‑93 (MQ=255)
                             cGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGc                                                            >  1:7998588/1‑93 (MQ=255)
                              gTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCg                                                           >  2:8083823/1‑93 (MQ=255)
                              gTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCg                                                           >  2:2515332/1‑93 (MQ=255)
                                ttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGAt                                                         <  1:3319502/93‑1 (MQ=255)
                                ttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGAt                                                         <  1:3775448/93‑1 (MQ=255)
                                ttGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGAt                                                         <  1:4011995/93‑1 (MQ=255)
                                 tGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATa                                                        <  1:5590572/93‑1 (MQ=255)
                                 tGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATa                                                        <  2:3457305/93‑1 (MQ=255)
                                  gTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATaa                                                       <  1:5851590/93‑1 (MQ=255)
                                    aTGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtt                                                     <  2:429322/93‑1 (MQ=255)
                                    aTGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtt                                                     <  2:4851415/93‑1 (MQ=255)
                                      gAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtttt                                                   <  1:3183412/93‑1 (MQ=255)
                                      gAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAAtttt                                                   >  2:7027408/1‑93 (MQ=255)
                                       aaGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTg                                                  >  1:2536039/1‑93 (MQ=255)
                                       aaGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTg                                                  >  2:262606/1‑93 (MQ=255)
                                         gCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGt                                                <  1:7179512/93‑1 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGTCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              >  1:2786574/1‑93 (MQ=255)
                                           aGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTaa                                              >  2:6786268/1‑93 (MQ=255)
                                              cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaa                                           <  2:5882360/93‑1 (MQ=255)
                                              cTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaa                                           <  2:4704832/93‑1 (MQ=255)
                                               tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGaaa                                          >  2:598563/1‑93 (MQ=255)
                                               tataTCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGGAAATATGTGTTGCGGATATGCCAGCGATAATTTTGGTAAGata                                          >  1:2896280/1‑91 (MQ=255)
                                                atatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAat                                         >  1:685708/1‑93 (MQ=255)
                                                 tatCACTTACATGCGTACCGACTCCACTAGCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        <  1:7250175/93‑1 (MQ=255)
                                                 tatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        <  2:3387689/93‑1 (MQ=255)
                                                 tatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        <  2:4004887/93‑1 (MQ=255)
                                                 tatCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAata                                        <  2:7534420/93‑1 (MQ=255)
                                                   tCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                      >  1:6659363/1‑93 (MQ=255)
                                                    cACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCt                                     >  1:2390591/1‑93 (MQ=255)
                                                     aCTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTTGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTg                                    <  2:2177503/93‑1 (MQ=255)
                                                     aCTTACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTg                                    >  2:8511587/1‑93 (MQ=255)
                                                       ttACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTCTATCAGCGATAATTTTGGTAAGAAATATCTGcc                                  >  1:842918/1‑93 (MQ=255)
                                                       ttACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCAGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGcc                                  <  2:4257129/93‑1 (MQ=255)
                                                        tACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCg                                 >  1:5748401/1‑93 (MQ=255)
                                                        tACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCg                                 >  2:3745251/1‑93 (MQ=255)
                                                        tACATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCg                                 <  1:5201337/93‑1 (MQ=255)
                                                         aCATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCgg                                >  2:8552228/1‑93 (MQ=255)
                                                          cATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGa                               <  1:6604387/93‑1 (MQ=255)
                                                          cATGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGa                               >  2:5761093/1‑93 (MQ=255)
                                                           aTGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaa                              >  2:3223362/1‑93 (MQ=255)
                                                            tGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             >  2:4194318/1‑93 (MQ=255)
                                                            tGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             <  2:6644183/93‑1 (MQ=255)
                                                            tGCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGaaa                             >  2:5932054/1‑93 (MQ=255)
                                                             gCGTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAg                            <  2:7468295/93‑1 (MQ=255)
                                                               gTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          <  2:3271184/93‑1 (MQ=255)
                                                               gTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          >  1:3420427/1‑93 (MQ=255)
                                                               gTACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          <  1:7675172/93‑1 (MQ=255)
                                                               gTACCGAATCCACTAACCTGAGTCAGGGAGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTc                          <  2:6344701/93‑1 (MQ=255)
                                                                tACCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTcc                         <  2:273125/93‑1 (MQ=255)
                                                                 aCCGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCg                        <  2:3072277/93‑1 (MQ=255)
                                                                  ccGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGa                       >  2:8812181/1‑93 (MQ=255)
                                                                   cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      >  1:1378626/1‑93 (MQ=255)
                                                                   cGACTCCACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                      <  1:8510752/93‑1 (MQ=255)
                                                                        ccACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATAAga                 >  1:7089075/1‑92 (MQ=255)
                                                                         cACTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTa                >  1:4136451/1‑93 (MQ=255)
                                                                          aCTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTAc               >  2:1422005/1‑93 (MQ=255)
                                                                          aCTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTAc               <  2:8182529/93‑1 (MQ=255)
                                                                           cTAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACg              >  1:1818908/1‑93 (MQ=255)
                                                                            tAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGc             <  2:785239/93‑1 (MQ=255)
                                                                            tAACCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGc             <  2:8342604/93‑1 (MQ=255)
                                                                             aaCCTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc            >  1:6348690/1‑93 (MQ=255)
                                                                                cTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGc         >  2:1727762/1‑93 (MQ=255)
                                                                                cTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGc         <  1:8194992/93‑1 (MQ=255)
                                                                                cTGAGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGc         >  1:739578/1‑93 (MQ=255)
                                                                                   aGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCaaa      <  2:6716781/93‑1 (MQ=255)
                                                                                   aGTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCaaa      <  1:4589726/93‑1 (MQ=255)
                                                                                    gTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAg     <  2:8308141/93‑1 (MQ=255)
                                                                                    gTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAg     <  2:713765/93‑1 (MQ=255)
                                                                                    gTCAGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAg     <  2:2000212/93‑1 (MQ=255)
                                                                                       aGGGCGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGaaa  >  1:1485277/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAA  >  NC_000913/1329970‑1330149

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