Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 786,031 T→C 5.3% I133I (ATT→ATC)  aroG → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,0310TC 5.3% 249.5 / 2.9 75I133I (ATT→ATC) aroG3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base C (2/2);  total (34/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 = 9.94e-01

GATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGG  >  NC_000913/785885‑786179
                                                                                                                                                  |                                                                                                                                                    
gATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTaa                                                                                                                                                    >  3:92283/1‑149 (MQ=255)
       tGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGc                                                                                                                                             >  4:569215/1‑149 (MQ=255)
              cGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGc                                                                                                                                      <  3:569215/149‑1 (MQ=255)
              cGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGc                                                                                                                                      <  2:139585/149‑1 (MQ=255)
               gTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGcc                                                                                                                                     <  1:53578/149‑1 (MQ=255)
                    gcgcGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCgg                                                                                                                                <  2:357845/149‑1 (MQ=255)
                      gcgTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCa                                                                                                                              >  1:6220/1‑149 (MQ=255)
                       cgTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAg                                                                                                                             <  4:51021/149‑1 (MQ=255)
                         tCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTACATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGt                                                                                                                           <  2:555227/149‑1 (MQ=255)
                         tCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGt                                                                                                                           <  1:105742/149‑1 (MQ=255)
                           tATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGa                                                                                                                         <  1:77235/149‑1 (MQ=255)
                             ttttGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATCAACGACAGCGGTCTGCCAGCGGCAGGTGAGt                                                                                                                       >  4:61884/1‑149 (MQ=255)
                             ttttGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATCAACGACAGCGGTCTGCCAGCGGCAGGTGAGt                                                                                                                       >  2:512883/1‑149 (MQ=255)
                              tttGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGtt                                                                                                                      >  3:534105/1‑149 (MQ=255)
                                 gAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTtct                                                                                                                   >  3:316023/1‑149 (MQ=255)
                                  aaaaGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTtctc                                                                                                                  <  2:74645/149‑1 (MQ=255)
                                   aaaGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCg                                                                                                                 >  1:465660/1‑149 (MQ=255)
                                   aaaGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCg                                                                                                                 >  4:56063/1‑149 (MQ=255)
                                     aGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGat                                                                                                               <  2:587686/149‑1 (MQ=255)
                                     aGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGat                                                                                                               <  3:561835/149‑1 (MQ=255)
                                      gCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGata                                                                                                              <  2:6220/149‑1 (MQ=255)
                                      gCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGata                                                                                                              >  3:390461/1‑149 (MQ=255)
                                         gcgTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCAGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGa                                                                                                           >  3:359741/1‑149 (MQ=255)
                                                aCGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCAcccc                                                                                                    <  4:534105/149‑1 (MQ=255)
                                                         tGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATAtct                                                                                           >  4:184861/1‑149 (MQ=255)
                                                           gAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCg                                                                                         <  1:329659/149‑1 (MQ=255)
                                                                 gCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGAcc                                                                                   <  1:504247/149‑1 (MQ=255)
                                                                     aTTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCtgat                                                                               <  4:337306/149‑1 (MQ=255)
                                                                     aTTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCtgat                                                                               <  1:505456/149‑1 (MQ=255)
                                                                            aTCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTgg                                                                        <  4:457401/149‑1 (MQ=255)
                                                                               cGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGgc                                                                     <  2:594938/149‑1 (MQ=255)
                                                                                gCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGgcg                                                                    >  1:98436/1‑149 (MQ=255)
                                                                                   tatGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCaa                                                                 >  1:306841/1‑149 (MQ=255)
                                                                                      ggATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTg                                                              >  1:393085/1‑149 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGa                                                                                                                         >  4:310804/1‑89 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGa                                                                                                                         <  3:310804/89‑1 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTgg                                                             >  1:205998/1‑149 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTgg                                                             >  1:301330/1‑149 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATCAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTgg                                                             <  1:512883/149‑1 (MQ=255)
                                                                                       gATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATCAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTgg                                                             <  3:61884/149‑1 (MQ=255)
                                                                                        ataataGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGgc                                                            <  4:370005/149‑1 (MQ=255)
                                                                                        ataataGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGgc                                                            >  1:572657/1‑149 (MQ=255)
                                                                                          aataGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGgcgc                                                          <  4:447312/149‑1 (MQ=255)
                                                                                                ttCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTac                                                    <  2:465660/149‑1 (MQ=255)
                                                                                                     gATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCg                                               >  2:71244/1‑149 (MQ=255)
                                                                                                      aTCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGa                                              >  3:336648/1‑149 (MQ=255)
                                                                                                         aaCGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTacca                                                  >  1:54006/1‑142 (MQ=255)
                                                                                                         aaCGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTacca                                                  <  2:54006/142‑1 (MQ=255)
                                                                                                           cgacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAACCGc                                         <  1:513329/149‑1 (MQ=255)
                                                                                                            gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGata                                                                                                              >  1:546371/1‑79 (MQ=255)
                                                                                                            gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGata                                                                                                              <  2:546371/79‑1 (MQ=255)
                                                                                                            gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCtga                                                                                >  4:541693/1‑109 (MQ=255)
                                                                                                            gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCtga                                                                                <  3:541693/109‑1 (MQ=255)
                                                                                                            gacgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCa                                        <  2:205998/149‑1 (MQ=255)
                                                                                                             acgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAg                                       <  1:177688/149‑1 (MQ=255)
                                                                                                              cgGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAgg                                      <  3:184861/149‑1 (MQ=255)
                                                                                                                      gTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACcgc                              >  3:108515/1‑149 (MQ=255)
                                                                                                                       tataGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACcgcg                             <  4:92283/149‑1 (MQ=255)
                                                                                                                        ataGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGa                            >  3:402401/1‑149 (MQ=255)
                                                                                                                            cccGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTg                        <  4:390461/149‑1 (MQ=255)
                                                                                                                            cccGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTg                        <  3:336512/149‑1 (MQ=255)
                                                                                                                            cccGTAAATCGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTg                        >  3:548130/1‑149 (MQ=255)
                                                                                                                             ccGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTgg                       >  4:279614/1‑149 (MQ=255)
                                                                                                                                  aaTTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCa                  >  1:560830/1‑149 (MQ=255)
                                                                                                                                  aaTTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCa                  <  4:530979/149‑1 (MQ=255)
                                                                                                                                   aTTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAg                 >  3:258285/1‑149 (MQ=255)
                                                                                                                                    ttGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAgg                <  4:108515/149‑1 (MQ=255)
                                                                                                                                    ttGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAgg                <  2:393085/149‑1 (MQ=255)
                                                                                                                                    ttGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAgg                >  4:507141/1‑149 (MQ=255)
                                                                                                                                     tgctgcTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAggg               <  4:316023/149‑1 (MQ=255)
                                                                                                                                        tgctTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCtt            >  2:379122/1‑149 (MQ=255)
                                                                                                                                          ctTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTc          >  4:125152/1‑149 (MQ=255)
                                                                                                                                           ttGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCt         <  3:125152/149‑1 (MQ=255)
                                                                                                                                            tGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCtt        <  4:402401/149‑1 (MQ=255)
                                                                                                                                             gATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTg       >  1:404511/1‑149 (MQ=255)
                                                                                                                                                 ttAACGACAGCGGTCTGCCGGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCg   >  3:344771/1‑149 (MQ=255)
                                                                                                                                                 ttAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCg   >  2:22092/1‑149 (MQ=255)
                                                                                                                                                  tAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTAGCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCgg  >  2:490370/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATAATAGCTTCCAGATCAACGACGGTCTGCGTATAGCCCGTAAATTGCTGCTTGATATTAACGACAGCGGTCTGCCAGCGGCAGGTGAGTTTCTCGATATGATCACCCCACAATATCTCGCTGACCTGATGAGCTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGG  >  NC_000913/785885‑786179

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