Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,051,723 C→G 6.4% D199H (GAT→CAT)  ubiI ← 2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,051,7230CG 6.4% 500.4 / 20.5 141D199H (GAT→CAT) ubiI2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (66/66);  new base G (7/2);  total (73/68)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 147 total aligned reads displayed.

GCGCGATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACATGCTGCCATCTTTC  >  NC_000913/3051575‑3051871
                                                                                                                                                    |                                                                                                                                                    
gcgcgATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGcatc                                                                                                                                                      <  3:138567/149‑1 (MQ=255)
gaacgATTAAATTCGTCTGCGCTTCACTGCTGCAGTAGCTGCGCTTCCTTTGGCGACAGTGACCTGACAATCGAGCAAAGCTGCGTATCGCTAAGCGGTAAAAAGGCCATAATTCCTTCGCCATGGAAAACCTGCCGCGCGACCGcatc                                                                                                                                                      >  2:800962/4‑149 (MQ=255)
 cgcgATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGcatca                                                                                                                                                     <  4:95351/149‑1 (MQ=255)
  gcgATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGcatcat                                                                                                                                                    <  3:590280/149‑1 (MQ=255)
    gATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGc                                                                                                                                                  <  3:469538/149‑1 (MQ=255)
     aTTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCg                                                                                                                                                 >  1:1071387/1‑149 (MQ=255)
      ttAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCgg                                                                                                                                                <  4:680966/149‑1 (MQ=255)
      ttAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCg                                                                                                                                                 >  1:1196362/1‑148 (MQ=255)
         aaTTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGttc                                                                                                                                             <  2:298190/149‑1 (MQ=255)
          aTTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGttct                                                                                                                                            >  2:540860/1‑149 (MQ=255)
           ttCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGttctt                                                                                                                                           >  2:47964/1‑149 (MQ=255)
           ttCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGttctt                                                                                                                                           >  4:420080/1‑149 (MQ=255)
                    aCTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGaa                                                                                                                                  <  2:299934/149‑1 (MQ=255)
                    aCTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGaa                                                                                                                                  >  1:524459/1‑149 (MQ=255)
                     cTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAAt                                                                                                                                 >  4:603437/1‑149 (MQ=255)
                         ccTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGGGAATGGTc                                                                                                                             >  2:405465/1‑149 (MQ=255)
                                cATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCa                                                                                                                      >  2:406295/1‑149 (MQ=255)
                                    cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAgcgc                                                                                                                  <  1:207634/149‑1 (MQ=255)
                                    cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAgcgc                                                                                                                  >  4:161983/1‑149 (MQ=255)
                                     gCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGcgcg                                                                                                                 >  2:1123219/1‑149 (MQ=255)
                                     gCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGcgcg                                                                                                                 >  1:1047962/1‑149 (MQ=255)
                                     gCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGcgcg                                                                                                                 >  1:366617/1‑149 (MQ=255)
                                       tGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTg                                                                                                               >  1:1273499/1‑149 (MQ=255)
                                       tGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTg                                                                                                               >  3:170390/1‑149 (MQ=255)
                                        gcgcTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAGGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGa                                                                                                              <  4:1390403/149‑1 (MQ=255)
                                        gcgcTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGa                                                                                                              <  1:983794/149‑1 (MQ=255)
                                          gcTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATg                                                                                                            >  1:701005/1‑149 (MQ=255)
                                              ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGa                                                                                                        <  3:1250484/149‑1 (MQ=255)
                                              ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGa                                                                                                        <  3:1318649/149‑1 (MQ=255)
                                              ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGa                                                                                                        <  1:220653/149‑1 (MQ=255)
                                              ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGa                                                                                                        <  3:1360454/149‑1 (MQ=255)
                                              ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGa                                                                                                        <  2:1461404/149‑1 (MQ=255)
                                               ctctGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGAt                                                                                                       >  4:1199204/1‑149 (MQ=255)
                                               ctctGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGAt                                                                                                       <  3:823371/149‑1 (MQ=255)
                                                tctGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATa                                                                                                      >  1:482602/1‑149 (MQ=255)
                                                tctGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATa                                                                                                      >  4:21101/1‑149 (MQ=255)
                                                  tGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAAt                                                                                                    >  3:678828/1‑149 (MQ=255)
                                                   ggCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATc                                                                                                   <  3:849499/149‑1 (MQ=255)
                                                   ggCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATc                                                                                                   <  1:475280/149‑1 (MQ=255)
                                                      gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCa                                                                                                <  2:1345295/149‑1 (MQ=255)
                                                      gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCa                                                                                                >  2:1294739/1‑149 (MQ=255)
                                                       aCAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAg                                                                                               >  4:464742/1‑149 (MQ=255)
                                                         aGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGaa                                                                                             >  2:1455891/1‑149 (MQ=255)
                                                          gTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGaaa                                                                                            <  2:711152/149‑1 (MQ=255)
                                                            gACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCa                                                                                                >  1:864827/1‑143 (MQ=255)
                                                            gACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCa                                                                                                <  2:864827/143‑1 (MQ=255)
                                                            gACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGaaa                                                                                            <  3:180217/147‑1 (MQ=255)
                                                              ccAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCa                                                                                        <  4:103741/149‑1 (MQ=255)
                                                                aGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATc                                                                                                   <  3:704141/136‑1 (MQ=255)
                                                                aGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGc                                                                                      <  3:216292/149‑1 (MQ=255)
                                                                        gAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATc                                                                              <  1:481460/149‑1 (MQ=255)
                                                                           cAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGc                                                                           >  1:1336512/1‑149 (MQ=255)
                                                                           cAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGc                                                                           <  1:347487/149‑1 (MQ=255)
                                                                             aaGATGCGGATCGCTAAGCGGGAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCtt                                                                         >  4:256451/1‑149 (MQ=255)
                                                                                  gCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTttgttg                                                                    <  2:482602/149‑1 (MQ=255)
                                                                                   cGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTgc                                                                   <  2:1336512/149‑1 (MQ=255)
                                                                                   cGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTgc                                                                   <  3:635762/149‑1 (MQ=255)
                                                                                   cGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTgc                                                                   <  2:588504/149‑1 (MQ=255)
                                                                                        cGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAAc                                                              <  1:1427471/149‑1 (MQ=255)
                                                                                        cGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAAc                                                              <  3:333080/149‑1 (MQ=255)
                                                                                          cTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCa                                                            <  4:675763/149‑1 (MQ=255)
                                                                                           tAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAg                                                           >  1:1145718/1‑149 (MQ=255)
                                                                                           tAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAg                                                           >  3:576721/1‑149 (MQ=255)
                                                                                           tAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAg                                                           <  4:1337334/149‑1 (MQ=255)
                                                                                           tAAGCGGTAAAAAGGCCAGAATGCCTTCGACATGGGAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAGTATCGGCTTTGTTGCGCAACCAg                                                           >  1:800150/1‑149 (MQ=255)
                                                                                            aaGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAgg                                                          >  2:746774/1‑149 (MQ=255)
                                                                                              gCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGaa                                                        <  1:1103479/149‑1 (MQ=255)
                                                                                                 gTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTc                                                                                                                             >  1:1133178/1‑77 (MQ=255)
                                                                                                 gTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAg                                                           >  4:352550/1‑143 (MQ=255)
                                                                                                 gTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTa                                                     >  4:197587/1‑149 (MQ=255)
                                                                                                  tAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAg                                                    >  2:1470581/1‑149 (MQ=255)
                                                                                                   aaaaaGGCCAGAATGCCTTCGCGATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAgc                                                   <  2:1241205/149‑1 (MQ=255)
                                                                                                       aGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCg                                               >  1:89916/1‑149 (MQ=255)
                                                                                                            aGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGaa                                                                                  >  3:1094027/1‑109 (MQ=255)
                                                                                                            aGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGaa                                                                                  <  4:1094027/109‑1 (MQ=255)
                                                                                                            aGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCcgc                                          >  4:378041/1‑149 (MQ=255)
                                                                                                              aaTGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCgcgc                                        >  2:1116570/1‑149 (MQ=255)
                                                                                                                   cTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATc                                   >  3:1452719/1‑149 (MQ=255)
                                                                                                                    ttCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCa                                  >  1:980745/1‑149 (MQ=255)
                                                                                                                    ttCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCa                                  <  3:407672/149‑1 (MQ=255)
                                                                                                                      cGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCAcc                                <  2:1154421/149‑1 (MQ=255)
                                                                                                                      cGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCAcc                                <  1:527268/149‑1 (MQ=255)
                                                                                                                              aaaaCCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACgcgc                        >  2:1240784/1‑149 (MQ=255)
                                                                                                                              aaaaCCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACgcgc                        <  3:609848/149‑1 (MQ=255)
                                                                                                                                 aCCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAAt                                                                                                                                 >  4:16121/1‑41 (MQ=255)
                                                                                                                                 aCCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAAt                                                                                                                                 <  3:16121/41‑1 (MQ=255)
                                                                                                                                  ccTGCTGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGtt                    >  4:845849/1‑149 (MQ=255)
                                                                                                                                  ccTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGtt                    <  1:608080/149‑1 (MQ=255)
                                                                                                                                     gCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGCCGCGCCGTTAAc                 >  1:1472760/1‑149 (MQ=255)
                                                                                                                                     gCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAAc                 >  3:1269346/1‑149 (MQ=255)
                                                                                                                                       ggaagaaCCGCATGATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAt               >  1:188328/7‑149 (MQ=255)
                                                                                                                                       ggaagaaCCGCATGATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAt               >  3:274701/7‑149 (MQ=255)
                                                                                                                                           gaaCCGCATGATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACa                >  3:224191/3‑144 (MQ=255)
                                                                                                                                           gaaCCGCATGATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAtgctg           >  1:627143/3‑149 (MQ=255)
                                                                                                                                            cACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAtgctgc          <  2:1190343/149‑1 (MQ=255)
                                                                                                                                            cACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAtgctgc          >  4:283939/1‑149 (MQ=255)
                                                                                                                                            aaCCGCATGATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACAtgctgc          >  4:632850/2‑149 (MQ=255)
                                                                                                                                            aaCCGCATGATGCGGTTCTTCCGCGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGtt                    >  3:307502/2‑139 (MQ=255)
                                                                                                                                              ccGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCCACCAGGACTTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACATGCTGCCa        >  2:1373876/1‑149 (MQ=255)
                                                                                                                                               cGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCATACGCGCCGTTAACATGCTGCCAt       <  3:1199204/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCGATTAAATTCGTCTTCACTTGCCTGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCCACCGCATCATGCGGTTCTTCCGTGCGAATGGTCGCTACCAGCGCGTGATGCTGATAATCCCAGAAAGTCAGCGGAATATCGGCTTTGTTGCGCAACCAGGAATTAGCGCCGTCCGCGCCAATCACCAGACGCGCCGTTAACATGCTGCCATCTTTC  >  NC_000913/3051575‑3051871

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