Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,234,967 T→C 6.4% G10G (GGT→GGC)  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,9670TC 6.4% 237.9 / 5.0 78G10G (GGT→GGC) fadRfatty acid metabolism regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

AAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGAGATGGC  >  NC_000913/1234819‑1235114
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                            >  6:448714/1‑128 (MQ=255)
                                              tCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGaaa                                                                                                                            <  5:448714/128‑1 (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                                                                                  >  6:160907/1‑149 (MQ=255)
                                                                   ccTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCg                                                                                  >  3:343118/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                                                        >  8:378250/1‑149 (MQ=255)
                                                                                             ttttGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGaattaatt                                                        <  7:378250/149‑1 (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                       >  4:281100/1‑149 (MQ=255)
                                                                                                                              ttAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTa                       >  1:42504/1‑149 (MQ=255)
                                                                                                                               tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt                           >  1:430389/6‑144 (MQ=255)
                                                                                                                               tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt                           <  2:430389/139‑1 (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)
                                                                                                                                            ccGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCAc          >  7:238236/1‑148 (MQ=255)
                                                                                                                                             cggcggGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACga        <  7:409441/149‑1 (MQ=255)
                                                                                                                                              ggcggGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCAc          >  2:109026/1‑146 (MQ=255)
                                                                                                                                                 gggTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGAGATg    <  1:228832/149‑1 (MQ=255)
                                                                                                                                                   gTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGAGATGGc  >  5:47662/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGAGATGGC  >  NC_000913/1234819‑1235114

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