Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,587,133 C→T A62V (GCG→GTG)  mrr → methylated adenine and cytosine restriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,1330CT98.3% 204.4 / NA 60A62V (GCG→GTG) mrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/36);  minor base A (1/0);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCT  >  NC_000913/4586996‑4587267
                                                                                                                                         |                                                                                                                                      
tCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGgtgc                                                                                                                                      >  1:1078494/1‑140 (MQ=255)
tCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGgtgc                                                                                                                                      >  1:584238/1‑140 (MQ=255)
tCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGgtgc                                                                                                                                      >  2:903988/1‑140 (MQ=255)
tCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGgtgc                                                                                                                                      >  2:1363534/1‑140 (MQ=255)
tCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGgtgc                                                                                                                                      >  2:125639/1‑140 (MQ=255)
  tGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCAt                                                                                                                                    <  2:584238/140‑1 (MQ=255)
     caacaaAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGAc                                                                                                                                 <  1:463917/140‑1 (MQ=255)
          aaaCCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGttt                                                                                                                            >  2:175220/1‑140 (MQ=255)
              cGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTaaaa                                                                                                                        >  2:1042280/1‑140 (MQ=255)
                       cAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCggg                                                                                                               >  2:168696/1‑140 (MQ=255)
                       cAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCggg                                                                                                               <  2:749252/140‑1 (MQ=255)
                              gcgTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGc                                                                                                        <  2:882563/140‑1 (MQ=255)
                               cgTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCa                                                                                                       >  2:797744/1‑140 (MQ=255)
                               cgTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCa                                                                                                       <  2:777184/140‑1 (MQ=255)
                               cgTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCa                                                                                                       <  1:1238727/140‑1 (MQ=255)
                               cgTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCa                                                                                                       <  2:256162/140‑1 (MQ=255)
                                        cATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTc                                                                                              >  1:362552/1‑140 (MQ=255)
                                          tGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTcgc                                                                                            <  1:168696/140‑1 (MQ=255)
                                          tGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTcgc                                                                                            <  1:841745/140‑1 (MQ=255)
                                             ggCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTg                                                                                         <  2:362552/140‑1 (MQ=255)
                                               cTGCCGCGGATGCATTAGGACTTGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGc                                                                                       >  1:401096/1‑140 (MQ=255)
                                                tGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCa                                                                                      >  2:546898/1‑140 (MQ=255)
                                                  ccGCGGATGCATTAGGACTGGATGACAGCCAGCTAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCaaa                                                                                    >  1:509378/1‑140 (MQ=255)
                                                  ccGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCaaa                                                                                    >  1:320816/1‑140 (MQ=255)
                                                  ccGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCaaa                                                                                    >  2:511043/1‑140 (MQ=255)
                                                  ccGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCaaa                                                                                    >  2:859622/1‑140 (MQ=255)
                                                   cgcgGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAAt                                                                                   <  1:59139/140‑1 (MQ=255)
                                                       gATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTg                                                                               >  1:1157647/1‑140 (MQ=255)
                                                          gCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCt                                                                            <  1:511043/140‑1 (MQ=255)
                                                          gCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCt                                                                            <  1:53174/140‑1 (MQ=255)
                                                             ttAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGAc                                                                         <  1:175220/140‑1 (MQ=255)
                                                                  aCTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTg                                                                    >  2:2121/1‑140 (MQ=255)
                                                                  aCTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTg                                                                    <  1:2121/140‑1 (MQ=255)
                                                                  aCTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTg                                                                    <  2:1157647/140‑1 (MQ=255)
                                                                     ggATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGAGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCgg                                                                 >  2:1230128/1‑140 (MQ=255)
                                                                      gATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCggg                                                                >  2:139591/1‑140 (MQ=255)
                                                                      gATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCggg                                                                <  1:903988/140‑1 (MQ=255)
                                                                      gATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCggg                                                                >  2:1164681/1‑140 (MQ=255)
                                                                       aTGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGt                                                               <  1:546898/140‑1 (MQ=255)
                                                                       aTGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGt                                                               <  1:631441/140‑1 (MQ=255)
                                                                           cagccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTg                                                           >  1:146951/1‑140 (MQ=255)
                                                                           cagccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTg                                                           >  1:307183/1‑140 (MQ=255)
                                                                             gccagcGAGCGAAAGTCATTACCAGCTGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGAc                                                         <  2:90492/140‑1 (MQ=255)
                                                                             gccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGAc                                                         <  1:859622/140‑1 (MQ=255)
                                                                             gccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGAc                                                         <  1:759389/140‑1 (MQ=255)
                                                                             gccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGAc                                                         <  1:734201/140‑1 (MQ=255)
                                                                             gccagcGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGAc                                                         <  1:1363534/140‑1 (MQ=255)
                                                                                agcgagcgAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTgg                                                      >  1:58235/1‑140 (MQ=255)
                                                                                 gcgagcgaAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTggg                                                     <  2:401096/140‑1 (MQ=255)
                                                                                  cgagcgaAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGt                                                    <  2:307183/140‑1 (MQ=255)
                                                                                     gcgaAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGc                                                 <  1:909560/140‑1 (MQ=255)
                                                                                          aGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGtctc                                            <  1:139591/140‑1 (MQ=255)
                                                                                          aGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGtctc                                            <  2:550889/140‑1 (MQ=255)
                                                                                               ttACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATccc                                       <  1:1164681/140‑1 (MQ=255)
                                                                                               ttACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATccc                                       >  2:816761/1‑140 (MQ=255)
                                                                                               ttACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATccc                                       <  1:1042280/140‑1 (MQ=255)
                                                                                               ttACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATccc                                       <  1:797744/140‑1 (MQ=255)
                                                                                                      cGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCaa                                >  2:654712/1‑140 (MQ=255)
                                                                                                                gTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCa                      <  2:1078494/140‑1 (MQ=255)
                                                                                                                gTTTATAAAAATCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCa                      >  1:706410/1‑140 (MQ=255)
                                                                                                                        aaaTCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGa              <  1:816761/140‑1 (MQ=255)
                                                                                                                          aTCGTGCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAAc            >  2:1362746/1‑140 (MQ=255)
                                                                                                                               gCAGGCTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCt       >  1:200151/1‑140 (MQ=255)
                                                                                                                                    cTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCt  <  2:380700/140‑1 (MQ=255)
                                                                                                                                    cTGGGTGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCt  <  2:965669/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TCTGGCAACAAAACCGGAAGGTGCAGCCGCGCGTGATGTTCATGAGGCTGCCGCGGATGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCT  >  NC_000913/4586996‑4587267

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