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,6910CA98.9% 290.3 / ‑4.5 88N48K (AAC→AAAhfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (33/53);  minor base G (1/0);  total (35/53)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGA  >  NC_000913_3_pae_pgi/4399544‑4399837
                                                                                                                                                    |                                                                                                                                                  
gagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGCTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaaac                                                                                                                                                  <  1:667490/150‑1 (MQ=255)
 agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                 <  2:424274/150‑1 (MQ=255)
 agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                 >  1:881853/1‑150 (MQ=255)
       gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGc                                                                                                                                           >  4:358413/1‑150 (MQ=255)
        cTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCACGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGcc                                                                                                                                          <  3:591812/150‑1 (MQ=255)
               gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGt                                                                                                                                   >  2:355427/1‑150 (MQ=255)
                  aTCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTa                                                                                                                                >  4:322264/1‑150 (MQ=255)
                       tACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGc                                                                                                                           <  3:533083/150‑1 (MQ=255)
                         cAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                         >  1:139117/1‑150 (MQ=255)
                         cAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                         <  1:269347/150‑1 (MQ=255)
                          aaGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAcg                                                                                                                        <  3:625447/150‑1 (MQ=255)
                           aGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAcgc                                                                                                                       <  1:886506/150‑1 (MQ=255)
                           aGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAcgc                                                                                                                       >  1:56557/1‑150 (MQ=255)
                             aTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGa                                                                                                                     <  3:107456/150‑1 (MQ=255)
                                 gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc                                                                                                                 >  4:9739/1‑150 (MQ=255)
                                  ttCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt                                                                                                                >  1:735420/1‑150 (MQ=255)
                                   tCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                                               <  2:684618/150‑1 (MQ=255)
                                     cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaaac                                                                                                                                                  <  3:628110/113‑1 (MQ=255)
                                     atGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCTGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                                                                             <  1:34280/149‑1 (MQ=255)
                                          cGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttgt                                                                                                        <  1:465217/150‑1 (MQ=255)
                                               tGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGt                                                                                                   <  4:416090/150‑1 (MQ=255)
                                               tGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGt                                                                                                   >  2:32578/1‑150 (MQ=255)
                                                gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                                                                  >  4:33512/1‑150 (MQ=255)
                                                gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCGTTTGATCAGTTCGTGAACCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                                                                  <  4:29986/150‑1 (MQ=255)
                                                    cGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGc                                                                                              <  1:355427/150‑1 (MQ=255)
                                                          cGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGtt                                                                                        >  2:659579/1‑150 (MQ=255)
                                                          cGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCACGGtt                                                                                        <  4:588335/150‑1 (MQ=255)
                                                           gtgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                       <  1:504017/150‑1 (MQ=255)
                                                            tgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCATTAGTTCAGTTCGTATTCTTTTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                                                      <  3:248059/150‑1 (MQ=255)
                                                             gtTTCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                     <  4:716794/150‑1 (MQ=255)
                                                              ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtctc                                                                                    >  2:889503/1‑150 (MQ=255)
                                                                ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                     >  4:450648/1‑147 (MQ=255)
                                                                ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                     <  3:450648/147‑1 (MQ=255)
                                                                ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                                                                  <  1:659579/150‑1 (MQ=255)
                                                                 cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatc                                                                                 >  3:144599/1‑150 (MQ=255)
                                                                    tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                              <  3:358413/150‑1 (MQ=255)
                                                                    tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                              <  2:56557/150‑1 (MQ=255)
                                                                    tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                              <  1:889503/150‑1 (MQ=255)
                                                                     ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAg                                                                             <  3:24820/150‑1 (MQ=255)
                                                                        tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                          >  4:237490/1‑150 (MQ=255)
                                                                        tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                          <  1:691643/150‑1 (MQ=255)
                                                                        tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGGTTTTGTATCAGTTCGTTAGCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                          <  1:167626/150‑1 (MQ=255)
                                                                            tatttGTTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGCTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaac                                                                      <  4:186572/150‑1 (MQ=255)
                                                                             atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACg                                                                     <  3:322264/150‑1 (MQ=255)
                                                                             atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACg                                                                     >  4:369138/1‑150 (MQ=255)
                                                                             atttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGGTTCTCATCACAGTAACAACg                                                                     >  4:30672/1‑150 (MQ=255)
                                                                                  gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGCTGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                                                         <  3:469101/125‑1 (MQ=255)
                                                                                  gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                                                         >  4:469101/1‑125 (MQ=255)
                                                                                  gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                                <  2:641793/150‑1 (MQ=255)
                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                                             <  1:101049/150‑1 (MQ=255)
                                                                                      aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGCAAAAAACGGTACGCCAGGTCGTTTACACGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                            <  3:359313/150‑1 (MQ=255)
                                                                                      aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                            <  3:669536/150‑1 (MQ=255)
                                                                                      aTGGGATTAAGCGGCAAGGGCAAATCGAGGCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                            <  4:511769/150‑1 (MQ=255)
                                                                                       tGGTATTAAGCTGCAAGGGCAAATCGGGTCTTTTGCTCAGTTCGTGATGCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                                                           <  4:238835/150‑1 (MQ=255)
                                                                                       tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                                <  2:953460/129‑1 (MQ=255)
                                                                                       tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                                >  1:953460/1‑129 (MQ=255)
                                                                                       tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                                                           >  3:124622/1‑150 (MQ=255)
                                                                                       tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                                                           >  2:230981/1‑150 (MQ=255)
                                                                                        ggTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGt                                                          <  1:480487/150‑1 (MQ=255)
                                                                                             tAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                     <  3:237490/150‑1 (MQ=255)
                                                                                             tAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                     >  3:378698/1‑150 (MQ=255)
                                                                                               aGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                                                   <  1:143659/150‑1 (MQ=255)
                                                                                                gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                                              <  1:95174/138‑1 (MQ=255)
                                                                                                gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                                              >  2:95174/1‑138 (MQ=255)
                                                                                                gCTGCAAGCGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag                                                  <  2:562256/150‑1 (MQ=255)
                                                                                                 cTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGt                                                 >  3:25205/1‑150 (MQ=255)
                                                                                                  tGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAGACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTa                                                >  3:542569/1‑150 (MQ=255)
                                                                                                        ggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAcc                                          >  3:492213/1‑150 (MQ=255)
                                                                                                        ggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAcc                                          <  3:656866/150‑1 (MQ=255)
                                                                                                          cAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat                                        >  1:822160/1‑150 (MQ=255)
                                                                                                          cAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat                                        <  2:735420/150‑1 (MQ=255)
                                                                                                               cGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATg                                    <  4:397815/150‑1 (MQ=255)
                                                                                                                gAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGt                                  >  2:917127/1‑150 (MQ=255)
                                                                                                                    cTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagca                              <  4:568375/150‑1 (MQ=255)
                                                                                                                         gATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcg                         <  4:430613/150‑1 (MQ=255)
                                                                                                                          aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                        <  4:144599/150‑1 (MQ=255)
                                                                                                                          aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                        >  3:341627/1‑150 (MQ=255)
                                                                                                                           tCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCa                       <  2:255586/150‑1 (MQ=255)
                                                                                                                              gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                                                   <  1:507679/119‑1 (MQ=255)
                                                                                                                              gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                                                   >  2:507679/1‑119 (MQ=255)
                                                                                                                              gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa                    >  4:542130/1‑150 (MQ=255)
                                                                                                                                    gATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACtt               >  2:820281/1‑149 (MQ=255)
                                                                                                                                    gATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACtt               <  1:820281/149‑1 (MQ=255)
                                                                                                                                      tCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCACTAACTACCATAATGTTAGCAACGCTTAGGATACTTCcg            <  2:41687/150‑1 (MQ=255)
                                                                                                                                        cTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgcg          <  2:881853/150‑1 (MQ=255)
                                                                                                                                         tGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCgcgc         >  3:318092/1‑150 (MQ=255)
                                                                                                                                         tGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                              >  2:531697/1‑113 (MQ=255)
                                                                                                                                         tGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                              <  1:531697/113‑1 (MQ=255)
                                                                                                                                         tGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCgcgc         >  2:941333/1‑150 (MQ=255)
                                                                                                                                           ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa       <  3:477581/150‑1 (MQ=255)
                                                                                                                                               aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca                                                                                   >  2:262969/1‑71 (MQ=255)
                                                                                                                                               aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca                                                                                   <  1:262969/71‑1 (MQ=255)
                                                                                                                                                 aaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGa  >  2:364440/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGA  >  NC_000913_3_pae_pgi/4399544‑4399837

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: