Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,555 T→C 5.1% Y168Y (TAT→TAC frlR → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,5550TC 5.1% 260.1 / 2.8 78Y168Y (TAT→TACfrlRputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (42/32);  new base C (2/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ATCCCCGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATCA  >  NC_000913/3504415‑3504701
                                                                                                                                            |                                                                                                                                                  
aTCCCCGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTAtt                                                                                                                                            <  4:572261/149‑1 (MQ=255)
   cccGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCa                                                                                                                                         >  3:731850/1‑149 (MQ=255)
   cccGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCa                                                                                                                                         >  3:208032/1‑149 (MQ=255)
    ccGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCTTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAg                                                                                                                                        <  1:697546/149‑1 (MQ=255)
                  gtgtTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCCGGAGCGTTTTGacac                                                                                                                          <  4:592856/149‑1 (MQ=255)
                     ttCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGt                                                                                                                                  <  2:571424/138‑1 (MQ=255)
                     ttCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGt                                                                                                                                  >  1:571424/1‑138 (MQ=255)
                     ttCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTACCAGTTATTTCAGGAGCGTTTTGACAcgcg                                                                                                                       >  2:725712/1‑149 (MQ=255)
                     ttCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTACCAGTTATTTCAGGAGCGTTTTGACAcgcg                                                                                                                       >  4:417117/1‑149 (MQ=255)
                         aTCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTg                                                                                                                   >  1:231903/1‑149 (MQ=255)
                          tctcTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTgg                                                                                                                  >  4:557053/1‑149 (MQ=255)
                              tGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAg                                                                                                              >  2:495605/1‑149 (MQ=255)
                               gCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGc                                                                                                             >  4:346801/1‑149 (MQ=255)
                                       aTGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACaaaaa                                                                                                     <  2:149960/149‑1 (MQ=255)
                                         gTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGa                                                                                                   >  1:495131/1‑149 (MQ=255)
                                         gTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGa                                                                                                   <  3:287143/149‑1 (MQ=255)
                                            tctcGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCa                                                                                                >  3:614417/1‑149 (MQ=255)
                                            tctcGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCa                                                                                                >  4:794107/1‑149 (MQ=255)
                                              tcGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATc                                                                                              <  1:519234/149‑1 (MQ=255)
                                              tcGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATc                                                                                              >  3:112740/1‑149 (MQ=255)
                                                gATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGt                                                                                                                    <  3:139745/125‑1 (MQ=255)
                                                gATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGt                                                                                                                    >  4:139745/1‑125 (MQ=255)
                                                gATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGa                                                                                            >  3:801965/1‑149 (MQ=255)
                                                gATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGa                                                                                            <  2:754701/149‑1 (MQ=255)
                                                 aTAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGat                                                                                           >  1:241959/1‑149 (MQ=255)
                                                       agCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCttt                                                                                     <  1:495605/149‑1 (MQ=255)
                                                         ccGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTgc                                                                                   <  4:423414/149‑1 (MQ=255)
                                                               tttATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCAc                                                                             <  4:536032/149‑1 (MQ=255)
                                                               tttATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCAc                                                                             <  3:740635/149‑1 (MQ=255)
                                                                   ttGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCAcccgc                                                                         >  1:553007/1‑149 (MQ=255)
                                                                     gATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCAcccgccc                                                                       >  4:412903/1‑149 (MQ=255)
                                                                             aTGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCaa                                                               >  2:711168/1‑149 (MQ=255)
                                                                             aTGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCaa                                                               >  2:785872/1‑149 (MQ=255)
                                                                             aTGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCaa                                                               >  1:682807/1‑149 (MQ=255)
                                                                               ggATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGa                                                                                                   >  1:29510/1‑111 (MQ=255)
                                                                               ggATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGa                                                                                                   <  2:29510/111‑1 (MQ=255)
                                                                                gATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAAt                                                            >  3:437532/1‑149 (MQ=255)
                                                                                 aTCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATg                                                           >  4:170627/1‑149 (MQ=255)
                                                                                    ccGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCt                                                        >  2:716664/1‑149 (MQ=255)
                                                                                    ccGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCt                                                        >  4:226940/1‑149 (MQ=255)
                                                                                    ccGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTACCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCt                                                        <  3:417117/149‑1 (MQ=255)
                                                                                    ccGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTACCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCt                                                        <  1:725712/149‑1 (MQ=255)
                                                                                     cGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTg                                                       >  2:54364/1‑149 (MQ=255)
                                                                                     cGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTg                                                       <  2:231903/149‑1 (MQ=255)
                                                                                         gTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAAt                                                   <  2:553007/149‑1 (MQ=255)
                                                                                         gTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAAt                                                   <  2:241959/149‑1 (MQ=255)
                                                                                             cgTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAGAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGa                                               <  1:711168/149‑1 (MQ=255)
                                                                                                       aCTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGaa                                     <  4:112740/149‑1 (MQ=255)
                                                                                                        cTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGaaa                                                    >  3:700668/1‑133 (MQ=255)
                                                                                                        cTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGaaa                                                    <  4:700668/133‑1 (MQ=255)
                                                                                                        cTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAAc                                    >  2:244395/1‑149 (MQ=255)
                                                                                                        cTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAAc                                    >  2:300278/1‑149 (MQ=255)
                                                                                                            gatgaGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgtt                                >  4:260424/1‑149 (MQ=255)
                                                                                                            gatgaGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgtt                                >  2:434249/1‑149 (MQ=255)
                                                                                                             atgaGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgttg                               >  1:239618/1‑149 (MQ=255)
                                                                                                              tgaGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgttgt                              <  3:170627/149‑1 (MQ=255)
                                                                                                                 gaTTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTc                           <  3:412903/149‑1 (MQ=255)
                                                                                                                 gaTTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTc                           <  2:239618/149‑1 (MQ=255)
                                                                                                                  aTTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCg                          >  2:216163/1‑149 (MQ=255)
                                                                                                                   tttACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgtt                                >  4:115171/1‑142 (MQ=255)
                                                                                                                   tttACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCgtt                                <  3:115171/142‑1 (MQ=255)
                                                                                                                       cGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCa                     <  3:794107/149‑1 (MQ=255)
                                                                                                                         tCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGc                   >  4:704677/1‑149 (MQ=255)
                                                                                                                           gaaggaagCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCaa                 >  3:271855/1‑149 (MQ=255)
                                                                                                                           gaaggaagCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCaa                 >  4:400829/1‑149 (MQ=255)
                                                                                                                            aaggaagCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCaaa                >  3:244310/1‑149 (MQ=255)
                                                                                                                                  gCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGcc          >  3:591308/1‑149 (MQ=255)
                                                                                                                                  gCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGcc          <  1:716664/149‑1 (MQ=255)
                                                                                                                                  gCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGcc          <  3:346801/149‑1 (MQ=255)
                                                                                                                                   cTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCt         <  1:785872/149‑1 (MQ=255)
                                                                                                                                    tCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCtt        >  2:536801/1‑149 (MQ=255)
                                                                                                                                      cACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGAGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGATGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAACATCGCCTTTg      <  1:421561/149‑1 (MQ=255)
                                                                                                                                       aCCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGa     <  4:801965/149‑1 (MQ=255)
                                                                                                                                        ccTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGAt    >  1:69859/1‑149 (MQ=255)
                                                                                                                                         cTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGGCCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATc   >  1:253158/1‑149 (MQ=255)
                                                                                                                                         cTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATc   >  4:721007/1‑149 (MQ=255)
                                                                                                                                          tATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATCa  <  4:437532/149‑1 (MQ=255)
                                                                                                                                          tATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATCa  <  4:614417/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
ATCCCCGGAAACAGCGAAGTGTTCCATCTCTGCCGGGTGATGTATCTCGATAAAGAGCCGCTGTTTATTGATAGTTCATGGATCCCGCTGTCGCGTTATCCTGACTTTGATGAGATTTACGTCGAAGGAAGCTCCACCTATCAGTTATTTCAGGAGCGTTTTGACACGCGAGTGGTCAGCGACAAAAAGACCATCGATATCTTTGCCGCCACCCGCCCGCAGGCAAAATGGCTGAAATGCGAACTGGGCGAACCGTTGTTTCGCATCAGCAAAATCGCCTTTGATCA  >  NC_000913/3504415‑3504701

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