Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,399,691 C→A N48K (AAC→AAA hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,399,6910CA96.4% 282.3 / 0.6 83N48K (AAC→AAAhfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (39/41);  total (40/43)
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.68e-01

AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGG  >  NC_000913_3_pae_pgi/4399545‑4399836
                                                                                                                                                  |                                                                                                                                                 
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                <  3:336534/150‑1 (MQ=255)
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                <  1:466335/150‑1 (MQ=255)
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                >  4:486453/1‑150 (MQ=255)
  aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGt                                                                                                                                              >  2:797836/1‑150 (MQ=255)
  aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTAGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGt                                                                                                                                              >  2:845279/1‑150 (MQ=255)
               cAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGtt                                                                                                                                 <  3:202897/150‑1 (MQ=255)
                  tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTAc                                                                                                                              <  1:50015/150‑1 (MQ=255)
                     ttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAg                                                                                                                           <  1:1014995/150‑1 (MQ=255)
                     ttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAg                                                                                                                           >  1:126025/1‑150 (MQ=255)
                       aCAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCa                                                                                                                         <  2:835624/150‑1 (MQ=255)
                        cAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                        <  1:212989/150‑1 (MQ=255)
                            aTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGa                                                                                                                    <  3:486453/150‑1 (MQ=255)
                                gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc                                                                                                                >  1:1111633/1‑150 (MQ=255)
                                    cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                                                                            >  3:578109/1‑150 (MQ=255)
                                    cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                                                                            >  2:1105509/1‑150 (MQ=255)
                                    cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                                                                            >  3:412561/1‑150 (MQ=255)
                                     tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  3:473272/1‑150 (MQ=255)
                                     tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  2:690381/1‑150 (MQ=255)
                                       aaCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgtt                                                                                                         >  2:947925/1‑150 (MQ=255)
                                           cACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc                                                                                                     >  3:429195/1‑150 (MQ=255)
                                               gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                                                                 <  3:216489/150‑1 (MQ=255)
                                                    gggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                                                                                                               >  4:222668/1‑131 (MQ=255)
                                                    gggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                                                                                                               <  3:222668/131‑1 (MQ=255)
                                                      gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                          <  2:126025/150‑1 (MQ=255)
                                                      gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                          <  2:1111633/150‑1 (MQ=255)
                                                       aaCGTGTACCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGCTCAGTTCGTGATCCTGTAGAAAAAAACGTTCCGCCAGCTGGCTTACCCGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCgg                                                                                         >  2:179689/1‑150 (MQ=255)
                                                          tttttCCAGTTTCTATTTATTTGGTTAATGTTATTAAGCTGCAAGCGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTAAAAAAAACGGTCATCCAGATGGTTTACAAGCACGCCATTTCTACTGTTGTTCCGTCTCGCCCGGttt                                                                                      <  2:580677/147‑1 (MQ=255)
                                                           tgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                                                     <  1:1105509/150‑1 (MQ=255)
                                                            gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                    <  4:458744/150‑1 (MQ=255)
                                                                 aGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                               <  1:690381/150‑1 (MQ=255)
                                                                 aGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                               >  2:853806/1‑150 (MQ=255)
                                                                 aGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                               <  4:429195/150‑1 (MQ=255)
                                                                   tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                             >  3:16945/1‑150 (MQ=255)
                                                                     tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGt                                                                           <  1:533279/150‑1 (MQ=255)
                                                                      cTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  2:196950/1‑117 (MQ=255)
                                                                      cTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           <  1:196950/117‑1 (MQ=255)
                                                                         tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaaca                                                                       >  1:63024/1‑150 (MQ=255)
                                                                            atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  2:262815/1‑111 (MQ=255)
                                                                            atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           <  1:262815/111‑1 (MQ=255)
                                                                            atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACg                                                                    >  4:660492/1‑150 (MQ=255)
                                                                             tttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCa                                                                                                                                        <  1:1039656/81‑1 (MQ=255)
                                                                             tttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCa                                                                                                                                        >  2:1039656/1‑81 (MQ=255)
                                                                                ggTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc                                                                                                     <  4:196176/113‑1 (MQ=255)
                                                                                ggTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc                                                                                                     >  3:196176/1‑113 (MQ=255)
                                                                                 gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                               <  3:531134/150‑1 (MQ=255)
                                                                                 gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                               >  1:869319/1‑150 (MQ=255)
                                                                                   gAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                                             <  1:414119/150‑1 (MQ=255)
                                                                                     aTGGTATTAGGCTGCAGTGGCAAATCGAGTCTTTTGATGAGTTCGTGATAATGTTAAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                           <  4:120992/150‑1 (MQ=255)
                                                                                      tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGa                                                                                                                                      <  1:265366/74‑1 (MQ=255)
                                                                                      tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGa                                                                                                                                      >  2:265366/1‑74 (MQ=255)
                                                                                           ttAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACca                                                     >  3:632582/1‑150 (MQ=255)
                                                                                            tAAGCTGCAAGGGCAAATCGAGTCTGTTGATCAGTGCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGATACcag                                                    <  2:350422/150‑1 (MQ=255)
                                                                                                cTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGt                                                >  3:56476/1‑150 (MQ=255)
                                                                                                 tGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTa                                               >  1:1097545/1‑150 (MQ=255)
                                                                                                  gCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCGGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa                                              >  3:647192/1‑150 (MQ=255)
                                                                                                  gCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa                                              <  2:63024/150‑1 (MQ=255)
                                                                                                   cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                                                                                                      >  4:237663/1‑93 (MQ=255)
                                                                                                   cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                                                                                                      <  3:237663/93‑1 (MQ=255)
                                                                                                          aaaTCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatc                                      >  1:804811/1‑150 (MQ=255)
                                                                                                              cGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCGGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATgg                                  <  4:647192/150‑1 (MQ=255)
                                                                                                                 gTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                               <  4:56476/150‑1 (MQ=255)
                                                                                                                    ttttGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagcag                            >  2:1035495/1‑150 (MQ=255)
                                                                                                                       tGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgc                         <  1:947925/150‑1 (MQ=255)
                                                                                                                         aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                       >  2:315852/1‑150 (MQ=255)
                                                                                                                         aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                       <  1:845279/150‑1 (MQ=255)
                                                                                                                             gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                    >  2:613289/1‑85 (MQ=255)
                                                                                                                             gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                    <  1:613289/85‑1 (MQ=255)
                                                                                                                              ttCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAAt                  >  3:80426/1‑150 (MQ=255)
                                                                                                                                   gATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtctc                                                                                   <  2:553506/80‑1 (MQ=255)
                                                                                                                                   gATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtctc                                                                                   >  1:553506/1‑80 (MQ=255)
                                                                                                                                    aTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                      <  4:113715/76‑1 (MQ=255)
                                                                                                                                    aTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                      >  3:113715/1‑76 (MQ=255)
                                                                                                                                    aTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc            <  1:854323/150‑1 (MQ=255)
                                                                                                                                    aTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc            <  4:473272/150‑1 (MQ=255)
                                                                                                                                      ccTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc          <  2:804811/150‑1 (MQ=255)
                                                                                                                                      ccTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc          <  1:334633/150‑1 (MQ=255)
                                                                                                                                      ccTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACAGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc          <  1:334908/150‑1 (MQ=255)
                                                                                                                                       cTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgcg         >  3:83836/1‑150 (MQ=255)
                                                                                                                                          ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                       <  4:234577/133‑1 (MQ=255)
                                                                                                                                          ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                       >  3:234577/1‑133 (MQ=255)
                                                                                                                                          ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa      >  1:908131/1‑150 (MQ=255)
                                                                                                                                          ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa      <  4:632582/150‑1 (MQ=255)
                                                                                                                                            gAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACa    >  4:658875/1‑150 (MQ=255)
                                                                                                                                             aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAg   >  1:91596/1‑150 (MQ=255)
                                                                                                                                              aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatcat                                    <  1:382519/116‑1 (MQ=255)
                                                                                                                                              aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatcat                                    >  2:382519/1‑116 (MQ=255)
                                                                                                                                              aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAgg  >  2:645007/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGG  >  NC_000913_3_pae_pgi/4399545‑4399836

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