Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,227 A→G 5.3% intergenic (‑117/+196) sapA ← / ← ymjA antimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2270AG 5.3% 208.9 / 2.8 78intergenic (‑117/+196)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base A (41/31);  new base G (2/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

TCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGT  >  NC_000913/1357080‑1357372
                                                                                                                                                   |                                                                                                                                                 
tCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGAAATGAGGGCCGTGTGGTTGGTCGTTCAACAGGTGAc                                                                                                                                                  >  2:2086/1‑149 (MQ=255)
   ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                               >  2:941399/1‑149 (MQ=255)
   ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                               >  2:946056/1‑149 (MQ=255)
   ccAAAAGAGACGATAATACCTGGCGCATAATATTAACGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                               >  1:275896/1‑149 (MQ=255)
    cAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                              <  2:451913/149‑1 (MQ=255)
         gagaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTATTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACTTCGCCATc                                                                                                                                         >  1:1059438/1‑149 (MQ=255)
          agaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTGAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCt                                                                                                                                        >  1:690155/1‑149 (MQ=255)
              gATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTAGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTc                                                                                                                                    >  1:950821/1‑149 (MQ=255)
                 aataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                 >  1:990928/1‑149 (MQ=255)
                 aataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTAGCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                                                                 <  2:1165039/149‑1 (MQ=255)
                  ataCCTGGCGCATACTATTAACGGATTTTATGTATAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                                                                <  1:259311/149‑2 (MQ=255)
                     ccTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTc                                                                                                                             <  2:923957/149‑1 (MQ=255)
                      cTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                                                                            >  2:108474/1‑149 (MQ=255)
                         gcgcATAATATTAAGGGATTTTATGTAAAGTGGCTATATTACTAATATTTAATGACAATTGACATTACCGTTTTTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTAGCCATCTGTTCAACATTCGTAc                                                                                                                         <  2:782960/149‑1 (MQ=255)
                         gcgcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTATTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAc                                                                                                                         <  1:881110/149‑1 (MQ=255)
                            cATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTGCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGACGCCATCTGTTCAACATTCGTACCCg                                                                                                                      >  2:759185/1‑149 (MQ=255)
                              taataTTAAGGGATTTTATCTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTATGTTCAGGGGTCGTAATGAGGCCCGTGTGGTTGGTCTTTGAACAGGTGACGTCGCCATCTGTTAAACATTCTTACACGtt                                                                                                                    >  2:441868/1‑148 (MQ=255)
                               aataTTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                                                                   <  1:590101/149‑1 (MQ=255)
                                atatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAt                                                                                                                  >  1:200727/1‑149 (MQ=255)
                                 tatTAAGGGATTTTATTTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGCAATGAGGGCCGTGTGGTTTGTCTTTGAACAGGTGACGTCGACATCTGTTAAACATTCGTACCCGTTAtt                                                                                                                 >  2:848514/1‑149 (MQ=255)
                                    tAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACAGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                              >  1:393415/1‑149 (MQ=255)
                                      tcgtaTTTTATGTAAAGAGTCGATCTTTCTTATATTTAATCACATTTGCCATTACCGTTTGTGTTCCGGGTGCGTATTGAGGGCGGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                                                            <  2:852359/145‑1 (MQ=255)
                                          aTTTTATGTAAAGAGGCTATCTTACTAATCTTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                                        <  2:198728/149‑1 (MQ=255)
                                          aTTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                                        >  2:536577/1‑149 (MQ=255)
                                           ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCTCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                       >  2:550073/1‑149 (MQ=255)
                                           ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                       >  2:594222/1‑149 (MQ=255)
                                           ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataaa                                                                                                       >  2:269057/1‑148 (MQ=255)
                                            tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCAGTTCAACATTCGTACACGTGATTTCTctataact                                                                                                      >  2:1164994/1‑149 (MQ=255)
                                             ttATGTACAGAGGCTATCATACTAAAATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                                                                     <  1:207151/149‑1 (MQ=255)
                                                    acgagGCTATCTTACTAATATTTAATGCCATTTGCCATTACCGATTGTGGTCCGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACCGGGGACGTCGCCATCTGTTCAACATTCGTACCCGGGATTTCTCTATAACTATAACTca                                                                                              <  2:198957/147‑1 (MQ=255)
                                                    aaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTAGCGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTTGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTAAAACTca                                                                                              <  2:458943/149‑1 (MQ=255)
                                                    aaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                                                                              >  1:979745/1‑149 (MQ=255)
                                                     agatGCTATATTACTAATATTTAATGACATTTGCCATTACGGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACTTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                                                             <  1:686422/149‑1 (MQ=255)
                                                      gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcaca                                                                                            >  2:979448/1‑149 (MQ=255)
                                                         ggtATCTTACTAATATTTAATGAAATTTGCCATTAACGTTTGTGTTGAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                         <  1:441868/147‑1 (MQ=255)
                                                         gcgaTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                         <  1:361316/146‑1 (MQ=255)
                                                         gCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                         >  2:729113/1‑149 (MQ=255)
                                                            aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaa                                                                                      >  1:533941/1‑149 (MQ=255)
                                                              cTTACTAATATTTAATTACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACt                                                                                    >  2:1123525/1‑149 (MQ=255)
                                                                 aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                                  <  2:508501/148‑1 (MQ=255)
                                                                 aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                                  >  1:508501/1‑148 (MQ=255)
                                                                 aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                                 >  2:120093/1‑149 (MQ=255)
                                                                 aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                                 >  1:572308/1‑149 (MQ=255)
                                                                    aataatTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTCATGAGGGCCGTGTGGTTGGTCGTTCAACACGTGTCGTCGCCATCTGTTCAACAGTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCga                                                                              <  2:903555/149‑1 (MQ=255)
                                                                    aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCga                                                                              >  2:467369/1‑149 (MQ=255)
                                                                            aTGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCaga                                                                      >  2:312119/1‑149 (MQ=255)
                                                                              gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                                    >  2:406926/1‑149 (MQ=255)
                                                                               aCATTTGCCATTACCGTTTTTGTGCAGGGGTCTTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCTTCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                   <  1:1154335/149‑1 (MQ=255)
                                                                               aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                   >  1:180043/1‑149 (MQ=255)
                                                                               aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                   <  1:754205/149‑1 (MQ=255)
                                                                                  tttGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTATTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTa                                                                >  2:601059/1‑149 (MQ=255)
                                                                                         ttACCGTTTGTGTTCAGGGGTCGTAATGTGGGCGGTGTGGTTGGTCGTTGAACTGGTGACGTCGCCATGTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGATCAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                         <  2:954279/149‑1 (MQ=255)
                                                                                         ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                         <  2:327590/149‑1 (MQ=255)
                                                                                         ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTCAACAGGTGACGTCGCCATCTGTTCAACATACGTACCCGTGATTTCTCTATAACTATATCTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                         >  2:921344/1‑149 (MQ=255)
                                                                                              gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                            >  1:1103302/1‑125 (MQ=255)
                                                                                              gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                            <  2:1103302/125‑1 (MQ=255)
                                                                                               tttGTGTTGAGGGGTCTTAATGATGGGGGTGTGTTATGGAGTTGAACAGTTGAAGTCTCCATCTGTTCAACATTCGTTCCCGTGATTTCTCTCTAACTAAAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa                                                   <  1:820946/149‑1 (MQ=255)
                                                                                                     ttCAGGGGTCGTAATGAGGGACGTGTTGTTGGTCGTTGAACAGGTGACGTCGACATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGAGGTACGATCGTGAGAAAAACAtt                                             >  2:1084781/1‑149 (MQ=255)
                                                                                                        tggggTCTTTATGATGGCCGTGTGGGTGGGCGTTGTACAGGTGTCGTCGCCATCTGTTCAACATTCGTACCGGTGATTTCTCTATTTCTATAACTCAGAGAACAACTTTGCGAGTTGGTGAGCGGTTCGATCGTGAGAAACGGATTaaa                                          <  2:1102608/148‑1 (MQ=255)
                                                                                                        aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa                                          <  1:941399/149‑1 (MQ=255)
                                                                                                        aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa                                          <  1:946056/149‑1 (MQ=255)
                                                                                                           ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                                       >  1:436659/1‑149 (MQ=255)
                                                                                                             tCGTAATGAGCGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACg                                     <  2:275896/149‑1 (MQ=255)
                                                                                                                   tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAGCCCATTAAAGCACGCTGCAg                               >  1:34702/1‑149 (MQ=255)
                                                                                                                   tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                               <  1:729113/149‑1 (MQ=255)
                                                                                                                       ggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCAATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCg                           >  1:93643/1‑149 (MQ=255)
                                                                                                                              tggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAGTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAtt                    >  2:994224/1‑149 (MQ=255)
                                                                                                                                gttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGt                  <  1:536577/149‑1 (MQ=255)
                                                                                                                                    gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGc              <  2:341889/149‑1 (MQ=255)
                                                                                                                                     tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc             >  2:687251/1‑149 (MQ=255)
                                                                                                                                     tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc             <  1:108474/149‑1 (MQ=255)
                                                                                                                                     gcGTTTAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTACAGCACGCTGCAGTGCGTGTCATTGTTAGcc             <  1:1114520/148‑1 (MQ=255)
                                                                                                                                       gTTGAACAGATGGCGACGTCACCTGTTCAACAGTGGTCCCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAg           >  1:695675/1‑149 (MQ=255)
                                                                                                                                        ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATACCTCCCAGAACAACTTAGCGAGGAGCAGAGCGGGACGCTCGTGAGAAACCCATTAACGCACGCTGCAGTGCGTGTCATTGTTAGCCAGa          <  2:871893/149‑1 (MQ=255)
                                                                                                                                            aCAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcg      >  1:467508/1‑149 (MQ=255)
                                                                                                                                             cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcgg     >  1:141132/1‑149 (MQ=255)
                                                                                                                                                gTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGt  <  2:572308/149‑1 (MQ=255)
                                                                                                                                                gTGACGTCGCCATCTGGTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACCGAACTACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGt  <  2:1168557/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGT  >  NC_000913/1357080‑1357372

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: