Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,236 T→C 5.2% intergenic (‑126/+187) 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,2360TC 5.2% 215.5 / 2.3 76intergenic (‑126/+187)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTAT  >  NC_000913/1357089‑1357383
                                                                                                                                                   |                                                                                                                                                   
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                                                                                             >  1:508501/1‑148 (MQ=255)
                                                        aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                                             <  2:508501/148‑1 (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                                                                                       <  2:1103302/125‑1 (MQ=255)
                                                                                     gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                                       >  1:1103302/1‑125 (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)
                                                                                                                                            gTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACg        >  2:1030214/1‑149 (MQ=255)
                                                                                                                                               gCCATCTGTTCAACATTCGTCCCCGTTTTTTGTCTATAACTATCACTCACTCACGAACTTAGCGTGGAGCAGAGCGGTACGTTCGTGAGAAACCGATTAAAGCACGATGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCtt     <  1:1091600/149‑1 (MQ=255)
                                                                                                                                                  aTCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTAt  <  1:687251/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTAT  >  NC_000913/1357089‑1357383

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: