Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 839,105 C→T 6.5% A35T (GCA→ACA)  ybiX ← Fe(II)‑dependent oxygenase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913839,1050CT 6.5% 2325.3 / 93.5 719A35T (GCA→ACA) ybiXFe(II)‑dependent oxygenase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (336/336);  new base T (27/20);  total (363/356)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 744 total aligned reads displayed.

AACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGA  >  NC_000913/838957‑839253
                                                                                                                                                    |                                                                                                                                                    
 aCAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTgcg                                                                                                                                                     >  2:4860408/1‑149 (MQ=255)
 aCAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTgcg                                                                                                                                                     >  1:4657149/1‑149 (MQ=255)
    gcggcgTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCg                                                                                                                                                  <  3:1768476/149‑1 (MQ=255)
     cggcgTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCgg                                                                                                                                                 <  4:4544488/149‑1 (MQ=255)
         gTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCggtggt                                                                                                                                             >  2:5108793/1‑149 (MQ=255)
          tGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGgtggtg                                                                                                                                            <  3:4049532/149‑1 (MQ=255)
            gAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGAc                                                                                                                                          >  4:6730944/1‑149 (MQ=255)
                     cgcgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGtccatc                                                                                                                                 <  2:4037157/149‑1 (MQ=255)
                      gcgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGtccatcc                                                                                                                                >  1:1080320/1‑149 (MQ=255)
                         gCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCAcc                                                                                                                             >  2:720526/1‑149 (MQ=255)
                              gCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTc                                                                                                                        >  2:1271419/1‑149 (MQ=255)
                                 gcgGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGc                                                                                                                     <  1:4143789/149‑1 (MQ=255)
                                  cgGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCt                                                                                                                    <  4:306587/149‑1 (MQ=255)
                                     cAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGt                                                                                                                 <  3:1500615/149‑1 (MQ=255)
                                      aaaGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGtt                                                                                                                <  1:1269203/149‑1 (MQ=255)
                                      aaaGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGtt                                                                                                                >  3:5236617/1‑149 (MQ=255)
                                      aaaGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGtt                                                                                                                >  4:5109209/1‑149 (MQ=255)
                                      aaaGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGtt                                                                                                                <  1:27691/149‑1 (MQ=255)
                                       aaGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTc                                                                                                               <  3:2782329/149‑1 (MQ=255)
                                         gAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCa                                                                                                             >  4:6297062/1‑149 (MQ=255)
                                         gAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCa                                                                                                             <  3:1371883/149‑1 (MQ=255)
                                         gAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCa                                                                                                             <  4:3485019/149‑1 (MQ=255)
                                          aaTAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAg                                                                                                            >  2:4215384/1‑149 (MQ=255)
                                          aaTAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAg                                                                                                            <  4:4228543/149‑1 (MQ=255)
                                              aaGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAgttgt                                                                                                        >  4:68947/1‑149 (MQ=255)
                                                gCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGTGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTc                                                                                                      >  3:987231/1‑149 (MQ=255)
                                                       ggTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGaaaa                                                                                               <  1:6871656/149‑1 (MQ=255)
                                                       ggTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGaaaa                                                                                               <  4:2025295/149‑1 (MQ=255)
                                                            aCCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGc                                                                                          <  4:4800092/149‑1 (MQ=255)
                                                            aCCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGc                                                                                          <  2:6865037/149‑1 (MQ=255)
                                                             ccGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTg                                                                                                                  >  1:5171726/1‑124 (MQ=255)
                                                             ccGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCg                                                                                         >  2:3599670/1‑149 (MQ=255)
                                                              cgcgTTCAGTACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTAGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCTAGTTGTTCGCGAAAACGAGCGa                                                                                        >  1:6240229/1‑149 (MQ=255)
                                                              cgcgTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGa                                                                                        >  4:2882371/1‑149 (MQ=255)
                                                              cgcgTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGa                                                                                        >  3:3408140/1‑149 (MQ=255)
                                                                  ttCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTc                                                                                    <  2:5995819/149‑1 (MQ=255)
                                                                    cAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAAc                                                                                              >  2:1532010/1‑137 (MQ=255)
                                                                       cACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCtgt                                                                                <  1:3048426/148‑1 (MQ=255)
                                                                         ccTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGc                                                                             >  1:1485262/1‑149 (MQ=255)
                                                                          cTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCg                                                                            >  3:1034472/1‑149 (MQ=255)
                                                                          cTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCg                                                                            >  2:1070124/1‑149 (MQ=255)
                                                                          cTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCg                                                                            >  2:214650/1‑149 (MQ=255)
                                                                          cTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACATGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGAGGTGACGCGTCCATACACCCATTCGGCTTGATCCAGTTGTTCGCGAAAACGAGCGACGTCCAGTGGCg                                                                            >  4:5929735/1‑149 (MQ=255)
                                                                            cATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGAt                                                                          <  1:1604286/149‑1 (MQ=255)
                                                                            cATTTTGCAGGGCGGCGTATAACGTGCTGCAGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGAt                                                                          <  3:6271586/149‑1 (MQ=255)
                                                                              ttttGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATaa                                                                        >  3:3608145/1‑149 (MQ=255)
                                                                              ttttGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATaa                                                                        <  3:6608392/149‑1 (MQ=255)
                                                                              ttttGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATaa                                                                        >  3:3166865/1‑149 (MQ=255)
                                                                                  gCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACg                                                                    >  3:4890312/1‑149 (MQ=255)
                                                                                  gCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACg                                                                    >  1:4901544/1‑149 (MQ=255)
                                                                                     gggCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCg                                                                 <  1:4429018/149‑1 (MQ=255)
                                                                                     gggCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCg                                                                 >  1:341504/1‑149 (MQ=255)
                                                                                     gggCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCg                                                                 <  2:6475040/149‑1 (MQ=255)
                                                                                      ggcggcGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCgg                                                                <  4:6504343/149‑1 (MQ=255)
                                                                                            gTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAAtgt                                                          <  2:2961351/149‑1 (MQ=255)
                                                                                             tataACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAAtgtg                                                         <  3:5635184/149‑1 (MQ=255)
                                                                                                aaCGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGAc                                                                                       >  3:6068830/1‑116 (MQ=255)
                                                                                                aaCGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTa                                                      >  4:2210344/1‑149 (MQ=255)
                                                                                                  cGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTAca                                                    <  2:4657149/149‑1 (MQ=255)
                                                                                                   gTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGTGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTAcat                                                   <  4:6113334/149‑1 (MQ=255)
                                                                                                     gctgcGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTAcatca                                                 >  2:322797/1‑149 (MQ=255)
                                                                                                      ctgcGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTAcatcat                                                >  4:2700276/1‑149 (MQ=255)
                                                                                                       tgcGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATa                                               <  3:6954383/149‑1 (MQ=255)
                                                                                                        gcGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATAc                                              <  2:4489310/149‑1 (MQ=255)
                                                                                                           ggTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAAt                                           >  4:248544/1‑149 (MQ=255)
                                                                                                             tgtCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCgg                                                                >  4:5528097/1‑126 (MQ=255)
                                                                                                              gtCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCTTTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATACCTCCCCGGTAACGTGGTACATCACTCATTTtc                                        <  4:4018003/149‑1 (MQ=255)
                                                                                                               tCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTtct                                       <  1:6649342/149‑1 (MQ=255)
                                                                                                                 gACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTcc                                     <  2:662735/149‑1 (MQ=255)
                                                                                                                  aCCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGa                                                                                            <  3:2903842/93‑1 (MQ=255)
                                                                                                                  aCCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCgg                                                                >  3:447105/1‑121 (MQ=255)
                                                                                                                   ccTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCaa                                   >  4:761396/1‑149 (MQ=255)
                                                                                                                      gTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCaaaaa                                <  4:3291811/149‑1 (MQ=255)
                                                                                                                       ttGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCt                                                                                  <  3:2344731/98‑1 (MQ=255)
                                                                                                                       ttGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCg                                                                 <  1:3905919/115‑1 (MQ=255)
                                                                                                                       ttGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAg                               >  1:5674776/1‑149 (MQ=255)
                                                                                                                         gATTGTTCTTAACTTGCGCACCTTGTGTGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTg                             >  3:5528000/1‑149 (MQ=255)
                                                                                                                         gATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTg                             >  3:2518503/1‑149 (MQ=255)
                                                                                                                         gATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTg                             <  2:6184406/149‑1 (MQ=255)
                                                                                                                          aTTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTgg                            <  4:5784711/149‑1 (MQ=255)
                                                                                                                          aTTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTgg                            >  3:1953526/1‑149 (MQ=255)
                                                                                                                          aTTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTgg                            >  2:3831555/1‑149 (MQ=255)
                                                                                                                           ttGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTggg                           <  1:3651038/149‑1 (MQ=255)
                                                                                                                           ttGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTggg                           <  1:603065/149‑1 (MQ=255)
                                                                                                                              ttCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGcc                        <  1:1841353/149‑1 (MQ=255)
                                                                                                                                 ttAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCt                                                                                                                    >  4:185580/1‑54 (MQ=255)
                                                                                                                                 ttAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGAt                                                                          >  1:5368597/1‑96 (MQ=255)
                                                                                                                                 ttAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGAt                                                                          <  2:5368597/96‑1 (MQ=255)
                                                                                                                                  tAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTgcg                    >  3:894145/1‑149 (MQ=255)
                                                                                                                                     cTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATa                                                                         <  1:2466779/93‑1 (MQ=255)
                                                                                                                                     cTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATa                                                                         >  2:2466779/1‑93 (MQ=255)
                                                                                                                                     cTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGccc                 >  2:5676703/1‑149 (MQ=255)
                                                                                                                                     cTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGccc                 <  1:6394973/149‑1 (MQ=255)
                                                                                                                                      ttGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGcccc                >  4:1258260/1‑149 (MQ=255)
                                                                                                                                       tGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCca               >  1:802231/1‑149 (MQ=255)
                                                                                                                                          gcACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACAt            >  4:3701117/1‑149 (MQ=255)
                                                                                                                                          gcACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACAt            >  1:2217243/1‑149 (MQ=255)
                                                                                                                                            aCCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCt          <  3:2210344/149‑1 (MQ=255)
                                                                                                                                              cTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGa        >  2:2828561/1‑149 (MQ=255)
                                                                                                                                               ttGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGaa       <  2:3396813/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCATTCGGCTTGTTCCAGTTGTTCGCGAAAACGAGCGACGTCCTGTGGCGATAACACGCCGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGA  >  NC_000913/838957‑839253

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: ‑