Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,859,142 G→A A684T (GCC→ACC)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,859,1420GA98.8% 246.1 / NA 81A684T (GCC→ACC) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (25/54);  minor base T (0/1);  total (25/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTA  >  NC_000913/2859008‑2859278
                                                                                                                                      |                                                                                                                                        
aTGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACaccacc                                                                                                                                     >  2:912373/1‑140 (MQ=255)
  gTACCGGCACAAACAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGa                                                                                                                                   <  2:118150/140‑1 (MQ=255)
  gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGTCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACTCCACCGa                                                                                                                                   <  2:917399/140‑1 (MQ=255)
  gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGa                                                                                                                                   <  2:917590/140‑1 (MQ=255)
   tACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGaa                                                                                                                                  >  2:478267/1‑140 (MQ=255)
    aCCGGCACAAAAAGTCGAGATTGGACCTATCGCTCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAAt                                                                                                                                 <  1:1337477/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTTCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:163024/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:1209158/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:1209353/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  2:478422/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:669698/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  2:827563/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:700333/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  2:977770/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:702954/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:1659348/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:164411/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:163001/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  1:162494/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                >  2:1283100/1‑140 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTGTACATAACACCACCGAATa                                                                                                                                <  1:669688/140‑1 (MQ=255)
     ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTAAACCTTTATGGTGGAGATGACTGATACAGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  2:827187/140‑1 (MQ=255)
          acaAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATTGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTc                                                                                                                           <  2:489343/140‑1 (MQ=255)
             aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTgg                                                                                                                        >  2:1442994/1‑140 (MQ=255)
             aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTgg                                                                                                                        >  2:1442703/1‑140 (MQ=255)
                 gTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGTGTTGTCTGTTTCCGCCAATATTTTTCTTAACACCACCTATTACATTGTGGTGtt                                                                                                                    <  1:328712/140‑1 (MQ=255)
                  tCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTa                                                                                                                   <  1:740042/140‑1 (MQ=255)
                   cGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTaa                                                                                                                  <  2:183482/140‑1 (MQ=255)
                   cGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTaa                                                                                                                  <  2:293613/140‑1 (MQ=255)
                            aCCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                <  2:1575914/117‑1 (MQ=255)
                            aCCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATa                                                                                                                                >  1:1575914/1‑117 (MQ=255)
                             ccTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGAt                                                                                                        <  1:681813/140‑1 (MQ=255)
                             ccTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGAt                                                                                                        <  1:1578689/140‑1 (MQ=255)
                                atcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATAGTTTACAAAACACCACCGAATACAGTCTGGTGTAAATGGAGGAGATCgg                                                                                                     >  2:887622/1‑140 (MQ=255)
                                      cgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTgg                                                                                               >  2:1181503/1‑140 (MQ=255)
                                          tCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACg                                                                                           <  1:574108/140‑1 (MQ=255)
                                          tCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACg                                                                                           <  1:574094/140‑1 (MQ=255)
                                           cTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt                                                                                          >  2:1041503/1‑140 (MQ=255)
                                            tttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATTACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc                                                                                         <  1:176108/140‑1 (MQ=255)
                                            tttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc                                                                                         <  2:1668728/140‑1 (MQ=255)
                                            tttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc                                                                                         <  2:1649428/140‑1 (MQ=255)
                                             ttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                        <  1:1117441/140‑1 (MQ=255)
                                               aCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                      >  2:934713/1‑140 (MQ=255)
                                               aCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                      >  2:937027/1‑140 (MQ=255)
                                                   gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCTCCTAc                                                                                  >  2:316672/1‑140 (MQ=255)
                                                    cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCGGGTGATAATGGATGAGATCGGGCGAGGAACGTCCACCTACg                                                                                 >  2:93036/1‑140 (MQ=255)
                                                          cgcgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTc                                                                           <  2:657257/140‑1 (MQ=255)
                                                              gCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTc                                                                           >  2:984309/1‑136 (MQ=255)
                                                              gCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTc                                                                           <  1:984309/136‑1 (MQ=255)
                                                              gCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                       >  2:71320/1‑140 (MQ=255)
                                                                aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGc                                                                     <  2:1255222/140‑1 (MQ=255)
                                                                aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGc                                                                     <  2:510871/140‑1 (MQ=255)
                                                                       cTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACCGTCTGGTGTTAATGGAGGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtg                                                              >  1:134963/1‑140 (MQ=255)
                                                                       cTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtg                                                              >  1:730416/1‑140 (MQ=255)
                                                                                gggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGTAACTTCAACCTACGATGGTCTGTCGATgtcgttggggttgta                                                     <  1:734142/140‑5 (MQ=255)
                                                                                               aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                      <  1:337165/140‑1 (MQ=255)
                                                                                               aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                      <  1:337237/140‑1 (MQ=255)
                                                                                               aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                      <  2:692456/140‑1 (MQ=255)
                                                                                               aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                      <  2:552018/140‑1 (MQ=255)
                                                                                               aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                      <  1:64582/140‑1 (MQ=255)
                                                                                                    ggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGa                                 >  1:443780/1‑140 (MQ=255)
                                                                                                     gagaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt                                >  1:1646553/1‑140 (MQ=255)
                                                                                                     gagaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt                                >  1:294343/1‑140 (MQ=255)
                                                                                                     gagaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt                                >  1:1461226/1‑140 (MQ=255)
                                                                                                     gagaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt                                >  1:747401/1‑140 (MQ=255)
                                                                                                             tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAt                        >  1:1403598/1‑140 (MQ=255)
                                                                                                              gAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAtt                       <  1:255290/140‑1 (MQ=255)
                                                                                                               aaaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTg                      <  1:618306/140‑1 (MQ=255)
                                                                                                                  ccGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACg                   <  1:542142/140‑1 (MQ=255)
                                                                                                                  ccGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACg                   <  1:934713/140‑1 (MQ=255)
                                                                                                                  ccGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACg                   <  1:937027/140‑1 (MQ=255)
                                                                                                                    gCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGtt                 <  1:978793/140‑1 (MQ=255)
                                                                                                                      cAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGATCGTCCACCTACGTTGGTCTGTCGCTGGCGTGGGAGTGATATTAAAATCTGGCGAATAAGATTATTGCATTgtggtttt               <  1:1167892/140‑8 (MQ=255)
                                                                                                                           ttttACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCt          <  2:1486261/140‑1 (MQ=255)
                                                                                                                            tttACATAACACCACCGAATACAGTCTGGTGTTACTGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTa         <  2:557843/140‑1 (MQ=255)
                                                                                                                                cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAATATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa     <  2:177489/140‑1 (MQ=255)
                                                                                                                                cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa     <  2:1164257/140‑1 (MQ=255)
                                                                                                                                cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa     <  1:535805/140‑1 (MQ=255)
                                                                                                                                cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCa     <  2:1164477/140‑1 (MQ=255)
                                                                                                                                   aaCACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTa  >  1:788602/1‑140 (MQ=255)
                                                                                                                                   aaCACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTa  >  1:788855/1‑140 (MQ=255)
                                                                                                                                   aaCACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTa  >  1:100303/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
ATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTA  >  NC_000913/2859008‑2859278

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