Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,331,930 A→G T295A (ACC→GCC)  topA → DNA topoisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,9300AG98.6% 230.0 / ‑4.5 69T295A (ACC→GCC) topADNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/32);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATC  >  NC_000913/1331798‑1332064
                                                                                                                                    |                                                                                                                                      
aaaGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCt                                                                                                                                 <  1:484070/140‑1 (MQ=255)
    aaCAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGa                                                                                                                             >  2:1052884/1‑140 (MQ=255)
     aCAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGAt                                                                                                                            >  2:875851/1‑140 (MQ=255)
              ggCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGa                                                                                                                   >  1:1392821/1‑140 (MQ=255)
                  gCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGaaaaa                                                                                                               <  1:1044958/140‑1 (MQ=255)
                   cGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGaaaaaa                                                                                                              >  1:90397/1‑140 (MQ=255)
                          tCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGATAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCatga                                                                                                       <  2:851218/140‑1 (MQ=255)
                             gctgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCatgatga                                                                                                    >  1:179591/1‑140 (MQ=255)
                               tgGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGCCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCatgatgatg                                                                                                  >  1:858886/1‑140 (MQ=255)
                                        cgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGt                                                                                         <  2:510999/140‑1 (MQ=255)
                                        cgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGt                                                                                         <  1:875851/140‑1 (MQ=255)
                                        cgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGt                                                                                         <  1:818359/140‑1 (MQ=255)
                                        cgcgcTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGt                                                                                         >  1:695842/1‑140 (MQ=255)
                                                  cGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGTCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAg                                                                               <  2:162043/140‑1 (MQ=255)
                                                  cGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAg                                                                               <  2:1392821/140‑1 (MQ=255)
                                                  cGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAg                                                                               <  2:858886/140‑1 (MQ=255)
                                                       tGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGc                                                                          >  2:1097880/1‑140 (MQ=255)
                                                              tGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTa                                                                                    <  2:88542/123‑1 (MQ=255)
                                                              tGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTa                                                                                    >  1:88542/1‑123 (MQ=255)
                                                              tGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCa                                                                   >  2:339701/1‑140 (MQ=255)
                                                                   aCAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTAc                                                              >  1:1319271/1‑140 (MQ=255)
                                                                         cGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGt                                                        <  2:1375654/140‑1 (MQ=255)
                                                                          gACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTa                                                       <  2:647589/140‑1 (MQ=255)
                                                                              aCCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGa                                                   <  1:222326/140‑1 (MQ=255)
                                                                              aCCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGa                                                   <  2:1319271/140‑1 (MQ=255)
                                                                              aCCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGAAGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGa                                                   <  2:951666/140‑1 (MQ=255)
                                                                               ccAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGAc                                                  >  2:1252565/1‑140 (MQ=255)
                                                                                  gTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTcc                                               >  1:954419/1‑140 (MQ=255)
                                                                                      aCCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCAc                                             >  2:512039/1‑138 (MQ=255)
                                                                                      aCCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCAc                                             <  1:512039/138‑1 (MQ=255)
                                                                                          ggCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCt                                       >  1:857230/1‑140 (MQ=255)
                                                                                          ggCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCt                                       >  1:495858/1‑140 (MQ=255)
                                                                                          ggCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCt                                       >  2:60658/1‑140 (MQ=255)
                                                                                           gcgcTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTg                                      >  2:237237/1‑140 (MQ=255)
                                                                                            cgcTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGa                                     >  1:276371/1‑140 (MQ=255)
                                                                                             gcTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAg                                    >  2:44993/1‑140 (MQ=255)
                                                                                                 cTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAg                                >  2:71528/1‑140 (MQ=255)
                                                                                                   tttATTACCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              <  1:1252565/140‑1 (MQ=255)
                                                                                                    ttattaCCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGAc                             <  1:71528/140‑1 (MQ=255)
                                                                                                      attaCCTCTACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGAcgc                           <  2:663919/140‑1 (MQ=255)
                                                                                                           ctctACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGGCTCCACTAACCTGAGTCAGGAcgc                           <  2:1244699/135‑1 (MQ=255)
                                                                                                           ctctACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              >  2:699157/1‑132 (MQ=255)
                                                                                                           ctctACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              <  1:699157/132‑1 (MQ=255)
                                                                                                           ctctACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGAcgc                           >  1:1244699/1‑135 (MQ=255)
                                                                                                              tACGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAGAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAata                   >  1:774746/1‑140 (MQ=255)
                                                                                                               aCGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAatat                  >  1:441094/1‑140 (MQ=255)
                                                                                                                cGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              >  1:756524/1‑127 (MQ=255)
                                                                                                                cGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              >  1:235098/1‑127 (MQ=255)
                                                                                                                cGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              <  2:756524/127‑1 (MQ=255)
                                                                                                                cGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGa                              <  2:235098/127‑1 (MQ=255)
                                                                                                                cGCTGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATg                 >  2:1105656/1‑140 (MQ=255)
                                                                                                                   tGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              <  2:954419/140‑1 (MQ=255)
                                                                                                                   tGCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              <  1:1052884/140‑1 (MQ=255)
                                                                                                                    gCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAGAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc             <  2:774746/140‑1 (MQ=255)
                                                                                                                    gCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc             <  1:60658/140‑1 (MQ=255)
                                                                                                                    gCAACAAGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc             <  1:1105656/140‑1 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTTGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTcg            >  2:1396976/1‑140 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              <  2:667811/138‑1 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              <  2:132242/138‑1 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              >  1:667811/1‑138 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGtt              >  1:132242/1‑138 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTcg            >  2:155789/1‑140 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTcg            >  2:93570/1‑140 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc             <  1:860448/139‑1 (MQ=255)
                                                                                                                     caacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc             >  2:860448/1‑139 (MQ=255)
                                                                                                                      aacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAata                   <  1:246731/132‑1 (MQ=255)
                                                                                                                      aacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAata                   >  2:246731/1‑132 (MQ=255)
                                                                                                                      aacaaGCTGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTcgc           <  2:276371/140‑1 (MQ=255)
                                                                                                                             tGCCAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTtata    >  1:693337/1‑140 (MQ=255)
                                                                                                                               ccAGCGCCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATc  >  2:715169/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                      
AAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATC  >  NC_000913/1331798‑1332064

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