Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,804,528 C→A 11.1% intergenic (+67/‑287) nrdF → / → proV ribonucleoside‑diphosphate reductase 2, beta subunit, ferritin‑like protein/glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,804,5280CA11.1% 198.8 / 4.3 74intergenic (+67/‑287)nrdF/proVribonucleoside‑diphosphate reductase 2, beta subunit, ferritin‑like protein/glycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (40/24);  new base A (2/6);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCGG  >  NC_000913/2804380‑2804675
                                                                                                                                                    |                                                                                                                                                     
tCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATc                                                                                                                                                       <  8:102332/149‑1 (MQ=255)
tCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATc                                                                                                                                                       <  1:179357/149‑1 (MQ=255)
    gATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACg                                                                                                                                                    >  4:362220/1‑148 (MQ=255)
     aTTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAATCGGTTGAAACAGAAGATGAAGACTGGAATTTATGAGGTTTTTATTTTCAAAAATATCACTACACGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACga                                                                                                                                                   >  3:46132/1‑147 (MQ=255)
      tttCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCa                                                                                                                                                  >  8:110034/1‑148 (MQ=255)
         cTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATAACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaat                                                                                                                                              >  4:32548/1‑148 (MQ=255)
          ttttCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaataa                                                                                                                                            >  8:236821/1‑149 (MQ=255)
          ttttCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaataa                                                                                                                                            >  2:249617/1‑149 (MQ=255)
            ttCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaata                                                                                                                                             <  6:232876/146‑1 (MQ=255)
               cGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgt                                                                                                                                      <  3:169160/149‑1 (MQ=255)
                     aGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGTAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGAt                                                                                                                                <  3:369500/149‑1 (MQ=255)
                     aGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGAt                                                                                                                                <  7:107736/149‑1 (MQ=255)
                       gCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCAAGCAAATAATTTGTGGTGATc                                                                                                                               >  5:25275/1‑148 (MQ=255)
                        ctcctcTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTa                                                                                                                             <  4:65201/149‑1 (MQ=255)
                             ctTATGAGAGGGGGAAAGCGGTTGCAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCCCAGCAGGGAAATAATTCCCGCCAAATTTTTTTTTTTATGTTAATTATATTTTGGTGATCTACACTg                                                                                                                        <  4:82460/149‑1 (MQ=255)
                              ttATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGa                                                                                                                       >  6:237082/1‑149 (MQ=255)
                                atgggctGGGGAAAGCGGTTGAAACACAACCTGCAGCCTGGAATTTCTGAGGGTGTTAATTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATTTTTTTTTTTATGTTAATTATTTTTTGGTGATCTACACTGa                                                                                                                       <  2:255531/141‑1 (MQ=255)
                                aTGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGa                                                                                                                       <  3:364495/147‑1 (MQ=255)
                                aTGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGAt                                                                                                                      <  7:236821/148‑1 (MQ=255)
                                   tgATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATActc                                                                                                                  >  1:163023/1‑149 (MQ=255)
                                    gATGGCGAAAGCGGGTGAAACAGAAGATGCCGACTGGAAGTTCTGCGGGTGGTAATTTCAAAAATATCACTACCCGCAGCATGGAAATAATACCCGCCAAATTTTTTTTTTTATGTTAATTATTTTTTGGAGATCTACACTGATActct                                                                                                                 <  4:224583/149‑1 (MQ=255)
                                     atcgggAAAGCGCTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATAACTACCCGCAGCAGGGAAATAATTACCGCCAAATAGATATTTATCACGCAAATAATTTGTGGTGATCTAAACTGATACTCTg                                                                                                                >  8:156795/4‑149 (MQ=255)
                                        gggAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGtt                                                                                                              <  5:58682/148‑1 (MQ=255)
                                          gAAAGCGGTGGAAACAGAAGATGATGACTGGAATTTCTGAGTCTGTTATGATCAAAAATATCACTACCCGCAGCAGGGAAATAAATCCCGCCTAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCa                                                                                                           <  7:194734/149‑1 (MQ=255)
                                           aaaGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTTTTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgtg                                                                                                                                     >  3:204462/1‑122 (MQ=255)
                                           aaaGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgtg                                                                                                                                     <  4:204462/122‑1 (MQ=255)
                                                 gTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATc                                                                                                                                                       >  5:49784/1‑100 (MQ=255)
                                                 gTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATc                                                                                                                                                       <  6:49784/100‑1 (MQ=255)
                                                  ttGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATTTTTTTTTTTATGTTAATTATATTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGc                                                                                                   <  2:190366/149‑1 (MQ=255)
                                                  gtGAAACAGAAGATGAAGACTCGAATGTCTGACGGTGTTATTTTCAAACATATCACTACCCGCAGCAGGGAAATAATTCCAGCCATTTTTTTTTTTTTATGTTAATTATTTTTTGGAGATCTACACTGATACTCTGTTGCATTATTCGc                                                                                                   <  8:192455/148‑1 (MQ=255)
                                                        cAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAATATATCACTACCCGCAGAATGTTAATAATTGCCGCCATATTTTTTTTTTTATTTTAATTATTTTTTTGAGATCTACACTGATACTCTGTTGCATTATTCGCCTGaaa                                                                                             <  6:211202/149‑1 (MQ=255)
                                                         agaagaTGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAAc                                                                                            >  8:171565/1‑149 (MQ=255)
                                                          gaagaTGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGTTGATCTACACTGATAc                                                                                                                    >  2:230340/1‑124 (MQ=255)
                                                          gaagaTGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATAc                                                                                                                    <  1:230340/124‑1 (MQ=255)
                                                          gaagaTGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAAGAGCTTTTTATAACGCAAATAATTTGTGGTGATCTACACTGATACTCTTTTGAATTATTCGCATGAAAc                                                                                            >  8:364215/1‑148 (MQ=255)
                                                               tGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACaa                                                                                       >  4:215025/1‑148 (MQ=255)
                                                                gAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATAACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAata                                                                                     >  1:162349/1‑149 (MQ=255)
                                                                 aaGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAatat                                                                                    >  4:78578/1‑149 (MQ=255)
                                                                   gACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTc                                                                                  <  5:237082/149‑1 (MQ=255)
                                                                    aCTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATAtt                                                                                   >  4:226966/1‑147 (MQ=255)
                                                                              cTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGttttt                                                                        <  7:171565/148‑1 (MQ=255)
                                                                                   ggTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTa                                                                  >  4:146663/1‑149 (MQ=255)
                                                                                   ggTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTa                                                                  >  7:237300/1‑149 (MQ=255)
                                                                                      gtTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCt                                                               >  7:266377/1‑149 (MQ=255)
                                                                                      gtTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCt                                                               >  3:323308/1‑149 (MQ=255)
                                                                                       ttATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgtg                                                                                                                                     >  4:285189/1‑78 (MQ=255)
                                                                                       ttATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgtg                                                                                                                                     <  3:285189/78‑1 (MQ=255)
                                                                                          ttttCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCAGTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGa                                                           >  5:114959/1‑149 (MQ=255)
                                                                                             tCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACa                                                         >  8:339020/1‑148 (MQ=255)
                                                                                                          tACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACtt                                           >  4:104844/1‑149 (MQ=255)
                                                                                                           aCCCGCAGCAGGGAAATAAATCCCGCCATATATTTTTTTATCTGCAAATTATTTTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACttt                                          <  6:254078/149‑1 (MQ=255)
                                                                                                            cccGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTc                                                                      <  2:391809/120‑1 (MQ=255)
                                                                                                            cccGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTc                                                                      >  1:391809/1‑120 (MQ=255)
                                                                                                             ccGCAGCATGGAAATAATTCCCGCCAAATATTTTTTTATATGTTAAATAATTTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTtct                                        <  4:46132/149‑1 (MQ=255)
                                                                                                               gcagcaGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAAGTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTtctc                                       >  5:396232/1‑148 (MQ=255)
                                                                                                                cagcagGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGa                                     >  3:136864/1‑149 (MQ=255)
                                                                                                                   cagGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGAttt                                  <  7:110034/149‑1 (MQ=255)
                                                                                                                     gggAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTtct                                        <  3:135426/141‑1 (MQ=255)
                                                                                                                     gggAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTtct                                        >  4:135426/1‑141 (MQ=255)
                                                                                                                      ggAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGc                                >  6:284651/1‑148 (MQ=255)
                                                                                                                      ggAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTCTGGTGATCTACACTGATACTCTGTTGCATTATTCACCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGc                                >  2:359495/1‑148 (MQ=255)
                                                                                                                         aataatTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGa                                     <  5:393406/140‑1 (MQ=255)
                                                                                                                         aataatTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGa                                     >  6:393406/1‑140 (MQ=255)
                                                                                                                          ataatTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCa                           <  3:104844/149‑1 (MQ=255)
                                                                                                                                cccGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCt                      <  3:78578/148‑1 (MQ=255)
                                                                                                                                 ccGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTta                    >  5:127434/1‑149 (MQ=255)
                                                                                                                                  cGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTta                    >  4:9040/1‑148 (MQ=255)
                                                                                                                                      aaaTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCa               >  2:105706/1‑149 (MQ=255)
                                                                                                                                       aaTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCTCCTGAAACCACAATATTCAGGCGATTTTTCGCAATCTTTAACAAAAAATATCAACTTTCTCGATTTGCTCTCATCCCTTATATAAc              >  6:33153/1‑149 (MQ=255)
                                                                                                                                       aaTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCAc              >  6:218183/1‑149 (MQ=255)
                                                                                                                                       aaTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCAc              <  8:266377/149‑1 (MQ=255)
                                                                                                                                        aTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACg             <  3:226966/149‑1 (MQ=255)
                                                                                                                                        aTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACg             >  4:320219/1‑149 (MQ=255)
                                                                                                                                           gCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACggg           >  4:210927/1‑148 (MQ=255)
                                                                                                                                            ttttttATAACTTAAATATTTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTtata                  <  2:44373/139‑1 (MQ=255)
                                                                                                                                            cTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTtata                  >  1:44373/1‑140 (MQ=255)
                                                                                                                                              ttttATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTATCCCTATCATTGACAAAAAATATCAACTTTCTCGATTTAATCTAATCCCATATATAACTGGAAAt       >  8:416496/1‑149 (MQ=255)
                                                                                                                                                   tCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCgg  <  3:210927/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCGG  >  NC_000913/2804380‑2804675

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