Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,236,299 T→A 5.7% L13* (TTA→TAA)  yjbF → lipoprotein YjbF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,236,2990TA 5.7% 281.6 / 4.6 88L13* (TTA→TAA) yjbFlipoprotein YjbF
Reads supporting (aligned to +/- strand):  ref base T (39/44);  new base A (2/3);  total (41/47)
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.95e-01

TATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTAT  >  NC_000913/4236154‑4236448
                                                                                                                                                 |                                                                                                                                                     
tatgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCTGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAg                                                                                                                                                   <  2:970489/150‑1 (MQ=255)
  tgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGc                                                                                                                                                 <  2:1308193/150‑1 (MQ=255)
   gtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGcc                                                                                                                                                <  2:787788/150‑1 (MQ=255)
   gtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGcc                                                                                                                                                <  2:1405404/150‑1 (MQ=255)
    tatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCt                                                                                                                                               >  2:1079234/1‑150 (MQ=255)
    tatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTAATAACAGGcag                                                                                                                                               <  2:13883/150‑3 (MQ=255)
     atCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTg                                                                                                                                              >  2:277364/1‑150 (MQ=255)
       ccccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGtt                                                                                                                                            >  2:893050/1‑150 (MQ=255)
        cccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTc                                                                                                                                           >  2:727694/1‑150 (MQ=255)
        cccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTc                                                                                                                                           >  2:374629/1‑150 (MQ=255)
        cccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTAATAACAGGcagataa                                                                                                                                           <  1:778205/150‑7 (MQ=255)
        cccAAAATAATTCGAGTCATTGCAGCTGTGGCTAGCAGTCTGAAGGGCGTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTAATAACAGGcagataa                                                                                                                                           <  1:842646/150‑7 (MQ=255)
              ataatTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCAc                                                                                                                                     >  1:1496668/1‑150 (MQ=255)
                 atTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAc                                                                                                                                  <  1:12412/150‑1 (MQ=255)
                 atTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAc                                                                                                                                  <  1:972383/150‑1 (MQ=255)
                 atTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAc                                                                                                                                  >  1:1329812/1‑150 (MQ=255)
                      aGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAag                                                                                                                             <  1:499751/150‑1 (MQ=255)
                         cATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCATATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAagagg                                                                                                                          >  1:1091615/1‑149 (MQ=255)
                           ttGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTa                                                                                                                                >  1:53651/1‑142 (MQ=255)
                           ttGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTa                                                                                                                                <  2:53651/142‑1 (MQ=255)
                           ttGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTg                                                                                                                        >  1:87845/1‑150 (MQ=255)
                            tGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGGGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTgg                                                                                                                       <  2:1029267/150‑1 (MQ=255)
                              cATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGc                                                                                                                     <  1:1539480/150‑1 (MQ=255)
                                tCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCaa                                                                                                                   <  1:335621/150‑1 (MQ=255)
                                tCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCaa                                                                                                                   <  2:806302/150‑1 (MQ=255)
                                  tgtgGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAAtt                                                                                                                 >  2:696676/1‑150 (MQ=255)
                                       cTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtg                                                                                                            >  1:159510/1‑150 (MQ=255)
                                        tAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGCCCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtgt                                                                                                           >  1:1020384/1‑150 (MQ=255)
                                        tAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtgt                                                                                                           <  1:696676/150‑1 (MQ=255)
                                          gAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTgg                                                                                                         <  1:746284/150‑1 (MQ=255)
                                          gAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTgg                                                                                                         <  1:1308819/150‑1 (MQ=255)
                                            aGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGGCCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGa                                                                                                       >  1:347246/1‑150 (MQ=255)
                                                 gAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTAACTGTGGGACAGTc                                                                                                  >  1:441490/1‑150 (MQ=255)
                                                 gAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCAGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTc                                                                                                  >  2:1498434/1‑150 (MQ=255)
                                                    gggATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTggg                                                                                                        >  2:1083201/1‑141 (MQ=255)
                                                    gggATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTggg                                                                                                        <  1:1083201/141‑1 (MQ=255)
                                                      gATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTc                                                                                             >  1:164598/1‑150 (MQ=255)
                                                        ttAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCgg                                                                                           <  2:1329812/150‑1 (MQ=255)
                                                          aaccataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCa                                                                                         <  1:407168/150‑1 (MQ=255)
                                                          aaccataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCa                                                                                         >  2:1367322/1‑150 (MQ=255)
                                                          aaccataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCa                                                                                         >  1:1352919/1‑150 (MQ=255)
                                                           accataaccaCACTCCGGTGTGGTTATTATGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTTCCTGTTATTACAGGCCTGTTAAGCCAAGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCAc                                                                                        >  2:595457/1‑150 (MQ=255)
                                                            ccataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGCCCTGCACTCATTCTTATCAGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACg                                                                                       >  1:31691/1‑150 (MQ=255)
                                                            ccataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACg                                                                                       <  1:501233/150‑1 (MQ=255)
                                                              ataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGcc                                                                                     <  1:804204/150‑1 (MQ=255)
                                                                  ccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGc                                                                                 >  2:1536418/1‑150 (MQ=255)
                                                                   cacaCTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCg                                                                                <  2:1496668/150‑1 (MQ=255)
                                                                   cacaCTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCg                                                                                <  1:1308475/150‑1 (MQ=255)
                                                                   cacaCTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATATGCCTGTTATTTCAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCg                                                                                >  1:210072/1‑150 (MQ=255)
                                                                      acTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTAc                                                                             >  1:659060/1‑150 (MQ=255)
                                                                                ggTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGga                                                                   <  2:165019/150‑1 (MQ=255)
                                                                                ggTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGga                                                                   <  1:1079234/150‑1 (MQ=255)
                                                                                 gTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGAc                                                                                                      >  2:1417115/1‑114 (MQ=255)
                                                                                 gTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGAc                                                                                                      <  1:1417115/114‑1 (MQ=255)
                                                                                 gTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAg                                                                                                    >  2:166776/1‑116 (MQ=255)
                                                                                 gTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAg                                                                                                    <  1:166776/116‑1 (MQ=255)
                                                                                      tCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGgatgatga                                                             >  1:1089170/1‑150 (MQ=255)
                                                                                      tCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGgatgatga                                                             <  1:307972/150‑1 (MQ=255)
                                                                                      tCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGgatgatga                                                             <  1:277364/150‑1 (MQ=255)
                                                                                      tCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACACCTGACGgatgatga                                                             >  1:230633/1‑150 (MQ=255)
                                                                                       cTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGgatgatgat                                                            <  2:659060/150‑1 (MQ=255)
                                                                                       cTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGgatgatgat                                                            <  1:1536418/150‑1 (MQ=255)
                                                                                        tGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGata                                                           >  2:494836/1‑150 (MQ=255)
                                                                                         gCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGatat                                                          >  2:1348326/1‑150 (MQ=255)
                                                                                           cccTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTc                                                        <  1:487338/150‑1 (MQ=255)
                                                                                                tGGAGATGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAAATCACTGTGTGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAat                                                   <  1:574325/150‑1 (MQ=255)
                                                                                                 ggAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATATGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGATCTGGGCAATTCAATGTCGGACAGTCTGTTCGGCCCGCCAGCCGTACAGCTTACGGATGATGATATTCAAAata                                                  >  2:655465/1‑150 (MQ=255)
                                                                                                   agaagaGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCaaa                                                     >  1:556427/1‑145 (MQ=255)
                                                                                                   agaagaGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCaaa                                                     <  2:556427/145‑1 (MQ=255)
                                                                                                    gaagaGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGc                                               >  2:1627243/1‑150 (MQ=255)
                                                                                                     aagaGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGcc                                              <  2:159510/150‑1 (MQ=255)
                                                                                                       gagTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCt                                            >  2:110016/1‑150 (MQ=255)
                                                                                                       gagGCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGGTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCt                                            <  2:379126/150‑1 (MQ=255)
                                                                                                         gTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTAc                                          <  2:140316/150‑1 (MQ=255)
                                                                                                           cGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgc                                        >  2:612338/1‑150 (MQ=255)
                                                                                                            gTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgcc                                       <  1:374629/150‑1 (MQ=255)
                                                                                                              gAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgccag                                     <  1:727694/150‑1 (MQ=255)
                                                                                                                    aCCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTa                               <  1:494836/150‑1 (MQ=255)
                                                                                                                     ccTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTAc                              <  2:347246/150‑1 (MQ=255)
                                                                                                                              aTTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCtt                     <  2:1622467/150‑1 (MQ=255)
                                                                                                                                tCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTaa                   <  2:31691/150‑1 (MQ=255)
                                                                                                                                  gaacaggCCTGTAATAACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATg                 >  2:1452629/7‑150 (MQ=255)
                                                                                                                                   aacaggCCTGTAATAACAGGCCTGTTCAGCCACGACTAATGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATgg                >  2:185520/6‑150 (MQ=255)
                                                                                                                                    aTCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGc               >  1:1192508/1‑150 (MQ=255)
                                                                                                                                         ccTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGcc          <  2:193651/150‑1 (MQ=255)
                                                                                                                                          cTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgccag                                     <  1:1204440/122‑1 (MQ=255)
                                                                                                                                          cTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgccag                                     >  2:1204440/1‑122 (MQ=255)
                                                                                                                                          cTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCg         >  1:1079152/1‑150 (MQ=255)
                                                                                                                                                 taCAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTAt  >  1:618752/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTAT  >  NC_000913/4236154‑4236448

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