Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,957,139 A→G 5.1% S332P (TCC→CCC)  torY ← TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,957,1390AG 5.1% 246.1 / 2.3 78S332P (TCC→CCC) torYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (2/2);  total (37/41)
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

GTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTGT  >  NC_000913/1956992‑1957275
                                                                                                                                                   |                                                                                                                                        
gTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCGATGGAc                                                                                                                                         >  5:371550/1‑149 (MQ=255)
  cTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACgg                                                                                                                                       >  7:215615/1‑149 (MQ=255)
     gTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGGAATGGACGGCCa                                                                                                                                    >  3:214894/1‑149 (MQ=255)
        aaTGTCATGACTTCCTCACTGATCAAAGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGATCATGCTGGTTCGTGCCCCCATTCCTATGGCAATGGACGGCCACg                                                                                                                                  >  8:68565/1‑148 (MQ=255)
         aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTCGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCa                                                                                                                                >  4:111200/1‑149 (MQ=255)
              aTAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa                                                                                                                            <  5:221012/148‑1 (MQ=255)
              aTAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGAc                                                                                                                           <  5:216478/149‑1 (MQ=255)
              aTAACGTCCTCCCTGATCAACGCGGCTCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa                                                                                                                            <  8:137223/148‑1 (MQ=255)
               tAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                                                                          >  3:164435/1‑149 (MQ=255)
               tAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                                                                          <  4:124549/149‑1 (MQ=255)
                aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg                                                                                                                         >  8:99174/1‑149 (MQ=255)
                       tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATg                                                                                                                  >  7:322422/1‑149 (MQ=255)
                         ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgtg                                                                                                                       <  2:331343/142‑1 (MQ=255)
                         ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgtg                                                                                                                       >  1:331343/1‑142 (MQ=255)
                         ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                                                                                >  8:312183/1‑149 (MQ=255)
                         ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGGGTAATGGt                                                                                                                >  3:125138/1‑149 (MQ=255)
                           tGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATg                                                                                                                  >  1:102411/1‑145 (MQ=255)
                           tGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATg                                                                                                                  <  2:102411/145‑1 (MQ=255)
                            gATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                                                                             >  1:66932/1‑149 (MQ=255)
                              tCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcg                                                                                                           <  8:302705/149‑1 (MQ=255)
                                aaCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcgg                                                                                                          <  7:195440/148‑1 (MQ=255)
                                 aCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCggcggc                                                                                                         >  3:354533/1‑148 (MQ=255)
                                    aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgt                                                                                                                        >  2:274519/1‑130 (MQ=255)
                                    aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgt                                                                                                                        <  1:274519/130‑1 (MQ=255)
                                      gATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAgg                                                                                                   >  2:23487/1‑149 (MQ=255)
                                       aTCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAggg                                                                                                  <  4:354533/149‑1 (MQ=255)
                                       aTCACTGTTTATCGTTAATATGTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAggg                                                                                                  <  2:384129/149‑1 (MQ=255)
                                          aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGGCGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                                                                               >  1:97541/1‑149 (MQ=255)
                                          aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGGCGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                                                                               >  1:92453/1‑149 (MQ=255)
                                           cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                                                                              <  2:143395/149‑1 (MQ=255)
                                           cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                                                                              <  6:202753/149‑1 (MQ=255)
                                           cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                                                                              <  1:255941/149‑1 (MQ=255)
                                              tttCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGAc                                                                                           <  8:215615/149‑1 (MQ=255)
                                                 ctcGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAAc                                                                                        >  6:18191/1‑149 (MQ=255)
                                                    ggTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAg                                                                                     >  7:101344/1‑149 (MQ=255)
                                                     gTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                                                                                    <  7:312183/149‑1 (MQ=255)
                                                     gTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                                                                                    <  5:18191/149‑1 (MQ=255)
                                                           tCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCGATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAAt                                                                               <  6:371550/148‑1 (MQ=255)
                                                           tCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCaa                                                                                >  7:265040/1‑147 (MQ=255)
                                                             tttGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCTGCGGCAATAGGGGCAAGCCACCCAGCGCAAGAgg                                                                            >  6:153755/1‑149 (MQ=255)
                                                               tGGCGTGGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAg                                                                             <  5:232317/146‑1 (MQ=255)
                                                                 gCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCGACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATc                                                                        <  7:68565/149‑1 (MQ=255)
                                                                       tACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCaa                                                                   <  3:111200/148‑1 (MQ=255)
                                                                       tACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCaa                                                                   <  8:322422/148‑1 (MQ=255)
                                                                         cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                                                                             <  4:206030/104‑1 (MQ=255)
                                                                         cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGt                                                                 <  1:384071/148‑1 (MQ=255)
                                                                         cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGAGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                                                                             <  4:201101/104‑1 (MQ=255)
                                                                         cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGAGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                                                                             >  3:206030/1‑104 (MQ=255)
                                                                         cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGAGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                                                                             >  3:201101/1‑104 (MQ=255)
                                                                            gAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt                                                             >  1:224944/1‑149 (MQ=255)
                                                                              aaTACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa                                                           <  8:64969/149‑1 (MQ=255)
                                                                                   cgcgTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                                                                              <  8:138789/109‑1 (MQ=255)
                                                                                   cgcgTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                                                                              >  7:138789/1‑109 (MQ=255)
                                                                                   cgcgTTAAAATGTCAAGTTCGTTTTCGCTCATGCTGGTTCGTGACCCCATTCCTTAGGCAATGGACGGCAACGCATTGACGGTGTAATGCTCGGCGCCAATATGGTAATGACAACCAGCGCAAAAGGTATCGTCAAGTTTTTCAGCGTa                                                      >  8:14208/1‑149 (MQ=255)
                                                                                         aaaaTGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGGCGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGcc                                                 <  2:97541/148‑1 (MQ=255)
                                                                                         aaaaTGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGGCGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGcc                                                 <  2:92453/148‑1 (MQ=255)
                                                                                          caaTGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCAt                                               <  4:214894/148‑1 (MQ=255)
                                                                                               tCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa                                                           >  1:199621/1‑132 (MQ=255)
                                                                                               tCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa                                                           <  2:199621/132‑1 (MQ=255)
                                                                                                 cAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGtt                                        >  7:135232/1‑149 (MQ=255)
                                                                                                 cAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGtt                                        <  4:164435/149‑1 (MQ=255)
                                                                                                    ttCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTAc                                     <  4:71461/149‑1 (MQ=255)
                                                                                                        ttttCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTc                                 <  8:265040/149‑1 (MQ=255)
                                                                                                            cGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGc                              <  7:99174/148‑1 (MQ=255)
                                                                                                             gCTCATGCTGGTTCGTGCCCCCATTCCTTTGGTAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGGTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCt                             >  8:102310/1‑148 (MQ=255)
                                                                                                               tCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAAt                          >  7:196826/1‑149 (MQ=255)
                                                                                                                  tGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCg                       <  5:152931/149‑1 (MQ=255)
                                                                                                                      ggTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAGGCATAGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCACCGTATTGCCATAGCGGTTTACAGTCCGGAAATGCGGGAt                   >  5:123602/1‑149 (MQ=255)
                                                                                                                       gTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGc                              <  7:390508/137‑1 (MQ=255)
                                                                                                                       gTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGc                              >  8:390508/1‑137 (MQ=255)
                                                                                                                       gTTCGGGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATc                  <  3:280359/149‑1 (MQ=255)
                                                                                                                          cGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCgg                <  2:224944/148‑1 (MQ=255)
                                                                                                                            tGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTaa             <  7:14208/149‑1 (MQ=255)
                                                                                                                               ccccATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCgg                      <  8:150708/137‑1 (MQ=255)
                                                                                                                               ccccATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCgg                      >  7:150708/1‑137 (MQ=255)
                                                                                                                                    ttCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCAccc     >  5:141431/1‑149 (MQ=255)
                                                                                                                                      ccTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTg   >  4:258612/1‑149 (MQ=255)
                                                                                                                                        tttGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTGt  >  7:93470/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
GTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTGT  >  NC_000913/1956992‑1957275

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