Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,234,957 G→A 6.3% S7N (AGC→AAC)  fadR → fatty acid metabolism regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,234,9570GA 6.3% 249.5 / 4.4 79S7N (AGC→AAC) fadRfatty acid metabolism regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (2/3);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTG  >  NC_000913/1234811‑1235102
                                                                                                                                                  |                                                                                                                                                 
tCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCGACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGcc                                                                                                                                                 >  8:343881/1‑149 (MQ=255)
    ggggAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCcggc                                                                                                                                             <  1:132008/149‑1 (MQ=255)
    ggggAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCcggc                                                                                                                                             <  2:48404/149‑1 (MQ=255)
       gAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCggg                                                                                                                                          <  4:214846/149‑1 (MQ=255)
        aaaGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGt                                                                                                                                         <  1:91164/149‑1 (MQ=255)
        aaaGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGt                                                                                                                                         <  1:368615/149‑1 (MQ=255)
            agaTGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCGACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTcg                                                                                                                                     <  7:343881/149‑1 (MQ=255)
                gCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCgg                                                                                                                                  >  7:197987/1‑148 (MQ=255)
                 cGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGa                                                                                                                                 <  2:308798/148‑1 (MQ=255)
                 cGGAAAATGAAGCCGTGATCCCTTTTTCTTCTTTTTGTCAGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGaa                                                                                                                                <  6:180541/149‑1 (MQ=255)
                  ggAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCgg                                                                                                                                  >  2:437832/1‑146 (MQ=255)
                   gAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAaga                                                                                                                              <  5:169550/149‑1 (MQ=255)
                   gAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCAATCAGCGTAGTTAGCCCTCTGGTATGAGGAGGCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAaga                                                                                                                              <  5:285461/149‑1 (MQ=255)
                     aaaTGAAGCCTTGATCCCTTTTTCTTCTTTTTGACTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGt                                                                                                                            <  1:98141/149‑1 (MQ=255)
                      aaTGAATCCTTGATCCCTTTTTCATCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATTGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTa                                                                                                                           <  6:106345/149‑1 (MQ=255)
                          aaGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACat                                                                                                                        >  1:152959/1‑148 (MQ=255)
                                   tCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTTGTATGATGAGGCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGt                                                                                                              <  5:50716/149‑1 (MQ=255)
                                     ccTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTTAAAGTAt                                                                                                            >  5:374701/1‑149 (MQ=255)
                                     ccTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTAt                                                                                                            >  1:363975/1‑149 (MQ=255)
                                       tttttCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCt                                                                                                          >  4:47945/1‑149 (MQ=255)
                                              tcttTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATaa                                                                                                   >  1:354593/1‑149 (MQ=255)
                                                tttttGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAAcc                                                                                                 >  4:300257/1‑149 (MQ=255)
                                                   ttGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGc                                                                                               <  8:197987/148‑1 (MQ=255)
                                                    tGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGGCCAACTTTGTTTTGCTGTGGTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTGGCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCt                                                                                              <  4:342827/148‑1 (MQ=255)
                                                      tCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGaaa                                                                                                                <  5:448714/128‑1 (MQ=255)
                                                      tCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGaaa                                                                                                                >  6:448714/1‑128 (MQ=255)
                                                                gCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGAc                                                                                 <  3:78050/149‑1 (MQ=255)
                                                                gCGTAGTTAGCCCTATGGTATGATGAGTCCAACTTTGTGTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATATGGAATAACCGCTTCCCTCCCGGGAc                                                                                 >  2:441391/1‑149 (MQ=255)
                                                                     gTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTAttt                                                                            <  1:437832/149‑1 (MQ=255)
                                                                      ttAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTAtttt                                                                           <  1:96950/149‑1 (MQ=255)
                                                                         gCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGaa                                                                                                      >  6:148915/1‑119 (MQ=255)
                                                                         gCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGaa                                                                                                      <  5:148915/119‑1 (MQ=255)
                                                                          cccTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGccc                                                                       >  1:246869/1‑149 (MQ=255)
                                                                           ccTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCg                                                                      >  3:343118/1‑149 (MQ=255)
                                                                           ccTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCg                                                                      >  6:160907/1‑149 (MQ=255)
                                                                            ctctGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCg                                                                      >  2:87819/1‑148 (MQ=255)
                                                                                 gTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAAc                                                                >  2:42013/1‑149 (MQ=255)
                                                                                 gTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAAc                                                                >  8:233991/1‑149 (MQ=255)
                                                                                    tgatgaGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAAc                                                                >  3:185689/1‑146 (MQ=255)
                                                                                    tgatgaGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAAc                                                                <  4:185689/146‑1 (MQ=255)
                                                                                    tgatgaGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGt                                                              <  3:300257/148‑1 (MQ=255)
                                                                                       tgaGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAAc                                                          <  5:160907/149‑1 (MQ=255)
                                                                                         aGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACtt                                                        >  8:409441/1‑149 (MQ=255)
                                                                                           tCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACttt                                                       >  6:277818/1‑148 (MQ=255)
                                                                                            ccAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCa                                                     <  6:154812/149‑1 (MQ=255)
                                                                                              aaCTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGAGTATATTGCCCGCAGAACGAGAACTTTAAg                                                    >  8:164530/1‑148 (MQ=255)
                                                                                                     ttttGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGaattaatt                                            <  7:378250/149‑1 (MQ=255)
                                                                                                     ttttGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGaattaatt                                            >  8:378250/1‑149 (MQ=255)
                                                                                                                 aTGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAAcgcg                                 >  5:208643/1‑148 (MQ=255)
                                                                                                                  tGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGta                               <  2:363975/149‑1 (MQ=255)
                                                                                                                  tGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGt                                <  4:343118/148‑1 (MQ=255)
                                                                                                                   ggAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCa                                                     >  2:345535/1‑126 (MQ=255)
                                                                                                                       aTCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTA‑TTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgt                         >  2:272029/1‑149 (MQ=255)
                                                                                                                        tctcACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgt                         >  5:48781/1‑149 (MQ=255)
                                                                                                                          tcACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCTGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCTCAGAACGTGAACTTTCAGAATTAATTGGCGTAAc                                    <  8:397955/136‑1 (MQ=255)
                                                                                                                          tcACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAAc                                    >  7:397955/1‑136 (MQ=255)
                                                                                                                            aCTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttacg                     >  2:228832/1‑149 (MQ=255)
                                                                                                                            aCTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgtt                        <  6:48781/146‑1 (MQ=255)
                                                                                                                            aCTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTA‑TTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                      <  1:272029/147‑1 (MQ=255)
                                                                                                                              tATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCTGAAGAGTACATTATTGAAAGTATCTGGAATAACAGCTTACCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTg                   <  7:164530/149‑1 (MQ=255)
                                                                                                                              tATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttacgt                    >  2:321812/1‑148 (MQ=255)
                                                                                                                                tGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGa                  <  1:87819/148‑1 (MQ=255)
                                                                                                                                tGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTATGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGaa                 <  6:374701/149‑1 (MQ=255)
                                                                                                                                  gTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt               <  1:441391/149‑1 (MQ=255)
                                                                                                                                   tCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgtg              >  1:121943/1‑149 (MQ=255)
                                                                                                                                      ttAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTa           >  1:42504/1‑149 (MQ=255)
                                                                                                                                      ttAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTa           >  4:281100/1‑149 (MQ=255)
                                                                                                                                       tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt               <  2:430389/139‑1 (MQ=255)
                                                                                                                                       tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt               >  1:430389/6‑144 (MQ=255)
                                                                                                                                         ttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                      <  3:402475/132‑1 (MQ=255)
                                                                                                                                         ttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                      >  4:402475/4‑135 (MQ=255)
                                                                                                                                         aGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACa         <  2:152959/148‑1 (MQ=255)
                                                                                                                                         aGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACa         <  7:233991/148‑1 (MQ=255)
                                                                                                                                           ccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTAc          <  2:436028/144‑1 (MQ=255)
                                                                                                                                            cgcAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGt     <  2:246869/149‑1 (MQ=255)
                                                                                                                                               aaaGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                      <  3:184615/129‑1 (MQ=255)
                                                                                                                                               aaaGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                      >  4:184615/1‑129 (MQ=255)
                                                                                                                                               aaaGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTg  <  2:42504/149‑1 (MQ=255)
                                                                                                                                               aaaGCCAGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTAGGTGAAGTGTTACAGCGGCTg  >  1:146304/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTG  >  NC_000913/1234811‑1235102

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