Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,964,373 C→T 6.0% P261P (CCG→CCA thyA ← thymidylate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,964,3730CT 6.0% 222.2 / 3.7 67P261P (CCG→CCAthyAthymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base T (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

TTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATAT  >  NC_000913/2964225‑2964521
                                                                                                                                                    |                                                                                                                                                    
tttCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCAcc                                                                                                                                                      >  2:1102581/1‑149 (MQ=255)
     tcGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCggtgc                                                                                                                                                 >  3:62477/1‑146 (MQ=255)
      cGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCt                                                                                                                                                <  3:200388/149‑1 (MQ=255)
             tCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGc                                                                                                                                         <  4:806501/149‑1 (MQ=255)
               ttCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGccc                                                                                                                                       >  2:884274/1‑149 (MQ=255)
                  cGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGa                                                                                                                                    >  3:281835/1‑149 (MQ=255)
                   ggACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGAt                                                                                                                                   <  4:1119934/149‑1 (MQ=255)
                    gACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGTGCTTTAATGCCCGGATg                                                                                                                                  <  4:62477/149‑1 (MQ=255)
                               ccGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGc                                                                                                                       <  3:950435/149‑1 (MQ=255)
                                 gTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCCCCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCt                                                                                                                     <  1:233362/149‑1 (MQ=255)
                                    gCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCa                                                                                                                  >  1:787503/1‑149 (MQ=255)
                                      aGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAAt                                                                                                                >  1:1036127/1‑149 (MQ=255)
                                              gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCgg                                                                                                                               <  2:814329/126‑1 (MQ=255)
                                              gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCgg                                                                                                                               >  1:814329/1‑126 (MQ=255)
                                              gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGt                                                                                                        >  2:904453/1‑149 (MQ=255)
                                              gCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGt                                                                                                        >  1:392347/1‑149 (MQ=255)
                                                aaaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCt                                                                                                      >  1:705751/1‑149 (MQ=255)
                                                 aaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCtt                                                                                                     <  4:277405/149‑1 (MQ=255)
                                                 aaaTTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCtt                                                                                                     <  4:487513/149‑1 (MQ=255)
                                                    tttCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGa                                                                                                  >  2:861513/1‑149 (MQ=255)
                                                      tCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGaaa                                                                                                <  2:787503/149‑1 (MQ=255)
                                                        gggAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACg                                                                                              >  1:1055028/1‑149 (MQ=255)
                                                          gAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGt                                                                                            >  1:364325/1‑149 (MQ=255)
                                                           aaGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTa                                                                                           <  1:1102581/149‑1 (MQ=255)
                                                                 tctcGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACTGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGaa                                                                                     >  4:322157/1‑149 (MQ=255)
                                                                 tctcGAAGAATTGAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGaa                                                                                     >  3:84465/1‑149 (MQ=255)
                                                                   tcGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGa                                                                                   <  2:1055028/149‑1 (MQ=255)
                                                                    cGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGAt                                                                                  <  2:364325/149‑1 (MQ=255)
                                                                     gaagaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATg                                                                                 >  2:919620/1‑149 (MQ=255)
                                                                       agaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGa                                                                               <  4:84465/149‑1 (MQ=255)
                                                                        gaaTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGAt                                                                              <  4:281835/149‑1 (MQ=255)
                                                                          aTTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTc                                                                            <  2:1036127/149‑1 (MQ=255)
                                                                           tttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCg                                                                           >  2:486677/1‑149 (MQ=255)
                                                                           tttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCg                                                                           >  3:230620/1‑149 (MQ=255)
                                                                            ttAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCgg                                                                          <  1:904453/149‑1 (MQ=255)
                                                                                     ggTAAAAAAACCGACGTACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGtt                                                                 >  3:523913/1‑149 (MQ=255)
                                                                                        aaaaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGa                                                              <  1:884274/149‑1 (MQ=255)
                                                                                         aaaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGAt                                                             >  3:136740/1‑149 (MQ=255)
                                                                                          aaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGc                                                                                                                                 >  3:520194/1‑80 (MQ=255)
                                                                                          aaaaaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGc                                                                                                                                 <  4:520194/80‑1 (MQ=255)
                                                                                             aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATc                                                                                                                            >  2:926528/1‑82 (MQ=255)
                                                                                             aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATc                                                                                                                            <  1:926528/82‑1 (MQ=255)
                                                                                             aaCCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATc                                                         >  4:810184/1‑149 (MQ=255)
                                                                                                cGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAAc                                                      >  1:820142/1‑149 (MQ=255)
                                                                                                  aCGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACtt                                                    >  1:963169/1‑149 (MQ=255)
                                                                                                    gCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCg                                                  <  1:486677/149‑1 (MQ=255)
                                                                                                      acacTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGc                                                >  3:1146157/1‑149 (MQ=255)
                                                                                                      acacTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGc                                                >  1:880353/1‑149 (MQ=255)
                                                                                                         cTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTTAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGc                                             <  2:820142/149‑1 (MQ=255)
                                                                                                                 ggCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCgg                                     <  4:1137216/149‑1 (MQ=255)
                                                                                                                 ggCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCgg                                     >  1:242610/1‑149 (MQ=255)
                                                                                                                  gCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGt                                    >  2:36131/1‑149 (MQ=255)
                                                                                                                   ctctGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGtt                                   >  1:798459/1‑149 (MQ=255)
                                                                                                                      tGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTcgc                                >  2:716005/1‑149 (MQ=255)
                                                                                                                        gCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCgg                              >  1:253608/1‑149 (MQ=255)
                                                                                                                            gATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTa                          >  3:887305/1‑149 (MQ=255)
                                                                                                                               gTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAAtt                       <  4:1146157/149‑1 (MQ=255)
                                                                                                                                  tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCTGACGCGGTTCGCGGCTTAATTGCa                    >  2:750432/1‑149 (MQ=255)
                                                                                                                                  tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCa                    <  4:136740/149‑1 (MQ=255)
                                                                                                                                     tAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTACAGAt                 <  2:392347/149‑1 (MQ=255)
                                                                                                                                       aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGa               <  1:919620/149‑1 (MQ=255)
                                                                                                                                       aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGCCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGa               <  2:290313/149‑1 (MQ=255)
                                                                                                                                           gATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTg          >  3:1066415/1‑149 (MQ=255)
                                                                                                                                           gATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGCCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGGTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGttt           >  1:655549/1‑149 (MQ=255)
                                                                                                                                              aGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGAt        <  2:880353/149‑1 (MQ=255)
                                                                                                                                               gCCACTGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATc       <  3:322157/149‑1 (MQ=255)
                                                                                                                                               gCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATc       <  1:861513/149‑1 (MQ=255)
                                                                                                                                                  aCCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGTGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCat    >  3:704831/1‑149 (MQ=255)
                                                                                                                                                    cGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGTGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCatat  <  4:704831/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAGGACGTTGCAAAATTTCGGGAAGGCGTCTCGAAGAATTTAACGGAGGGTAAAAAAACCGACGCACACTGGCGTCGGCTCTGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATAT  >  NC_000913/2964225‑2964521

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