Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,037:1 +T 8.1% coding (282/1323 nt) atoE → short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,325,0371.T 8.1% 226.4 / 16.5 62coding (282/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base . (29/28);  new base T (3/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

GCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCT  >  NC_000913/2324889‑2325185
                                                                                                                                                     |                                                                                                                                                    
gCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑g                                                                                                                                                     <  1:7457/150‑1 (MQ=255)
     gTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTAccg                                                                                                                                                >  1:1164866/1‑148 (MQ=255)
       gAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAggg                                                                                                                                              <  1:406150/150‑1 (MQ=255)
           aTGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTGGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTc                                                                                                                                          <  2:43160/150‑1 (MQ=255)
                gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCCTCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCt                                                                                                                                     <  1:555096/150‑1 (MQ=255)
                gggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTctttgc                                                                                                                                     <  2:913399/150‑6 (MQ=255)
                  gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTctttgcgg                                                                                                                                   <  2:216858/150‑8 (MQ=255)
                  gTGACGGTTTCTGGAACTTGCTGGCGGTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTgg                                                                                                                                   <  1:163843/150‑1 (MQ=255)
                      cGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCAc                                                                                                                               >  2:342435/1‑150 (MQ=255)
                       ggTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACt                                                                                                                              >  1:611390/1‑150 (MQ=255)
                           tCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttct                                                                                                                          >  2:1132816/1‑150 (MQ=255)
                            cTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttctt                                                                                                                         >  2:1037771/1‑150 (MQ=255)
                             tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                        <  1:687036/150‑1 (MQ=255)
                               gAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCgg                                                                                                                      <  1:783235/150‑1 (MQ=255)
                                aagttGCTGGCGTTTGGTAGGCAGATGGCGCTTATCATGGTTACCGGTCATGCCCTTGCCAGCTCTGCTCGGGTGATAAGTTTGCTGCGTACTGCCGCCTGCGCCGTAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGt                                                                                                                     <  2:76706/147‑1 (MQ=255)
                                      cTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTc                                                                                                               <  2:1289760/150‑1 (MQ=255)
                                      cTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTc                                                                                                               >  2:1082004/1‑150 (MQ=255)
                                               ggTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTGGTGtc                                                                                                      >  1:810312/1‑150 (MQ=255)
                                                  aTGCAGATGGCGCTTATCATCGTTACCGGTCCTGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatc                                                                                                   <  2:545111/150‑1 (MQ=255)
                                                   tGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatca                                                                                                  <  2:426942/150‑1 (MQ=255)
                                                              cTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCGCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAttt                                                                                       <  1:1061637/150‑1 (MQ=255)
                                                                  tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtt                                                                                   >  2:1067143/1‑150 (MQ=255)
                                                                  tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtt                                                                                   >  2:256067/1‑150 (MQ=255)
                                                                   catcGTTACCGGTCATGCCCTTGCCAGCTCTGCGCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                  >  1:707575/1‑150 (MQ=255)
                                                                    aacGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTAAGTCGCTTGTGTCATCAACTGGGGAtttggtatt                                                                                 >  1:477686/3‑147 (MQ=255)
                                                                      cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATtttgtttggt                                                                               <  1:1037316/150‑1 (MQ=255)
                                                                      cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATttggtttggt                                                                               <  2:357246/150‑1 (MQ=255)
                                                                      cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGGGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATttggtttggt                                                                               <  2:334907/150‑1 (MQ=255)
                                                                        ttACCGGTCAGGCCCTTGCCAGCTCTGCGCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTAATGCTGGTCACTTTCTTCGGTGCAGTCGCTTGGGTCATCAACTGGGGATTTGGTTTGGTTg                                                                             <  2:96940/150‑1 (MQ=255)
                                                                          aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                           <  1:311273/150‑1 (MQ=255)
                                                                          aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                           <  1:1227841/150‑1 (MQ=255)
                                                                           ccGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAGTTGGTGTGGTTGTCg                                                                          >  1:634548/1‑150 (MQ=255)
                                                                           ccGGTCATGCCCTTGCCAGCGCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCg                                                                          <  2:263415/150‑1 (MQ=255)
                                                                            cGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCGCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCATCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCgg                                                                         <  1:538940/150‑1 (MQ=255)
                                                                               gcATGCCCTTGCCAGCGCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgcgc                                                                      <  1:1113222/149‑1 (MQ=255)
                                                                                cATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCa                                                                     >  1:1242676/1‑150 (MQ=255)
                                                                                cATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCa                                                                     <  2:10963/150‑1 (MQ=255)
                                                                                  tGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAAt                                                                   >  1:707735/1‑150 (MQ=255)
                                                                                    cccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGt                                                                 >  1:800399/1‑150 (MQ=255)
                                                                                     ccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTGGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGtt                                                                >  1:32256/1‑150 (MQ=255)
                                                                                     ccTGGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGtt                                                                <  1:563464/150‑1 (MQ=255)
                                                                                       tggCCAGCTCTGCTCCGTTGAAAAGGTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTg                                                              <  1:1082868/148‑1 (MQ=255)
                                                                                        tGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGc                                                             <  2:1211235/150‑1 (MQ=255)
                                                                                           cAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCg                                                          >  2:222109/1‑150 (MQ=255)
                                                                                                  gCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTc                                                   <  1:121642/150‑1 (MQ=255)
                                                                                                   cTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCg                                                  >  1:507338/1‑150 (MQ=255)
                                                                                                        gTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCgg                                             >  2:189982/1‑150 (MQ=255)
                                                                                                            aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACTGGGCGTCATGCTGGTCACTTTCTTCGGTTGAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                         >  2:650988/1‑150 (MQ=255)
                                                                                                            aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCAGTGAAGTCGCCCGGCGAg                                         >  2:1128613/1‑150 (MQ=255)
                                                                                                                    ctgcGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATAAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGt                                 >  1:947808/1‑150 (MQ=255)
                                                                                                                      gcGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGGTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCGCCGGttg                               >  1:491358/1‑149 (MQ=255)
                                                                                                                            tGCCGCCTCCGCCGCAAAGACGCCC‑GTCCAGGGCGTCATGCTGGTCACTTTCGTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCATTGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTa                         >  1:186201/1‑150 (MQ=255)
                                                                                                                                       ccgcAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCa              >  2:1032092/1‑150 (MQ=255)
                                                                                                                                       ccgcAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGCCGAGTCCCCGGTTCTGATTATCCGTTGCTCa              >  1:800481/1‑150 (MQ=255)
                                                                                                                                          cAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTg           <  2:707575/150‑1 (MQ=255)
                                                                                                                                           catgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc          >  1:439407/4‑150 (MQ=255)
                                                                                                                                           catgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc          >  1:667013/4‑150 (MQ=255)
                                                                                                                                            atgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc         >  2:346505/3‑150 (MQ=255)
                                                                                                                                            aaGACGCCC‑GTACAGGGCGGCATGCTGGTCACGTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGcc         <  1:36156/150‑1 (MQ=255)
                                                                                                                                              gACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTg       <  1:1132816/150‑1 (MQ=255)
                                                                                                                                               aCGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgc      >  2:388435/1‑150 (MQ=255)
                                                                                                                                                cGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg     >  1:385802/1‑150 (MQ=255)
                                                                                                                                                 gCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCACGGTTCTGATTATCCGTTGCTCATTGCCTgcgc    >  1:672316/1‑150 (MQ=255)
                                                                                                                                                  ccc‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc   >  1:996017/1‑150 (MQ=255)
                                                                                                                                                  ccc‑GTACAGGGCGTCATGCTGGGCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGCTGTTGGCGAAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCACCGGTTCTGATTATCCGTTGCTCATTGCCTGCGcc   <  1:557344/150‑1 (MQ=255)
                                                                                                                                                   cc‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCt  >  2:1287165/1‑150 (MQ=255)
                                                                                                                                                   cc‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGAGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCt  >  1:701736/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCT  >  NC_000913/2324889‑2325185

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