Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,234,962 G→C 6.5% A9P (GCG→CCG)  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,9620GC 6.5% 255.8 / 8.7 77A9P (GCG→CCG) fadRfatty acid metabolism regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base C (2/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGA  >  NC_000913/1234815‑1235108
                                                                                                                                                   |                                                                                                                                                  
ggggAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCcggc                                                                                                                                                   <  1:132008/149‑1 (MQ=255)
ggggAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCcggc                                                                                                                                                   <  2:48404/149‑1 (MQ=255)
   gAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCggg                                                                                                                                                <  4:214846/149‑1 (MQ=255)
    aaaGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGt                                                                                                                                               <  1:368615/149‑1 (MQ=255)
    aaaGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGt                                                                                                                                               <  1:91164/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                                                                            >  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                 >  1:42504/1‑149 (MQ=255)
                                                                                                                                  ttAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTa                 >  4:281100/1‑149 (MQ=255)
                                                                                                                                   tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt                     >  1:430389/6‑144 (MQ=255)
                                                                                                                                   tcttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAgt                     <  2:430389/139‑1 (MQ=255)
                                                                                                                                     ttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                            >  4:402475/4‑135 (MQ=255)
                                                                                                                                     ttccgcgAAACCCGCCGGGCTTTGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACtacgttac                            <  3:402475/132‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTCTGGTATGATGAGTCCAACTTTGTTTTGCTGTGTTATGGAAATCTCACTATGGTCATTAAGGCGCAAAGCCCGGCGGGTTTCGCGGAAGAGTACATTATTGAAAGTATCTGGAATAACCGCTTCCCTCCCGGGACTATTTTGCCCGCAGAACGTGAACTTTCAGAATTAATTGGCGTAACGCGTACTACGTTACGTGAAGTGTTACAGCGTCTGGCACGA  >  NC_000913/1234815‑1235108

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