Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 779,661 T→G 5.9% W22G (TGG→GGG)  ybgF → periplasmic TolA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913779,6610TG 5.9% 302.1 / 8.8 85W22G (TGG→GGG) ybgFperiplasmic TolA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (43/37);  new base G (2/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTG  >  NC_000913/779514‑779808
                                                                                                                                                   |                                                                                                                                                   
aaaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCt                                                                                                                                                     >  2:343290/1‑148 (MQ=255)
 aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTg                                                                                                                                                    <  2:529866/148‑1 (MQ=255)
 aaGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAACGGCCCCCTg                                                                                                                                                    >  1:243438/1‑148 (MQ=255)
  aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTgg                                                                                                                                                   >  1:258085/1‑148 (MQ=255)
  aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGACATACTCCAAAAACCGTCGTGCGGTACTGCTTTACTAAGAGAATTACATGAGCAGTAACTTCAGACATCAACTATTGAGTATGTCGTTACTGGTTGGTATAGCGGCCCCCTggg                                                                                                                                                  >  4:806508/1‑149 (MQ=255)
   gAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                 >  4:469806/1‑149 (MQ=255)
      aaaCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCg                                                                                                                                               <  2:637183/148‑1 (MQ=255)
            gCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTg                                                                                                                                         >  2:379866/1‑148 (MQ=255)
            gCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGc                                                                                                                                        >  3:801262/1‑149 (MQ=255)
                tACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAg                                                                                                                                    >  3:809884/1‑149 (MQ=255)
                   tGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCa                                                                                                                                 <  3:266502/149‑1 (MQ=255)
                   tGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCa                                                                                                                                 >  3:381449/1‑149 (MQ=255)
                     ggTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCAc                                                                                                                                >  2:910999/1‑148 (MQ=255)
                         aTGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGGACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                           >  4:856393/1‑149 (MQ=255)
                            aCGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagta                                                                                                                       <  3:799768/149‑1 (MQ=255)
                                aGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAgtg                                                                                                                    >  4:46717/1‑149 (MQ=255)
                                     cATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGc                                                                                                               <  4:883144/149‑1 (MQ=255)
                                         cTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAg                                                                                                           <  3:420825/149‑1 (MQ=255)
                                          tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCt                                                                                                                                                     <  3:374246/106‑1 (MQ=255)
                                          tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCt                                                                                                                                                     >  4:374246/1‑106 (MQ=255)
                                          tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAgg                                                                                                          >  4:121674/1‑149 (MQ=255)
                                           ccAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagt                                                                                                                        <  1:772615/134‑1 (MQ=255)
                                           ccAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagt                                                                                                                        >  2:772615/1‑134 (MQ=255)
                                              aagaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCg                                                                                                      <  3:598927/146‑1 (MQ=255)
                                              aaaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTc                                                                                                       <  2:258085/148‑1 (MQ=255)
                                                    gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGaa                                                                                                <  3:15630/149‑1 (MQ=255)
                                                        tGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGAcc                                                                                            <  4:809884/149‑1 (MQ=255)
                                                        tGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCta                                                                                                                                            >  1:80372/1‑100 (MQ=255)
                                                        tGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCta                                                                                                                                            <  2:80372/101‑2 (MQ=255)
                                                             tACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTc                                                                                       >  3:373092/1‑149 (MQ=255)
                                                                   tttACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCa                                                                                  <  1:379866/148‑1 (MQ=255)
                                                                        tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                            <  4:175558/149‑1 (MQ=255)
                                                                        tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                            <  3:121674/149‑1 (MQ=255)
                                                                         aaGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCTATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGa                                                                            >  2:750722/1‑148 (MQ=255)
                                                                          agagAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                          >  4:715806/1‑149 (MQ=255)
                                                                           gagaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                         >  3:189098/1‑149 (MQ=255)
                                                                            agaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                         <  2:636776/148‑1 (MQ=255)
                                                                            agaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCg                                                                         >  1:836738/1‑148 (MQ=255)
                                                                            agaATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGt                                                                        >  4:114765/1‑149 (MQ=255)
                                                                                    aTGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                                <  3:46717/149‑1 (MQ=255)
                                                                                     tGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTaa                                                                <  1:567646/148‑1 (MQ=255)
                                                                                       aGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCa                                                                                      <  1:285010/124‑1 (MQ=255)
                                                                                       aGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCa                                                                                      >  2:285010/1‑124 (MQ=255)
                                                                                         cAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCt                                                            >  1:645736/1‑148 (MQ=255)
                                                                                             aaCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCa                                                                                                                          >  1:419020/1‑82 (MQ=255)
                                                                                             aaCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCa                                                                                                                          <  2:419020/82‑1 (MQ=255)
                                                                                                ttCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcca                                                    >  4:159816/1‑149 (MQ=255)
                                                                                                   aGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCt                                                 >  3:899990/1‑149 (MQ=255)
                                                                                                     aCATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCtt                                                <  2:836738/148‑1 (MQ=255)
                                                                                                       aTCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCtttt                                              >  1:252445/1‑148 (MQ=255)
                                                                                                        tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                            >  3:524448/1‑149 (MQ=255)
                                                                                                        tCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                            >  3:250050/1‑149 (MQ=255)
                                                                                                         cAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                            >  2:56076/1‑148 (MQ=255)
                                                                                                          aaCTATTGGGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                          >  4:78274/1‑149 (MQ=255)
                                                                                                           aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAccc                                         >  4:244803/1‑149 (MQ=255)
                                                                                                           aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAcc                                          >  2:353728/1‑148 (MQ=255)
                                                                                                            cTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAccc                                         >  2:817725/1‑148 (MQ=255)
                                                                                                             tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                       <  4:381449/149‑1 (MQ=255)
                                                                                                             tATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                       <  4:801262/149‑1 (MQ=255)
                                                                                                              aTTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCaa                                       <  2:583165/148‑1 (MQ=255)
                                                                                                                   gTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCa                                  >  1:503570/1‑148 (MQ=255)
                                                                                                                    tCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGc                                <  4:373092/149‑1 (MQ=255)
                                                                                                                    tCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAA‑‑CAACTCCAGcaa                              >  3:424208/1‑149 (MQ=255)
                                                                                                                      tGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGca                               <  2:645736/148‑1 (MQ=255)
                                                                                                                            tACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACt                         <  1:750722/148‑1 (MQ=255)
                                                                                                                              cTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACttt                       <  2:503570/148‑1 (MQ=255)
                                                                                                                              cTGGTTGGTATAACGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACttt                       <  2:243438/148‑1 (MQ=255)
                                                                                                                               tggttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCt                     >  3:272008/1‑149 (MQ=255)
                                                                                                                                    ggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaa                >  3:439504/1‑149 (MQ=255)
                                                                                                                                    ggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaa                <  3:159816/149‑1 (MQ=255)
                                                                                                                                      tataGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatc              <  3:469806/149‑1 (MQ=255)
                                                                                                                                       ataGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACttt                       >  4:317406/1‑139 (MQ=255)
                                                                                                                                       ataGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACttt                       <  3:317406/139‑1 (MQ=255)
                                                                                                                                        aaGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                            <  4:574846/116‑1 (MQ=255)
                                                                                                                                        aaGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTaa                                            >  3:574846/2‑117 (MQ=255)
                                                                                                                                        aaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTc                      >  1:389912/2‑139 (MQ=255)
                                                                                                                                        aaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTc                      <  2:389912/138‑1 (MQ=255)
                                                                                                                                         aGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                        <  2:365887/136‑1 (MQ=255)
                                                                                                                                         aGCGGCCCCCTGGGCCGCTTTTGCCCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                        >  1:365887/1‑136 (MQ=255)
                                                                                                                                             gCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGa       <  3:715806/149‑1 (MQ=255)
                                                                                                                                              cccccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGa       >  2:508969/1‑148 (MQ=255)
                                                                                                                                                cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat    <  4:899990/149‑1 (MQ=255)
                                                                                                                                                cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat    <  4:439504/149‑1 (MQ=255)
                                                                                                                                                cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGatat    >  3:298904/1‑149 (MQ=255)
                                                                                                                                                 ccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt   <  4:189098/149‑1 (MQ=255)
                                                                                                                                                 ccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt   <  3:244803/149‑1 (MQ=255)
                                                                                                                                                  cTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt   >  1:60312/1‑148 (MQ=255)
                                                                                                                                                  cTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTg  >  3:833930/1‑149 (MQ=255)
                                                                                                                                                   tGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTg  >  2:204160/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTG  >  NC_000913/779514‑779808

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