Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,487 C→T 6.4% L542L (CTG→TTG)  recQ → ATP‑dependent DNA helicase RecQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,4870CT 6.4% 296.4 / 6.5 93L542L (CTG→TTG) recQATP‑dependent DNA helicase RecQ
Reads supporting (aligned to +/- strand):  ref base C (40/47);  new base T (4/2);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTT  >  NC_000913/4007340‑4007636
                                                                                                                                                   |                                                                                                                                                     
cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTg                                                                                                                                                     <  1:1004267/150‑1 (MQ=255)
cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTg                                                                                                                                                     <  1:209057/150‑1 (MQ=255)
  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCg                                                                                                                                                   >  1:1163914/1‑150 (MQ=255)
  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCg                                                                                                                                                   <  1:1212913/150‑1 (MQ=255)
    tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTa                                                                                                                                                 <  1:442936/150‑1 (MQ=255)
    tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCATTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTa                                                                                                                                                 <  1:1431935/150‑1 (MQ=255)
      aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTaaa                                                                                                                                               >  2:371033/1‑150 (MQ=255)
       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGCAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAAt                                                                                                                                              >  2:1426833/1‑150 (MQ=255)
       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGCTCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAAt                                                                                                                                              >  1:1443912/1‑150 (MQ=255)
        agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                             <  1:1360726/150‑1 (MQ=255)
        agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                             <  1:99221/150‑1 (MQ=255)
        agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGGTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                             >  1:1076371/1‑150 (MQ=255)
           ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttggc                                                                                                                                          <  1:31684/150‑7 (MQ=255)
           ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGAt                                                                                                                                          >  2:1105775/1‑150 (MQ=255)
           ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGAt                                                                                                                                          <  2:590061/150‑1 (MQ=255)
            gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATa                                                                                                                                         <  1:694565/150‑1 (MQ=255)
             ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAg                                                                                                                                        <  1:586312/150‑1 (MQ=255)
             ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAg                                                                                                                                        >  2:1512125/1‑150 (MQ=255)
              cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGc                                                                                                                                       >  1:810321/1‑150 (MQ=255)
                cgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCg                                                                                                                                     >  1:1437605/1‑150 (MQ=255)
                  cgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGAt                                                                                                                                   <  2:1078502/150‑1 (MQ=255)
                    cccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGa                                                                                                                                 <  2:1210536/150‑1 (MQ=255)
                     ccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGaa                                                                                                                                <  2:590664/150‑1 (MQ=255)
                        gTGCGGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGt                                                                                                                             <  2:981422/150‑1 (MQ=255)
                            tgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATg                                                                                                                         >  2:829518/1‑150 (MQ=255)
                            tgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATg                                                                                                                         >  1:1148822/1‑150 (MQ=255)
                             gcgcGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGt                                                                                                                        >  1:809095/1‑150 (MQ=255)
                                cggcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTccc                                                                                                                     >  2:593129/1‑150 (MQ=255)
                                  gcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCccgc                                                                                                                   >  1:600212/1‑150 (MQ=255)
                                   cgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCccgcc                                                                                                                  <  2:810321/150‑1 (MQ=255)
                                      aTCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTa                                                                                                               <  1:371033/150‑1 (MQ=255)
                                         ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgt                                                                                                            >  2:557720/1‑150 (MQ=255)
                                         ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgt                                                                                                            >  1:545840/1‑150 (MQ=255)
                                           ctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtgg                                                                                                          <  2:1483231/150‑1 (MQ=255)
                                            tttGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtggt                                                                                                         >  1:323239/1‑150 (MQ=255)
                                             ttGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtggtg                                                                                                        >  1:358390/1‑150 (MQ=255)
                                               gCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGtt                                                                                                      <  2:358390/150‑1 (MQ=255)
                                                cAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGttt                                                                                                     <  2:378590/150‑1 (MQ=255)
                                                   cTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAac                                                                                                  >  2:1129144/1‑150 (MQ=255)
                                                   cTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAac                                                                                                  <  1:593129/150‑1 (MQ=255)
                                                     tgccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAAATATGATCGCAAACTGTTCGCCAAATTACGGAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacga                                                                                                <  1:806901/150‑1 (MQ=255)
                                                         gtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCa                                                                                            <  1:1143949/150‑1 (MQ=255)
                                                             cgcgTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCt                                                                                        >  1:79490/1‑150 (MQ=255)
                                                             cgcaTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCt                                                                                        <  1:1273765/150‑1 (MQ=255)
                                                                 tATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgat                                                                                    >  2:654628/1‑150 (MQ=255)
                                                                   tCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgattg                                                                                  <  1:1105775/150‑1 (MQ=255)
                                                                       ggCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGAt                                                                              <  2:528445/150‑1 (MQ=255)
                                                                       ggCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGAt                                                                              <  1:90040/150‑1 (MQ=255)
                                                                        gcgcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATg                                                                             <  1:557720/150‑1 (MQ=255)
                                                                          gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                           <  1:829518/150‑1 (MQ=255)
                                                                          gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                           >  1:241622/1‑150 (MQ=255)
                                                                          gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                           <  1:266213/150‑1 (MQ=255)
                                                                           cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                           >  2:1043543/1‑149 (MQ=255)
                                                                           cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                           <  1:1043543/149‑1 (MQ=255)
                                                                           cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCt                                                                          >  1:1494240/1‑150 (MQ=255)
                                                                            tCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTg                                                                         >  1:604616/1‑150 (MQ=255)
                                                                             cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                        <  1:1216802/150‑1 (MQ=255)
                                                                             cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                        >  1:398774/1‑150 (MQ=255)
                                                                             cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                        >  1:873425/1‑150 (MQ=255)
                                                                             cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACACAACCTTGATTGAGATGGCTGa                                                                        <  2:425585/150‑1 (MQ=255)
                                                                              aaaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTgaga                                                                               >  2:1060244/1‑142 (MQ=255)
                                                                              aaaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTgaga                                                                               <  1:1060244/142‑1 (MQ=255)
                                                                              aaaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCAAGATTGAGATGGCTCaa                                                                       >  1:672903/1‑150 (MQ=255)
                                                                                    aaaGCGATGCAGAAATCGTTCGGCGGCAACTATGCTCGCAAACTGTTGGCCAAATTACGCAAACTGCGTAAATCGATACCCGATGTAAGTAATGGCCCGCCGTACGTGGTGTTTAACGCCGGAACCTTGATTGAGATGGCTGAACAGATg                                                                 <  2:857320/150‑1 (MQ=255)
                                                                                     aaGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGc                                                                <  2:1310338/150‑1 (MQ=255)
                                                                                          aTGCAGAAATAGTTCGGCGGCAACTATGATCGCAAAATGTTAGCCAAATTACGCAAACTGAGTAAATAAATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATc                                                           >  2:496441/1‑150 (MQ=255)
                                                                                           tGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCa                                                          >  1:25579/1‑150 (MQ=255)
                                                                                             cAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCAcc                                                        >  1:1273208/1‑150 (MQ=255)
                                                                                                aaaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGcc                                                     <  2:1084449/150‑1 (MQ=255)
                                                                                                aaaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGcc                                                     <  1:1157604/150‑1 (MQ=255)
                                                                                                        cggcggCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAAt                                             >  1:1252877/1‑150 (MQ=255)
                                                                                                                aCTATGATCTCAATCTGTTCGCCAAATTTCGCAAACTGCGTAAATCGATATAGGATTAAAGTAATGTCCCGCCGTAGGTGGTGTTTTACGACGCAAGCTTTATTGAGATGTATTATCAGATTCCGATAACCGCCAGCGAAATGCTCAGCg                                     <  2:1116167/150‑1 (MQ=255)
                                                                                                                 cTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGATAGTAATGTCCCGCCGTACGTGGGGTTTAACGACTCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGt                                    <  2:820845/150‑1 (MQ=255)
                                                                                                                      aTCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACg                               <  2:323239/150‑1 (MQ=255)
                                                                                                                          cAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGt                           >  2:260174/1‑150 (MQ=255)
                                                                                                                          cAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGAAGGCTGAACAGATGCAGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGt                           >  2:1068481/1‑150 (MQ=255)
                                                                                                                           aaaCTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGtt                          <  2:25579/150‑1 (MQ=255)
                                                                                                                              cTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTggg                       <  2:359499/150‑1 (MQ=255)
                                                                                                                                gTTCTCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGAt                     >  2:1253827/1‑150 (MQ=255)
                                                                                                                                   cGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcg                  <  2:949071/150‑1 (MQ=255)
                                                                                                                                   cGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcg                  >  1:28507/1‑150 (MQ=255)
                                                                                                                                    gCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcgc                 >  1:1032777/1‑150 (MQ=255)
                                                                                                                                     tcgatttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCa                >  2:289833/6‑150 (MQ=255)
                                                                                                                                     ccAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATAACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCa                >  2:21703/1‑150 (MQ=255)
                                                                                                                                      cgatttACGCAGTTTGCGTAATTCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCaa               >  1:316681/5‑150 (MQ=255)
                                                                                                                                      cAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCGCAGCGTTAACGGCGTTGGGATGCGCaa               <  2:873425/150‑1 (MQ=255)
                                                                                                                                       gatttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAg              >  1:522401/4‑150 (MQ=255)
                                                                                                                                        atttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGc             >  1:1099023/3‑150 (MQ=255)
                                                                                                                                        aaTTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGc             <  2:79490/150‑1 (MQ=255)
                                                                                                                                         aTTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCt            <  1:571329/150‑1 (MQ=255)
                                                                                                                                           tACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTgg          <  1:1263853/150‑1 (MQ=255)
                                                                                                                                            aCGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGa         >  1:422188/1‑150 (MQ=255)
                                                                                                                                                aaaCTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGc     <  1:1512125/150‑1 (MQ=255)
                                                                                                                                                   cTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCttt  >  2:1322600/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTT  >  NC_000913/4007340‑4007636

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑