Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,804,528 C→A 11.9% 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.9% 148.3 / 5.6 69intergenic (+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 (38/22);  new base A (2/7);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 2.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCGG  >  NC_000913/2804381‑2804675
                                                                                                                                                   |                                                                                                                                                   
cACGATTTCTTTTCCGGTTCAGGCTGCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCa                                                                                                                                                    <  2:1445822/149‑1 (MQ=255)
   gATTTCTTTTCCGGTTCAGGCTCCTCTTAAGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGGGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGc                                                                                                                                                 <  1:972962/149‑1 (MQ=255)
   gATTTCTTTTCCGGTTCAGGCGCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGTAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGc                                                                                                                                                 <  2:1178902/149‑1 (MQ=255)
    aTTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCa                                                                                                                                                >  2:760315/1‑149 (MQ=255)
        cTTGTCCGGTTCAGGCTCCTGTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaata                                                                                                                                            >  1:1515558/1‑149 (MQ=255)
         ttttCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGGAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAaataa                                                                                                                                           <  2:373030/149‑1 (MQ=255)
               ggTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTtgtg                                                                                                                                     <  1:566907/149‑1 (MQ=255)
                     ggCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATc                                                                                                                               <  1:239298/149‑1 (MQ=255)
                     ggCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACAGCAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGaaa                                                                                                                               >  1:679327/1‑147 (MQ=255)
                          ctctTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCAAGAAAATAATTTGTCGTGATCTacac                                                                                                                          >  2:93556/1‑149 (MQ=255)
                                tgtgATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTTATAc                                                                                                                    >  1:933613/1‑149 (MQ=255)
                                 gtgtTGGCGAAAGCGGATGAAACAGCAGATGAAGACTGGAATTTCTGAGGGTGGTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATTTTTTTTTTTATGTTAATTATTTTTTGGTGATCTACACTGATAct                                                                                                                   <  2:368042/149‑1 (MQ=255)
                                 gtgATGGGGAAAGATGTTGAAACAGAAGATGAAGACTGGAATTTATGAGTGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGTAATAATTGCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATAct                                                                                                                   <  1:1152356/149‑1 (MQ=255)
                                  tgATGGGGAAAGCGGTTGAAGCAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATActc                                                                                                                  >  1:23144/1‑149 (MQ=255)
                                     tGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGt                                                                                                               >  2:1363555/1‑149 (MQ=255)
                                      ggggAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGtt                                                                                                              >  2:902645/1‑149 (MQ=255)
                                         gAAAGCGGTTGAAACAGAAGATGAAGACTTGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCAGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCa                                                                                                           >  1:53242/1‑149 (MQ=255)
                                         gAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCa                                                                                                           >  2:1343791/1‑149 (MQ=255)
                                         gAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTAACACGCAAATAATAAGTGGTGATCTACACTGATACTCTGTTGCa                                                                                                           >  1:47177/1‑149 (MQ=255)
                                          aaaGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCat                                                                                                          >  2:1281731/1‑149 (MQ=255)
                                           aaGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATACCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCatt                                                                                                         <  1:610849/149‑1 (MQ=255)
                                             gCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCACCAAATAGCTTTTTATCACGCAAATAATTTGTGTTGATCTACACTGATACTCTGTTGCattat                                                                                                       >  1:1497736/1‑149 (MQ=255)
                                                 ttGAAACATAAGATCAATACTTGAATTTCTGAGGGTTTTATTTTCAAAAAAATCACTACCAGCAGCAGGGAAATAATTACCGACAAATATTTTTTTATATGTTAATTATTTTTTCGTGATCTACACTGATACTCTGTTGCATTATTCGc                                                                                                   <  2:906328/149‑1 (MQ=255)
                                                  tGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGcc                                                                                                  <  1:138247/149‑1 (MQ=255)
                                                   gAAACAGAAGATGAAGACTGGAATTTCTGAGTGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCt                                                                                                 <  1:1219461/149‑1 (MQ=255)
                                                    aaaCAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTg                                                                                                >  1:1333689/1‑149 (MQ=255)
                                                    aaaCAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTAGCCTg                                                                                                >  1:820019/1‑149 (MQ=255)
                                                      aCAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTCCCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCCCGCAAATAATTTGTGGTGATCTACCCTGATACTCTGTTGCATTATTCGCCTGaa                                                                                              >  1:1081482/1‑149 (MQ=255)
                                                      aCAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTACAAAAATATCACTACCCGCAGCAGGGAAATAAATCCCGCCAAATAGCTTTTTATAACGCAAATAATTTGTGGTCATCTACACTGCTACTCTGTTGCATTATTCGCCTGaa                                                                                              >  1:425448/1‑149 (MQ=255)
                                                        agaagaTGAAGTCTGGAATTTCTGAGGGTTTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAAc                                                                                            <  1:308025/149‑1 (MQ=255)
                                                             tgcaaGACTGGTTTTTCTGACGCTGTTATATTCAAAAATATCACTAGCGGCAGCAGCGTAATAATTCCCTCCATATTTTTTTTTTTATGTTAATTATTTTTTGGTGATCTAGACTGATACTCTGTTGCATTATTCGCCTGAAACCACaa                                                                                       <  2:1561915/146‑1 (MQ=255)
                                                               gAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGCGAAATAATTCCCGCCAAATAGCATTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCATGAAACCAAAata                                                                                     >  1:1521376/1‑149 (MQ=255)
                                                                aaGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAatat                                                                                    >  2:717315/1‑149 (MQ=255)
                                                                     tGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCATGAAACCACAATATTCAgg                                                                               >  1:511805/1‑149 (MQ=255)
                                                                      ggAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAATTAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACGAGAATATTCAGGc                                                                              >  2:1321483/1‑149 (MQ=255)
                                                                      ggAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAATTAATTTGTGTTGATATACACTGATACTCTGTTGCATTATTCTCCTGAAACCACAATATTCAAGc                                                                              >  2:1038393/1‑149 (MQ=255)
                                                                            tCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGAt                                                                                                                                >  2:858736/1‑93 (MQ=255)
                                                                            tCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGAt                                                                                                                                <  1:858736/93‑1 (MQ=255)
                                                                                 gggTGTTATTTTAAAAAATATCACTACCCTCTGCAGGGAAATAATTCCCGCCTAATTTTTTTTTTTAAGTTAATTTTTTTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCt                                                                   <  2:679327/149‑1 (MQ=255)
                                                                                           ttCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAATTAGCTTTTTATCACGCAAATAATTTGTGGTTATCTACACTGATACTCTGTTGCATTATTCACATGAAACAACAATATTCAGTCGTTTATTCGCTATCtttttca                                                         >  1:1185139/1‑145 (MQ=255)
                                                                                            tCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATAACGCAAATAATTTGTGCTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACaa                                                        >  1:1482729/1‑149 (MQ=255)
                                                                                              aaaaaTATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCTCATGAAACCAAAATATTCAGGCGTTTTTTCGCTATCTTTGACaaaa                                                      >  1:317213/1‑149 (MQ=255)
                                                                                               aaaaTATCACTACCCGCAGCAGGGAAATAATTACCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACaaaaa                                                     <  1:97240/149‑1 (MQ=255)
                                                                                                   tatCACTACCCTCAGCAGGGAAATAATTCCCGTCAAATATCTTTTTATCACGCATATAATTTGTGGTTATCTACAATGATACTCTGATGCATTATTCGCCTGAATCCACAATATTCATGCGTTTTTTCGCTATCTTTGACAAATAatat                                                 <  1:102274/149‑1 (MQ=255)
                                                                                                    atCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATc                                                >  2:54033/1‑149 (MQ=255)
                                                                                                     tCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTTTATCGAAATTATTTTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCa                                               <  2:920851/149‑1 (MQ=255)
                                                                                                     tCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCa                                               <  2:1333689/149‑1 (MQ=255)
                                                                                                      cACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTATTTATCACGCAAATAATTTGAGGTGATCTACACTGATACTCTGTTGCATTATTAGACTGAAACCACAATATTCAGGCTTTTATTCGCTATCTTTGACAATAAATATCaa                                              >  2:767434/1‑149 (MQ=255)
                                                                                                      cACTACCCGCAGCAGGGAAATAATTCCAGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCaa                                              <  1:1023513/149‑1 (MQ=255)
                                                                                                         tACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCCCTATCTTTGACAAAAAATATCAACtt                                           >  1:544127/1‑149 (MQ=255)
                                                                                                            ccGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTtct                                        <  1:524830/149‑1 (MQ=255)
                                                                                                              gcagcaGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGATATCTTTTACAAAAAATATCAACTTactcc                                      >  2:669584/1‑148 (MQ=255)
                                                                                                                agcagGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGAt                                    >  1:130637/1‑149 (MQ=255)
                                                                                                                  cagGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGAAACTCTGTTGCATTATTCGCCTGAAACCACAAAATTCAGGCGGTTTTTCGCTATCTTTGACAAAAAAAATCAACTTTCTCGAttt                                  >  1:1237166/1‑149 (MQ=255)
                                                                                                                   agGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTg                                 <  1:503014/149‑1 (MQ=255)
                                                                                                                   agGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTg                                 >  2:800538/1‑149 (MQ=255)
                                                                                                                      gAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGctc                              <  1:717315/149‑1 (MQ=255)
                                                                                                                             ttCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTTGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGccc                       <  1:760315/149‑1 (MQ=255)
                                                                                                                                ccGCAAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTATTCGCTATCTTTCACAAAAAATATCAACTTTCTCGATTTGCTCTCAGACCtaa                    >  2:871744/1‑147 (MQ=255)
                                                                                                                                 cGCCAAATAGTTTTTTATAACGAAATTATTTTTTGGTGATCTACACTGATACTCTGTTGCATTATTAGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTtat                   <  2:276276/149‑1 (MQ=255)
                                                                                                                                 cGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTtat                   >  1:248154/1‑149 (MQ=255)
                                                                                                                                   ccAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTatat                 >  1:818413/1‑149 (MQ=255)
                                                                                                                                     ataTAGCTTTTTTTCTCGGAAATATTTTTTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCa               <  2:154925/147‑1 (MQ=255)
                                                                                                                                       aTAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGTCCTTTTTTCGCTATCTTTGAAAAAAAATATCAACTTTCTCCATTTGCACTCAACCCTTATATAACg             >  2:918351/1‑149 (MQ=255)
                                                                                                                                        tAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACgg            >  1:1111954/1‑149 (MQ=255)
                                                                                                                                           cTTTTTATCACGCAAATAATTTGTGGTTATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGaa         <  1:902645/149‑1 (MQ=255)
                                                                                                                                             ttttATCCCGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATAAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAAt       >  2:284412/1‑149 (MQ=255)
                                                                                                                                              tttATAACTTAAATAATTTGTGGTGATATACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAAtt      <  2:1321566/149‑1 (MQ=255)
                                                                                                                                                  tCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCgg  >  2:1554760/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACGATTTCTTTTCCGGTTCAGGCTCCTCTTATGTGATGGGGAAAGCGGTTGAAACAGAAGATGAAGACTGGAATTTCTGAGGGTGTTATTTTCAAAAATATCACTACCCGCAGCAGGGAAATAATTCCCGCCAAATAGCTTTTTATCACGCAAATAATTTGTGGTGATCTACACTGATACTCTGTTGCATTATTCGCCTGAAACCACAATATTCAGGCGTTTTTTCGCTATCTTTGACAAAAAATATCAACTTTCTCGATTTGCTCTCAGCCCTTATATCACGGGAAATTCCGG  >  NC_000913/2804381‑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: