Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 3,539,413 C→T 8.3% L677L (CTG→TTG)  yhgF → putative transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,539,4130CT 8.3% 161.2 / 3.5 62L677L (CTG→TTG) yhgFputative transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGATCCGCGTCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTG  >  NC_000913_3_pae_pgi/3539265‑3539560
                                                                                                                                                    |                                                                                                                                                   
cgATCCGCGTCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGcc                                                                                                                                                     >  1:250795/1‑149 (MQ=255)
cgATCCGCGTCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGcc                                                                                                                                                     <  1:877314/149‑1 (MQ=255)
 gATCCGCGTCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCCAGTTTGTCGATATTGGCGTGCATCAGGACGGccg                                                                                                                                                    >  2:893802/1‑148 (MQ=255)
         ttcGGAATTTAAACCCGCTCAGGTTGCCGATGGCGTCGAGACAATGATCGGCCTGTAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTcaca                                                                                                                                            <  2:159517/147‑1 (MQ=255)
         tCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTcaca                                                                                                                                            <  2:159594/149‑1 (MQ=255)
         gccGGGATTTAAAACCGCTCAGGGTGCCGATGGCGTCGAGACCATGAAAGACCTCCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTcaca                                                                                                                                            <  2:801855/148‑1 (MQ=255)
            ggAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGAGCGATTGTCAATATTCGCGTGCATCAGGAAGGACTGGTTCACAtct                                                                                                                                         >  1:938561/1‑149 (MQ=255)
             gAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACAtctc                                                                                                                                        <  1:356311/149‑1 (MQ=255)
              aaTTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGAGCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCCATACTGACTTGGATCAGGACGGCCTGGTTCACATctct                                                                                                                                       <  2:720294/149‑1 (MQ=255)
                    aaaCCGCTCAGTTTGCCGATGGCGACGAGACAATGAACGACCAGCAACCGGGTATGACCCTCGAAGGCGCAGTGACCAACGTCCCCAACTTTGGCGCGTATGTCGATATGGGCGTGCAGAAGGACCGGCTGGTTCACATCTCTTCAtgg                                                                                                                                 >  1:109519/1‑147 (MQ=255)
                    aaaCCGACCAGGTCGCTGATCGCGTCGAGACACTGAACTACCTTCAACCGGGTATGATCCTCGAAGGCGCAGTGCCCAACGACAACAACTTCGGCGCCTTTGTCGATATTGTCGTGCATCAGGACTTCCGGGTTCACATCTCTTCATTg                                                                                                                                 <  2:164733/149‑1 (MQ=255)
                     aaCCGCTCAGTTTGCCGATGTCGTCGAGACGATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGATCAGATCTCTTCATTGt                                                                                                                                >  1:976741/1‑149 (MQ=255)
                      aCCGCTCAGTTTGCCGATTGCGTCAAGACAATGACCGACCTGCAACCTTGTATGAGCCTCGAACGCGCAGTGACCAATGGCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTc                                                                                                                               <  1:893802/149‑1 (MQ=255)
                          cTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTCGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAAc                                                                                                                           >  1:894628/1‑149 (MQ=255)
                                    cGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTATGGCGCGTTTGTCGATATTGGAGAGCATCAGGACGGCCTGGTTCACATCACTTCATTGTCGAACAAGTTTGTgg                                                                                                                 >  1:354837/1‑149 (MQ=255)
                                     gATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGCCCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCACGACGGCCTGGTTCACATCTCTTCATTGTCTACCCAGTGTGTGGa                                                                                                                >  2:78265/1‑149 (MQ=255)
                                      aTGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGaa                                                                                                               <  1:913504/149‑1 (MQ=255)
                                           gTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATcc                                                                                                          >  2:514021/1‑149 (MQ=255)
                                           gTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACAAACTTTGGCGCGGTCGTCGAGATTGGCGTGCATCAGGACGCCCTGGTTTACATCACTTCATCATCGAACAGATTTGTGGAACATcc                                                                                                          >  2:918300/1‑149 (MQ=255)
                                                   aaTGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACAGCCTGGTTCACATCTATTCATTGTCGAACAAGTTTGTGAAACATCGCCATCCCg                                                                                                  >  2:15144/1‑149 (MQ=255)
                                                      gatcGACCTGCAAACGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATGGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCgtgg                                                                                               <  1:269969/146‑1 (MQ=255)
                                                              tGCAACCGGTTATGAGCCTCGAAGGCGCAGTGACCAACGTCACAAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCg                                                                                       <  2:793672/149‑1 (MQ=255)
                                                                cAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCAtt                                                                                                                                  >  1:839979/1‑104 (MQ=255)
                                                                cAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCAtt                                                                                                                                  <  2:839979/104‑1 (MQ=255)
                                                                       gTATGATCCGCGAAGGTGGCGGGACCATAGTCACAAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACAtt                                                                              <  2:591623/149‑1 (MQ=255)
                                                                           gctcgtCGAAGTGGCAGGGGCCAACGTGACCAACTTTGGCGCGTTTGTCGATAGTAGTGTGCATCAGGACGGACAGATTCACAGCTCTTCATTGTCGAACCAGTTTGTGGAAGATCCGCATACCGTGTTGAAATCGGGCGACATTGTGa                                                                          <  2:985580/144‑1 (MQ=255)
                                                                           gATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGa                                                                          >  1:507222/1‑149 (MQ=255)
                                                                            aTCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGaa                                                                         >  1:6006/1‑149 (MQ=255)
                                                                            aTCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCCTGTTGAAAGCGGGCGACACTGTGaa                                                                         >  1:729800/1‑149 (MQ=255)
                                                                             tCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCGTTGTCGAACAAGGTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGGGACGTTGTGAAg                                                                        >  1:909713/1‑149 (MQ=255)
                                                                             tCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCg                                                                                       >  2:833753/1‑134 (MQ=255)
                                                                             tCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCg                                                                                       <  1:833753/134‑1 (MQ=255)
                                                                              gcgcGAAGGCGGAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCTTGGCTCACATCTCTTCTTTTGCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAgg                                                                       <  2:416853/146‑1 (MQ=255)
                                                                               cTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGGCGGCCTGGTTCACATCTCTTCATTGGCGAACAAGTTTGTGGGAGCTTGACATCCCTGGGTGAAAGCGGGCGACGTtgcagaggt                                                                      >  1:779813/1‑142 (MQ=255)
                                                                                       cgcAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGc                                                              <  1:1004071/149‑1 (MQ=255)
                                                                                         cAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATCTTGGCTTGCATCAGGACGGCTTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTg                                                            >  1:144944/1‑149 (MQ=255)
                                                                                         cAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCTTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAAGTGAAAGTGCTg                                                            >  1:428987/1‑149 (MQ=255)
                                                                                            tGACCAACGTCACCAACGTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACACTTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGaa                                                         <  1:484659/149‑1 (MQ=255)
                                                                                                 aaCGTCACCAACTGTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAACCGCTGGAAGTgca                                                    <  2:492648/149‑3 (MQ=255)
                                                                                                  aCGGCACCAACTTTGGCGCGGTGGGCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGTAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGAt                                                   <  2:296125/149‑1 (MQ=255)
                                                                                                   cGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGTCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATc                                                  <  1:753092/149‑1 (MQ=255)
                                                                                                     tCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCtt                                                >  1:370325/1‑149 (MQ=255)
                                                                                                         cAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGc                                            <  1:184630/149‑1 (MQ=255)
                                                                                                             tttGGCGCGTTTGTCGATATTGGCATGCATCAGGACGGCCTGGTTCACAGCTCTTCATTGTCGAACACGGTTGTGGAAGATCCGCATACGGTGGTGAAAGCGGCCGACATTGTGAAGGTGAAAGTGCTGGGAGGGGACCCTCAGCGTaa                                        >  2:579091/1‑149 (MQ=255)
                                                                                                               ttgcgcGTTTGTCGATATTGGCGTGCATCAGGACGGCTTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAAc                                      <  2:144944/147‑1 (MQ=255)
                                                                                                               tGGCGGGTTTGTCGATATTGGCGTGCATCAGGACGGCTTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAAc                                      <  2:428987/149‑1 (MQ=255)
                                                                                                                  cgcgTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTCAAAACGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGCTCTTCGGCGTGAACATa                                   >  2:264946/1‑149 (MQ=255)
                                                                                                                   gggTTTGTCGATTTTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTAt                                  <  2:860956/147‑1 (MQ=255)
                                                                                                                     gTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCg                                >  1:50872/1‑149 (MQ=255)
                                                                                                                         gTCGATATTGGCGTGCATCAGGACGGCATGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGGCCt                            >  2:732781/1‑149 (MQ=255)
                                                                                                                         gTCAATATTGGCTTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCt                            <  1:465283/149‑1 (MQ=255)
                                                                                                                               atTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGTTCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTAt                      <  1:81375/149‑1 (MQ=255)
                                                                                                                                ttGGTGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATg                     <  2:176204/149‑1 (MQ=255)
                                                                                                                                   gCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATgcgc                  >  2:924272/1‑149 (MQ=255)
                                                                                                                                   gCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTATGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGGCCTGACTATgcgc                  >  2:850329/1‑149 (MQ=255)
                                                                                                                                     gTGCATCAGGACGGTCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCt                <  1:58323/149‑1 (MQ=255)
                                                                                                                                     gTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCCCCt                >  2:934492/1‑149 (MQ=255)
                                                                                                                                      tGCATCAGGAGGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTg               <  2:709325/149‑1 (MQ=255)
                                                                                                                                         agcAGGACTGCTTGGTTAACATCTTTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGATATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGAt            <  2:768028/147‑1 (MQ=255)
                                                                                                                                          tCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGGGGAGAAAGACGGCGACATTGGGCAGGAGAGAGTGCAGTAAGTGGGTCTTCAGCGTAACAGTATCGCCGTGGCTAGTGTCCTGGATg           >  2:57436/1‑149 (MQ=255)
                                                                                                                                             ggACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGagc        >  2:851549/1‑149 (MQ=255)
                                                                                                                                              gACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGagca       >  1:536852/1‑149 (MQ=255)
                                                                                                                                              gACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCATGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATTagca       >  1:536630/1‑149 (MQ=255)
                                                                                                                                               aCGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGagcag      >  2:229762/1‑149 (MQ=255)
                                                                                                                                                cGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGTGTAAACGTATCACCCTGACTATGAGCCTGGATGagcagc     >  1:126770/1‑149 (MQ=255)
                                                                                                                                                   ccTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTTAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTg  <  1:256302/149‑1 (MQ=255)
                                                                                                                                                   ccTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAACGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGTTCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGAAGCCTg  >  1:167145/1‑149 (MQ=255)
                                                                                                                                                   ccTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACAtt                                                                              <  2:927081/73‑1 (MQ=255)
                                                                                                                                                   ccTGGTTCACATCTCTTCATTGTAGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACAtt                                                                              >  1:927081/1‑73 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATCCGCGTCCGGAATTTAAAACCGCTCAGTTTGCCGATGGCGTCGAGACAATGAACGACCTGCAACCGGGTATGATCCTCGAAGGCGCAGTGACCAACGTCACCAACTTTGGCGCGTTTGTCGATATTGGCGTGCATCAGGACGGCCTGGTTCACATCTCTTCATTGTCGAACAAGTTTGTGGAAGATCCGCATACCGTGGTGAAAGCGGGCGACATTGTGAAGGTGAAAGTGCTGGAAGTGGATCTTCAGCGTAAACGTATCGCCCTGACTATGCGCCTGGATGAGCAGCCTG  >  NC_000913_3_pae_pgi/3539265‑3539560

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