Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,329,911 T→C P280P (CCT→CCC 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,329,9110TC99.1% 323.7 / ‑4.6 106P280P (CCT→CCCtopADNA topoisomerase type I omega protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (48/57);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

Only 100 of 109 total aligned reads displayed.

CAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGC  >  NC_000913/1329821‑1329999
                                                                                          |                                                                                        
caaAGAACAAACTCAGGCTGCGGTAAGTCTGCTGAAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCgg                                                                                        <  1:4264360/93‑1 (MQ=255)
 aaaGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCTCGTTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGgc                                                                                       <  1:816842/93‑1 (MQ=255)
 aaaGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTTCAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGgc                                                                                       <  1:3414672/93‑1 (MQ=255)
  aaGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGgcg                                                                                      >  2:699607/1‑93 (MQ=255)
   aGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGgcgc                                                                                     >  1:2462461/1‑93 (MQ=255)
     aaCAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTc                                                                                   >  2:2865830/1‑93 (MQ=255)
      aCAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTcc                                                                                  <  1:623648/93‑1 (MQ=255)
      aCAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTcc                                                                                  >  2:2285217/1‑93 (MQ=255)
       cAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCt                                                                                 <  1:1156801/93‑1 (MQ=255)
       cAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCt                                                                                 <  1:4122931/93‑1 (MQ=255)
        aaaCTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCtt                                                                                <  2:2840528/93‑1 (MQ=255)
        aaaCTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCtt                                                                                >  2:3625917/1‑93 (MQ=255)
         aaCTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCttt                                                                               <  2:2665999/93‑1 (MQ=255)
         aaCTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCttt                                                                               <  1:2691057/93‑1 (MQ=255)
         aaCTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCttt                                                                               >  2:4356098/1‑93 (MQ=255)
             cAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTttatt                                                                           <  2:2117340/93‑1 (MQ=255)
               ggCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTAc                                                                         >  2:1210316/1‑93 (MQ=255)
                gCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTAcc                                                                        <  1:2541274/93‑1 (MQ=255)
                gCTGCGGTAAGCCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTAcc                                                                        <  2:573829/93‑1 (MQ=255)
                 cTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTTCTTGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACct                                                                       <  2:2110984/93‑1 (MQ=255)
                 cTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACct                                                                       >  1:262889/1‑93 (MQ=255)
                  tGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACctc                                                                      <  1:4452392/93‑1 (MQ=255)
                   gCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACctct                                                                     >  2:3208891/1‑93 (MQ=255)
                     gttAAGTCTGCTGGAAAAAGCGCGCTACAGCGTTCTGTAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTAc                                                                   <  1:170856/91‑1 (MQ=255)
                     ggTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTAc                                                                   >  2:204787/1‑93 (MQ=255)
                     ggTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTAc                                                                   <  2:4668247/93‑1 (MQ=255)
                      gTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACg                                                                  >  2:1620548/1‑93 (MQ=255)
                      gTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACg                                                                  >  1:2875719/1‑93 (MQ=255)
                      gTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACg                                                                  >  2:3536316/1‑93 (MQ=255)
                      gTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACg                                                                  >  1:365762/1‑93 (MQ=255)
                       tAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGc                                                                 <  2:1137876/93‑1 (MQ=255)
                        aaGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCt                                                                >  1:2407993/1‑93 (MQ=255)
                              tctgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGcaaca                                                          <  2:3106361/92‑1 (MQ=255)
                              gctgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGcaaca                                                          >  2:3436197/1‑93 (MQ=255)
                              gctgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGcaaca                                                          >  2:4271630/1‑93 (MQ=255)
                               ctgGAAAAAGCGCGCTACAGCGTGCTTGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGcaacaa                                                         <  1:1852052/93‑1 (MQ=255)
                               ctgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGcaacaa                                                         <  2:2521227/93‑1 (MQ=255)
                                tgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAg                                                        >  1:4193504/1‑93 (MQ=255)
                                 ggAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGc                                                       <  1:84563/93‑1 (MQ=255)
                                    aaaaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAGCCGACAACCAGTAACCCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGc                                                    <  1:208046/93‑1 (MQ=255)
                                    aaaaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGc                                                    <  1:2072039/93‑1 (MQ=255)
                                    aaaaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGc                                                    <  1:2681261/93‑1 (MQ=255)
                                     aaaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGcc                                                   >  2:1486078/1‑93 (MQ=255)
                                     aaaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGcc                                                   <  2:4016826/93‑1 (MQ=255)
                                      aaGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCa                                                  <  1:1691112/93‑1 (MQ=255)
                                       aGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAg                                                 >  1:193930/1‑93 (MQ=255)
                                        gcgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGc                                                >  1:1011537/1‑93 (MQ=255)
                                        gcgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCAAGc                                                >  2:2606358/1‑93 (MQ=255)
                                         cgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCa                                               <  1:4144139/93‑1 (MQ=255)
                                          gcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCAc                                              >  2:1021977/1‑93 (MQ=255)
                                            gcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCAccc                                            <  1:515906/93‑1 (MQ=255)
                                             cTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCg                                           >  2:1097326/1‑93 (MQ=255)
                                              tACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGt                                          >  1:1781010/1‑93 (MQ=255)
                                              tACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGt                                          >  2:148654/1‑93 (MQ=255)
                                               aCAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTc                                         >  1:767754/1‑93 (MQ=255)
                                                cAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCt                                        <  2:4082659/93‑1 (MQ=255)
                                                 aGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCtt                                       <  1:3155408/93‑1 (MQ=255)
                                                 aGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCtt                                       <  2:256023/93‑1 (MQ=255)
                                                  gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTg                                      >  2:2238035/1‑93 (MQ=255)
                                                  gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTg                                      <  2:2473836/93‑1 (MQ=255)
                                                  gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTg                                      >  1:2463904/1‑93 (MQ=255)
                                                  gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTg                                      >  2:4035904/1‑93 (MQ=255)
                                                  gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTg                                      <  1:3796601/93‑1 (MQ=255)
                                                    gTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGa                                    <  2:656265/93‑1 (MQ=255)
                                                     tGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGAt                                   <  1:2166198/93‑1 (MQ=255)
                                                     tGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGAt                                   >  1:3332145/1‑93 (MQ=255)
                                                      gCTGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGAtt                                  <  2:198782/93‑1 (MQ=255)
                                                        tGGAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTg                                >  1:3576234/1‑93 (MQ=255)
                                                          gAACGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTCCCTCTACGCTGCAACAAGCTGCCAGCACCCGGCTTGGATTTGGc                              <  2:72282/93‑1 (MQ=255)
                                                           aaCGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCg                             >  2:790436/1‑93 (MQ=255)
                                                             cGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTg                           <  1:4227325/93‑1 (MQ=255)
                                                             cGTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTg                           <  2:1628056/93‑1 (MQ=255)
                                                              gTGAAGACAAACCGACAACCAGTAACCCCGGCGCTCCTTTTATTGCCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGa                          <  2:2724975/93‑1 (MQ=255)
                                                              gTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGa                          <  1:2280842/93‑1 (MQ=255)
                                                              gTGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGa                          >  1:4128414/1‑93 (MQ=255)
                                                               tGAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGaa                         >  1:1206412/1‑93 (MQ=255)
                                                                gAAGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGaaa                        >  2:1446223/1‑93 (MQ=255)
                                                                  aGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGaaaaa                      <  2:648276/93‑1 (MQ=255)
                                                                  aGACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGaaaaa                      >  2:2635853/1‑93 (MQ=255)
                                                                   gACAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGaaaaaa                     >  2:1889888/1‑93 (MQ=255)
                                                                     cAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAAc                   <  1:1420627/93‑1 (MQ=255)
                                                                     cAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAAc                   <  1:2185282/93‑1 (MQ=255)
                                                                     cAAACCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAAc                   >  1:3203757/1‑93 (MQ=255)
                                                                       aaCCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCa                 <  1:230675/93‑1 (MQ=255)
                                                                        aCCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCat                >  2:3908379/1‑93 (MQ=255)
                                                                        aCCGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCat                <  1:3985047/93‑1 (MQ=255)
                                                                          cGACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCatga              <  2:2682243/93‑1 (MQ=255)
                                                                           gACAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCatgat             <  2:3000657/93‑1 (MQ=255)
                                                                             cAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCatgatga           <  2:1990116/93‑1 (MQ=255)
                                                                             cAACCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCatgatga           <  2:1652454/93‑1 (MQ=255)
                                                                              aaCCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCatgatgat          <  1:2303058/93‑1 (MQ=255)
                                                                              aaCCAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGAGTTGGCGTGAAAAAAACCatgatgat          >  2:506351/1‑93 (MQ=255)
                                                                                ccAGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATgg        <  2:713658/93‑1 (MQ=255)
                                                                                  aGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGgcg      >  1:557578/1‑93 (MQ=255)
                                                                                  aGTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGgcg      >  1:140312/1‑93 (MQ=255)
                                                                                   gTAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGgcgc     <  2:4177434/93‑1 (MQ=255)
                                                                                   ctAAACCCGCCGCTCTTCTTATTACCTCTACGCTGCAACAAGCTGCAAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGgcgc     <  1:2822855/92‑1 (MQ=255)
                                                                                    tAAACCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCa    >  2:1197239/1‑93 (MQ=255)
                                                                                     aaaCCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAg   >  1:19681/1‑93 (MQ=255)
                                                                                      aaCCCGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGc  <  2:3390097/93‑1 (MQ=255)
                                                                                          |                                                                                        
CAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGC  >  NC_000913/1329821‑1329999

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