Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,486 A→T 6.4% K541N (AAA→AAT) ‡ 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,4860AT 6.4% 282.5 / 6.0 92K541N (AAA→AAT) ‡recQATP‑dependent DNA helicase RecQ
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base T (4/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGC  >  NC_000913/4007337‑4007633
                                                                                                                                                     |                                                                                                                                                   
cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCaaa                                                                                                                                                     <  1:1173677/150‑1 (MQ=255)
   cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTg                                                                                                                                                  <  1:209057/150‑1 (MQ=255)
   cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTg                                                                                                                                                  <  1:1004267/150‑1 (MQ=255)
     aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCg                                                                                                                                                <  1:1212913/150‑1 (MQ=255)
     aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCg                                                                                                                                                >  1:1163914/1‑150 (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:99221/150‑1 (MQ=255)
           agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                          <  1:1360726/150‑1 (MQ=255)
           agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGGTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                          >  1:1076371/1‑150 (MQ=255)
              ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttggc                                                                                                                                       <  1:31684/150‑7 (MQ=255)
              ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGAt                                                                                                                                       <  2:590061/150‑1 (MQ=255)
              ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGAt                                                                                                                                       >  2:1105775/1‑150 (MQ=255)
               gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATa                                                                                                                                      <  1:694565/150‑1 (MQ=255)
                ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAg                                                                                                                                     >  2:1512125/1‑150 (MQ=255)
                ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAg                                                                                                                                     <  1:586312/150‑1 (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                                                                                                                      >  1:1148822/1‑150 (MQ=255)
                               tgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATg                                                                                                                      >  2:829518/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                                                                                                         >  1:545840/1‑150 (MQ=255)
                                            ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgt                                                                                                         >  2:557720/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:241622/1‑150 (MQ=255)
                                                                             gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                        <  1:266213/150‑1 (MQ=255)
                                                                             gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                        <  1:829518/150‑1 (MQ=255)
                                                                              cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                        <  1:1043543/149‑1 (MQ=255)
                                                                              cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                        >  2:1043543/1‑149 (MQ=255)
                                                                              cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCt                                                                       >  1:1494240/1‑150 (MQ=255)
                                                                               tCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTg                                                                      >  1:604616/1‑150 (MQ=255)
                                                                                cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                     >  1:873425/1‑150 (MQ=255)
                                                                                cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                     >  1:398774/1‑150 (MQ=255)
                                                                                cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                     <  1:1216802/150‑1 (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                                                  <  1:1157604/150‑1 (MQ=255)
                                                                                                   aaaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGcc                                                  <  2:1084449/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               >  1:28507/1‑150 (MQ=255)
                                                                                                                                      cGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcg               <  2:949071/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
CTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGC  >  NC_000913/4007337‑4007633

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: