Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,454,209 T→C 5.3% E397G (GAA→GGA)  gspA ← general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,454,2090TC 5.3% 235.2 / 2.9 75E397G (GAA→GGA) gspAgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base C (2/2);  total (35/40)
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.08e-01

TTGAAATTGCTTTATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAATTTTTTGTC  >  NC_000913/3454067‑3454354
                                                                                                                                              |                                                                                                                                                 
ttGAAATTGCTTTATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCgg                                                                                                                                               <  2:53748/147‑1 (MQ=255)
   aaaTTGCTTTATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGa                                                                                                                                            <  8:237953/147‑1 (MQ=255)
      ttGCTTTATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGATTCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACAtt                                                                                                                                        >  4:112106/1‑148 (MQ=255)
           ttATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCca                                                                                                                                   >  1:189351/1‑148 (MQ=255)
            aaTTTTATCTACCAACAAGCGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTAAACATCGATTCCAGCCAGTGAATCTGTTCCTCGCTTCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCca                                                                                                                                   <  6:108837/146‑1 (MQ=255)
               tttttCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAAGATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCacaa                                                                                                                                >  4:236047/1‑147 (MQ=255)
               tttttCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCacaa                                                                                                                                >  1:210453/1‑147 (MQ=255)
                   tCTACCAACAAGCGACGCCACTCCGCCGAAGGTTCTGGCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACaga                                                                                                                            <  6:285775/147‑1 (MQ=255)
                        caacaaGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATaa                                                                                                                      <  3:375551/148‑1 (MQ=255)
                         aacaaGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAAt                                                                                                                     >  7:324349/1‑148 (MQ=255)
                         aacaaGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAAt                                                                                                                     >  7:238818/1‑148 (MQ=255)
                          acaaGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAAt                                                                                                                     >  6:6681/1‑147 (MQ=255)
                             aGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCAcc                                                                                                                 <  4:255597/148‑1 (MQ=255)
                                gACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGt                                                                                                              <  3:253341/148‑1 (MQ=255)
                                    ccACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCa                                                                                                           >  5:197596/1‑147 (MQ=255)
                                     cACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATa                                                                                                         >  6:11658/1‑148 (MQ=255)
                                        tCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATAcc                                                                                                       >  3:311838/1‑147 (MQ=255)
                                         ccgccgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGa                                                                                                     >  6:187030/1‑148 (MQ=255)
                                         ccgccgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCACTTCAGGCGACATTTTCAACAAAATAAGATAATCACCCGTCCATACCGa                                                                                                     >  2:65457/1‑148 (MQ=255)
                                          cgccgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCg                                                                                                      <  5:176245/146‑1 (MQ=255)
                                          cgccgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGAc                                                                                                    <  6:197596/148‑1 (MQ=255)
                                             cgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTTCTATCACGCGTAATGGTACTCTCCCCTTCCTGCGACATTTTCCACAACAGAAGATAATCACCCGTCAATACCGACTc                                                                                                  >  5:82383/1‑147 (MQ=255)
                                             cgAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTc                                                                                                  >  1:305078/1‑147 (MQ=255)
                                                    tCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCAtt                                                                                          >  1:82925/1‑148 (MQ=255)
                                                       gTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTAt                                                                                       >  5:142526/1‑148 (MQ=255)
                                                        tCGAAATATGCAATGCGCGGTTTAACATCGTTTCCATCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATTGTACTCTCCCCTTCCTGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCAttt                                                                                         <  3:195053/145‑1 (MQ=255)
                                                        tCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCAttt                                                                                         <  1:293205/145‑1 (MQ=255)
                                                        tCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTAt                                                                                       <  1:289333/147‑1 (MQ=255)
                                                          gAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTAtgtg                                                                                    <  3:351352/148‑1 (MQ=255)
                                                            aaTATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAATATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTc                                                                                  <  2:83531/148‑1 (MQ=255)
                                                            aaTATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTc                                                                                  <  5:6681/148‑1 (MQ=255)
                                                                    aTGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCa                                                                          <  2:210453/148‑1 (MQ=255)
                                                                    aTGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCa                                                                          >  7:224263/1‑148 (MQ=255)
                                                                     tGCGCGGTTTTACATCGTTTCCAGCCAGATAATCTCTTCCTCGCTGCTATCACGAGTAATGGTACTCTCACCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCAc                                                                         <  8:357120/148‑1 (MQ=255)
                                                                      gcgcGGTTTAATATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTcc                                                                                                            >  1:54697/1‑112 (MQ=255)
                                                                      gcgcGGTTTAATATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTcc                                                                                                            <  2:54697/112‑1 (MQ=255)
                                                                         cgGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGttt                                                                     >  3:45219/1‑148 (MQ=255)
                                                                         cgGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGttt                                                                     <  8:224263/148‑1 (MQ=255)
                                                                          ggTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCACCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTg                                                                    <  5:239521/148‑1 (MQ=255)
                                                                           gTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGAGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCCGCGTCCACGTTTGc                                                                   >  3:64600/1‑148 (MQ=255)
                                                                               aaCATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGa                                                                <  5:11658/147‑1 (MQ=255)
                                                                                    cGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTAc                                                           >  5:126452/1‑147 (MQ=255)
                                                                                    cGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTCCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTAc                                                           >  5:244192/1‑147 (MQ=255)
                                                                                    cGTNTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTCCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTAc                                                           >  5:249855/1‑147 (MQ=255)
                                                                                     gtatCCATCAAGATAATCTCTTACTCGATGATATTAAGAGTAATGGTAATCTCCACTTACTGCGACATTTTCCACAACATAATATAATAACCAGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACg                                                          <  1:317840/144‑1 (MQ=255)
                                                                                     gTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACg                                                          <  3:236047/147‑1 (MQ=255)
                                                                                          cagccagAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACa                                                    <  8:281302/148‑1 (MQ=255)
                                                                                             ccagAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCa                                                 <  4:208157/148‑1 (MQ=255)
                                                                                             ccagAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTCCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCa                                                 <  6:244192/148‑1 (MQ=255)
                                                                                             ccagAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTCCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCa                                                 <  6:249855/148‑1 (MQ=255)
                                                                                                gagaATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACg                                              <  6:325370/148‑1 (MQ=255)
                                                                                                  taaTCTCTTCCTCGATGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCt                                            <  1:65457/147‑1 (MQ=255)
                                                                                                         ttCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGc                                         <  5:187030/144‑1 (MQ=255)
                                                                                                          tCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCAc                                    >  8:239472/1‑148 (MQ=255)
                                                                                                            cTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTACATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCg                                  >  8:140569/1‑148 (MQ=255)
                                                                                                                   ctATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGAcgcg                             <  8:238818/146‑1 (MQ=255)
                                                                                                                   ctATCACGCGTAATGGAACTCTCCCCTTCCGGCGACATTTTCCACACCAGAAGATAATCACCCGTACATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGa                            <  2:89804/147‑1 (MQ=255)
                                                                                                                       cACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTa                                                            <  1:292600/111‑1 (MQ=255)
                                                                                                                       cACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTa                                                            >  2:292600/1‑111 (MQ=255)
                                                                                                                              aaTGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGa                            >  1:233643/1‑136 (MQ=255)
                                                                                                                              aaTGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGa                            <  2:233643/136‑1 (MQ=255)
                                                                                                                              aaTGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCacaa                  >  1:13887/1‑146 (MQ=255)
                                                                                                                               aTGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCAATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACTAGGACATAAACGCTAGCATCACCGACGCGAACGACCacaac                 >  4:6406/1‑146 (MQ=255)
                                                                                                                                 ggTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAgg              >  8:167064/1‑147 (MQ=255)
                                                                                                                                   tACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCaa           <  7:167064/148‑1 (MQ=255)
                                                                                                                                    aCTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAt          >  1:335991/1‑148 (MQ=255)
                                                                                                                                     ctctCCCCTTCCGGCGACATTTTCCACAACATAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAt          <  2:49342/147‑1 (MQ=255)
                                                                                                                                       ctcCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAtttt       >  3:167275/1‑148 (MQ=255)
                                                                                                                                       ctcCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAtttt       >  1:311648/1‑148 (MQ=255)
                                                                                                                                        tcCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAttt        <  2:305078/146‑1 (MQ=255)
                                                                                                                                         ccccTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTAtg                                                                                      >  5:305988/1‑67 (MQ=255)
                                                                                                                                         ccccTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTAtg                                                                                      <  6:305988/67‑1 (MQ=255)
                                                                                                                                          cccTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAttttt      <  3:315180/146‑1 (MQ=255)
                                                                                                                                          cccTTCCGGCGACATTTTCCACAACAGAAGATAATCACACGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAAtttt       <  2:82925/145‑1 (MQ=255)
                                                                                                                                            cTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAATTTTTTGTc  >  8:1973/1‑148 (MQ=255)
                                                                                                                                              tCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCaa                                                <  3:184225/100‑1 (MQ=255)
                                                                                                                                              tCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCaa                                                >  4:184225/1‑100 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TTGAAATTGCTTTATTTTTTCTACCAACAAGGGACGCCACTCCGCCGAAGGTTCTGTCGAAATATGCAATGCGCGGTTTAACATCGTTTCCAGCCAGAGAATCTCTTCCTCGCTGCTATCACGCGTAATGGTACTCTCCCCTTCCGGCGACATTTTCCACAACAGAAGATAATCACCCGTCCATACCGACTCAAACCATTTATGTGTCAGCGTCCACGTTTGCTGACCTACGAGGACATCAACGCTAGCCTCACCGACGCGAACGACCACAACAGGCAATTTTTTGTC  >  NC_000913/3454067‑3454354

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