Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,309 G→C 6.6% D300H (GAC→CAC)  mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,398,3090GC 6.6% 288.6 / 7.5 85D300H (GAC→CAC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base C (2/4);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCG  >  NC_000913/4398160‑4398457
                                                                                                                                                     |                                                                                                                                                    
tGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTgg                                                                                                                                                      <  2:873575/150‑1 (MQ=255)
tGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTgg                                                                                                                                                      <  1:1307405/150‑1 (MQ=255)
 ggCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGa                                                                                                                                                     >  1:1010760/1‑150 (MQ=255)
 ggCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGa                                                                                                                                                     >  2:1117356/1‑150 (MQ=255)
 ggCAGAAATTCAGTATTGCTACGTGAACCGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGa                                                                                                                                                     >  2:834572/1‑150 (MQ=255)
  gCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGAc                                                                                                                                                    >  1:324117/1‑150 (MQ=255)
  gCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGAc                                                                                                                                                    <  2:1029728/150‑1 (MQ=255)
   cAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACg                                                                                                                                                   <  1:460501/150‑1 (MQ=255)
   cAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACg                                                                                                                                                   >  1:90337/1‑150 (MQ=255)
     gAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTc                                                                                                                                                 >  1:8432/1‑150 (MQ=255)
       aaTTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCaa                                                                                                                                               <  1:484182/150‑1 (MQ=255)
        aTTCCGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAAc                                                                                                                                              >  1:263002/1‑150 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                          <  1:450762/150‑1 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                          <  1:1555731/150‑1 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                          <  1:1290500/150‑1 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                          <  2:1511750/150‑1 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                          <  2:1603791/150‑1 (MQ=255)
            aGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGCACGTTGACGtcc                                                                                                                                          >  1:50098/1‑148 (MQ=255)
             gTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCa                                                                                                                                         <  1:999431/150‑1 (MQ=255)
               aTTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGCACGTTGACGTCCAcc                                                                                                                                       <  2:237002/150‑1 (MQ=255)
                ttGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGCACGTTGACGTCCAcct                                                                                                                                      <  1:313661/150‑2 (MQ=255)
                ttGCTACGTGAAATGTCGCATGATGCGCGATCTCCTTATCAATCACGCGATACGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGCACGTTGACGTCCAcct                                                                                                                                      <  2:50098/150‑2 (MQ=255)
                 tGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCAcccc                                                                                                                                     <  1:514366/150‑1 (MQ=255)
                   cTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGc                                                                                                                                   >  2:996770/1‑150 (MQ=255)
                   cTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGCACGTTGACGTCCAcctgat                                                                                                                                   <  1:24438/150‑5 (MQ=255)
                          aaCGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAacac                                                                                                                            >  2:382834/1‑150 (MQ=255)
                                 gCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGc                                                                                                                     <  1:1117356/150‑1 (MQ=255)
                                      tttcgctATCGCGTGATCATTAACGCGATCCGCCATTCCTTCGAAGAGATACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTc                                                                                                                <  2:254969/147‑1 (MQ=255)
                                      atgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGttta                                                                                                                <  1:834572/150‑2 (MQ=255)
                                      atgCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTc                                                                                                                >  1:78566/1‑150 (MQ=255)
                                         cgcgATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCAt                                                                                                             >  2:907671/1‑150 (MQ=255)
                                           ggATCGCGTGATTGATCAGGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCa                                                                                                           >  2:356721/2‑150 (MQ=255)
                                              tCGCGTGATTGATCAGGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTc                                                                                                        <  1:356721/150‑1 (MQ=255)
                                               cGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTcg                                                                                                       <  1:1106527/150‑1 (MQ=255)
                                                      tcaatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTgg                                                                                                <  2:263002/150‑1 (MQ=255)
                                                       caatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTggg                                                                                               >  2:192932/1‑149 (MQ=255)
                                                        aatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTg                                                                                              <  1:680011/150‑1 (MQ=255)
                                                        aatcaCGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTg                                                                                              <  1:977645/150‑1 (MQ=255)
                                                            aCGCGATCCGCCAGGCCTGCGAAGACAAACTGGCGGCCGATGAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATg                                                                                          <  2:688754/150‑1 (MQ=255)
                                                              gcgATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGAt                                                                                        <  2:1343635/150‑1 (MQ=255)
                                                                gATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGAttt                                                                                      <  2:609961/150‑1 (MQ=255)
                                                                 aTCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCGATCAGTCGCGTCTGGTTCATGAtttt                                                                                     <  2:773283/150‑1 (MQ=255)
                                                                    cGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTtatc                                                                                  <  2:496609/150‑1 (MQ=255)
                                                                          gCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAg                                                                            <  2:1010760/150‑1 (MQ=255)
                                                                                gAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGTAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTg                                                                      <  2:419705/150‑1 (MQ=255)
                                                                                 aaGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGc                                                                     <  1:595441/150‑1 (MQ=255)
                                                                                  aGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCt                                                                    <  1:762663/150‑1 (MQ=255)
                                                                                       aaCTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCg                                                               <  1:382834/150‑1 (MQ=255)
                                                                                          tGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGc                                                            <  2:49379/150‑1 (MQ=255)
                                                                                          tGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGc                                                            <  1:978668/150‑1 (MQ=255)
                                                                                           gggggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCt                                                           <  2:1120890/150‑1 (MQ=255)
                                                                                            ggggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTa                                                          >  2:721170/1‑150 (MQ=255)
                                                                                             gggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTac                                                         <  2:185629/150‑1 (MQ=255)
                                                                                             gggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTac                                                         >  1:882001/1‑150 (MQ=255)
                                                                                                 cGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTacaaca                                                     >  2:736080/1‑150 (MQ=255)
                                                                                                  gATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGCCCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAg                                                    >  1:424892/1‑150 (MQ=255)
                                                                                                  gATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAg                                                    >  2:616834/1‑150 (MQ=255)
                                                                                                  gATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAg                                                    <  2:90337/150‑1 (MQ=255)
                                                                                                  gATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAg                                                    >  1:810288/1‑150 (MQ=255)
                                                                                                       gCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACt                                               >  1:1505780/1‑150 (MQ=255)
                                                                                                        cAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTg                                              >  2:514940/1‑150 (MQ=255)
                                                                                                           ccGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGCCGTCAACGTGCCCCCCGCCAAACACGATGTGCGTTTCCATCAGTCGCGTCTGGTGCATGTTTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGaa                                           >  2:1545670/1‑150 (MQ=255)
                                                                                                                    gtgtTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTa                                  >  1:398091/1‑150 (MQ=255)
                                                                                                                      gttgtATCTGGAGATCGACCCACATCAGGTGGACTTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTAcc                                <  1:1158983/150‑1 (MQ=255)
                                                                                                                      gttgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTAcc                                >  1:711861/1‑150 (MQ=255)
                                                                                                                      gttgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTAcc                                <  1:996770/150‑1 (MQ=255)
                                                                                                                      gttgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTAcc                                >  2:470219/1‑150 (MQ=255)
                                                                                                                        tgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTtatctatc                                                                              >  1:414577/1‑102 (MQ=255)
                                                                                                                        tgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTtatctatc                                                                              <  2:414577/102‑1 (MQ=255)
                                                                                                                        tgtATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGc                              >  1:1234724/1‑150 (MQ=255)
                                                                                                                          tATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTg                            >  2:73599/1‑150 (MQ=255)
                                                                                                                             cTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTacaaca                                                     <  1:103254/122‑1 (MQ=255)
                                                                                                                             cTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTacaaca                                                     >  2:103254/1‑122 (MQ=255)
                                                                                                                             cTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGAc                         <  2:8432/150‑1 (MQ=255)
                                                                                                                              tGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACg                        >  2:978310/1‑150 (MQ=255)
                                                                                                                                 agaTCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATg                     <  2:398091/150‑1 (MQ=255)
                                                                                                                                  gaTCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGa                    >  2:1544272/1‑150 (MQ=255)
                                                                                                                                     cGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAAcc                 >  2:754024/1‑150 (MQ=255)
                                                                                                                                        cccACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCa              >  1:206230/1‑150 (MQ=255)
                                                                                                                                            caTCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGTGCGTGCTGAGCGTGCTACAACAGCAACTCGAAACGACGCTACCGCTGGACGATGAACCCCAAcca          >  2:1378195/1‑149 (MQ=255)
                                                                                                                                            caTCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAAcgccgc                                    <  2:1131328/124‑1 (MQ=255)
                                                                                                                                            caTCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAAcgccgc                                    >  1:1131328/1‑124 (MQ=255)
                                                                                                                                            caTCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTCGAAACGCCGCTACCGCTGGACGATGAACCCCAACCt          >  1:1234977/1‑150 (MQ=255)
                                                                                                                                             aTCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTg         >  2:1359174/1‑150 (MQ=255)
                                                                                                                                              tCAGGTGCACGTCAACGTGCACCCCGCCAAACACGAAGTGAGTTTCCATCAGTAGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCCTGCTACAACATCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCACCCTGc        >  2:767778/1‑150 (MQ=255)
                                                                                                                                                 ggTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCAcc     >  1:18853/1‑150 (MQ=255)
                                                                                                                                                  gTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACcg    <  2:78566/150‑1 (MQ=255)
                                                                                                                                                    ggACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACcgcg  >  1:1187949/1‑150 (MQ=255)
                                                                                                                                                    ggACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACcgcg  >  2:592747/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCG  >  NC_000913/4398160‑4398457

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: