Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,042 Δ1 bp 7.7% coding (287/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,0420A. 7.7% 182.5 / 10.1 65coding (287/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base . (3/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTACATTG  >  NC_000913/2324896‑2325191
                                                                                                                                                   |                                                                                                                                                     
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:256067/1‑150 (MQ=255)
                                                           tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtt                                                                                         >  2:1067143/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:1227841/150‑1 (MQ=255)
                                                                   aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTc                                                                                 <  1:311273/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)
                                                                                                                                                   aGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCAACATTg  >  2:1205955/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCTACATTG  >  NC_000913/2324896‑2325191

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