Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,119,997 T→C 5.4% A54A (GCA→GCG wzxC ← putative colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,119,9970TC 5.4% 232.0 / 2.6 74A54A (GCA→GCGwzxCputative colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base C (2/2);  total (42/32)
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.06e-01

GCCGATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGCGCCGCTGAT  >  NC_000913/2119850‑2120137
                                                                                                                                                   |                                                                                                                                            
gCCGATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCt                                                                                                                                              <  1:24342/148‑1 (MQ=255)
gCCGATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCt                                                                                                                                              <  5:317444/148‑1 (MQ=255)
 ccGATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGc                                                                                                                                            <  8:249064/149‑1 (MQ=255)
   gATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATGCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAAAAAGCGTATCTGCCa                                                                                                                                          <  1:340626/149‑1 (MQ=255)
         aTCACTCAACAAAAACACCGCCACGCCCACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGCTGCCGACGTCAGAAAGCGTATCTGCCAGCGCGa                                                                                                                                    >  4:82224/1‑149 (MQ=255)
         aTCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGATCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGa                                                                                                                                    >  5:3744/1‑149 (MQ=255)
           cACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCCGCCAgc                                                                                                                                        <  4:71083/143‑1 (MQ=255)
           cACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCCGCCAgc                                                                                                                                        >  3:71083/1‑143 (MQ=255)
            aCTCAACAAAAACACCGCCACGCACACCATGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATa                                                                                                                                  <  7:308507/148‑1 (MQ=255)
             cTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAAt                                                                                                                                >  1:229277/1‑149 (MQ=255)
              tCAACAAAAACACCGCCACGCACCCCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATc                                                                                                                               >  4:354978/1‑149 (MQ=255)
                aacaaAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCAc                                                                                                                             >  8:122567/1‑149 (MQ=255)
                 acaaAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGACACTGCTTTCTTTTCGCTGAATAATCGACTTAGCGATACCGAAGTCAGAAAGCGAATCTGCCAGCGCGATAATCAcc                                                                                                                            >  1:232837/1‑149 (MQ=255)
                  caaAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCa                                                                                                                           <  8:128499/149‑1 (MQ=255)
                   aaaaaCACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAg                                                                                                                          >  2:74063/1‑149 (MQ=255)
                   aaaaaCACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAg                                                                                                                          <  1:74063/149‑1 (MQ=255)
                   aaaaaCACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAg                                                                                                                          >  8:62115/1‑149 (MQ=255)
                                   acacCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAgg                                                                                                          <  7:291661/149‑1 (MQ=255)
                                     accacGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCa                                                                                                                                          >  6:224406/1‑115 (MQ=255)
                                     accacGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCa                                                                                                                                          <  5:224406/115‑1 (MQ=255)
                                       cacGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCg                                                                                                       >  7:333130/1‑148 (MQ=255)
                                         cGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAAc                                                                                                    <  7:122567/149‑1 (MQ=255)
                                            tCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt                                                                                                                 <  6:131027/133‑1 (MQ=255)
                                            tCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt                                                                                                                 >  5:131027/1‑133 (MQ=255)
                                            tCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtgg                                                                                                 >  8:482086/1‑149 (MQ=255)
                                            tCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtgg                                                                                                 >  4:25900/1‑149 (MQ=255)
                                               ccagccagACGATCAGCCAGTACAACTTGTAGATTTCAAGGTGACTAATTTATTTTAGCTGAATAATAGAGATAGCGATACCGAACTCAGAAAGAGTATCTCCCACCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtcgtag                                                                                              >  8:247567/1‑147 (MQ=255)
                                               ccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtggtgg                                                                                              <  8:97812/149‑1 (MQ=255)
                                                   cccGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGt                                                                                          <  3:192463/149‑1 (MQ=255)
                                                    cctaCGTTCAGCCAGTACAACGTTGTGAGTTCAAGGTGACTTAGATCTTGTCGCGGAATAATCGAGTGAGGGATACCGAACTCCGAAAGCGTATCTGCCAGCGCTATAATCACCCGAGACACGGTAAGCAGGCCGAACTGGTGGTTGTc                                                                                         <  2:232837/146‑1 (MQ=255)
                                                    ccGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTc                                                                                         <  7:62115/149‑1 (MQ=255)
                                                        cGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATa                                                                                     <  3:354978/149‑1 (MQ=255)
                                                             aGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCTTATCTGCCAGCGCGATAATCACAAGCGACACGGTAAGAAGGCCGAACTGGTGGTTGTCGATAAAccc                                                                                >  8:448974/1‑148 (MQ=255)
                                                              gCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTTGTGGTTGTCGATAATCcgc                                                                               <  8:333130/149‑1 (MQ=255)
                                                              gCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcgc                                                                               <  7:482086/149‑1 (MQ=255)
                                                              gCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcgc                                                                               >  4:47341/1‑149 (MQ=255)
                                                              gCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGAAACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGCCCGAACTGGTGGTTGTCGATAATCcgc                                                                               >  2:243436/1‑149 (MQ=255)
                                                               ccAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAAc                                                                                                    >  1:312283/1‑127 (MQ=255)
                                                               ccAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAAc                                                                                                    <  2:312283/127‑1 (MQ=255)
                                                               ccAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcgcg                                                                              <  5:458545/149‑1 (MQ=255)
                                                                      aaCGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCAc                                                                       >  3:425834/1‑149 (MQ=255)
                                                                         gtggtgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGt                                                                    >  2:104029/1‑149 (MQ=255)
                                                                            gtgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCAt                                                                 <  2:229277/149‑1 (MQ=255)
                                                                            gtgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCAt                                                                 >  3:248110/1‑149 (MQ=255)
                                                                             tgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATc                                                                <  3:25900/149‑1 (MQ=255)
                                                                                gTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTg                                                              >  3:459021/1‑148 (MQ=255)
                                                                                gTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCCACCAGcgcgta                                                                                                                                   >  6:314986/1‑77 (MQ=255)
                                                                                gTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCCACCAGcgcgta                                                                                                                                   <  5:314986/79‑3 (MQ=255)
                                                                                    aaGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCAc                                                           >  8:170185/1‑147 (MQ=255)
                                                                                    aaGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCa                                                         <  6:3744/149‑1 (MQ=255)
                                                                                            tGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtg                                                                                                  <  4:327487/100‑1 (MQ=255)
                                                                                            tGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtg                                                                                                  >  3:327487/1‑100 (MQ=255)
                                                                                              aTTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCCGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAgg                                                 <  8:174175/147‑1 (MQ=255)
                                                                                               tttCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACgg                                                                     <  2:377314/126‑1 (MQ=255)
                                                                                               tttCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACgg                                                                     >  1:377314/1‑126 (MQ=255)
                                                                                                        gCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCAc                                                           <  4:67832/127‑1 (MQ=255)
                                                                                                        gCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCAc                                                           >  3:67832/1‑127 (MQ=255)
                                                                                                        gCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCgatgatgat                                      >  7:197821/1‑148 (MQ=255)
                                                                                                          tGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGAGGATCAc                                   >  2:418236/1‑149 (MQ=255)
                                                                                                          tGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCAAGCGCGATAATCACCAGCGACACGGTAAGCAGTACGAACTGGTGGTTGTCGATAAGCCGAGCCAGCACGGTCATCTGCCCCAGCCCCAGGCCgaggatgatagc                                   >  5:477249/1‑146 (MQ=255)
                                                                                                             ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGt                                <  4:425834/149‑1 (MQ=255)
                                                                                                              taatCGAGTTAGCGATACAGAAGTCAGAAAGCTTATCTGACAGCGCGATAATCACCTGCGACACGTTAAGCAGGCCGAACTGGTGTTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACAAGCCCGAGGCCGATGATGATCACCGTg                               <  8:161914/149‑1 (MQ=255)
                                                                                                                    aGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATAACCGTGGCAATc                         >  5:383906/1‑149 (MQ=255)
                                                                                                                      ttAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGc                       >  2:23234/1‑149 (MQ=255)
                                                                                                                        aGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGcc                      <  1:104029/148‑1 (MQ=255)
                                                                                                                         gCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGa                    >  1:311323/1‑149 (MQ=255)
                                                                                                                         gCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAAGCGCCGa                    >  6:135810/1‑149 (MQ=255)
                                                                                                                             tACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCAc                <  2:311323/149‑1 (MQ=255)
                                                                                                                              aCCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACt               >  7:73174/1‑149 (MQ=255)
                                                                                                                                 gAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTcg            >  3:215749/1‑149 (MQ=255)
                                                                                                                                     tCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGcgcc        >  4:417562/1‑149 (MQ=255)
                                                                                                                                       aGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGcgccgc      >  7:73011/1‑149 (MQ=255)
                                                                                                                                       aGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGcgccgc      >  6:187133/1‑149 (MQ=255)
                                                                                                                                       aGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGcgccgc      >  1:32476/1‑149 (MQ=255)
                                                                                                                                        gAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGCGCCGCt     >  3:416210/1‑149 (MQ=255)
                                                                                                                                           aGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGCGCCGCTGAt  <  7:247567/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCCGATGAGATCACTCAACAAAAACACCGCCACGCACACCACGATCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGCACCAGCCCGAGGCCGATGATGATCACCGTGGCAATCGCCGACCACTTCGCGCCGCTGAT  >  NC_000913/2119850‑2120137

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