Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,797 T→G 6.4% P338P (CCT→CCG) ‡ potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,7970TG 6.4% 355.1 / 11.7 106P338P (CCT→CCG) ‡potFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (49/50);  new base G (4/3);  total (53/53)
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.99e-01

Only 100 of 107 total aligned reads displayed.

AGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCG  >  NC_000913/894650‑894945
                                                                                                                                                   |                                                                                                                                                    
aGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTg                                                                                                                                                     >  2:660768/1‑149 (MQ=255)
  aCGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCGGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCg                                                                                                                                                   >  4:578552/1‑149 (MQ=255)
  aCGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTACCCGCCTGCg                                                                                                                                                   >  4:587318/1‑149 (MQ=255)
    gAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGa                                                                                                                                                 >  2:342357/1‑149 (MQ=255)
    gAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGa                                                                                                                                                 >  4:391707/1‑149 (MQ=255)
     aaGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGAt                                                                                                                                                <  3:132220/149‑1 (MQ=255)
         cTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTc                                                                                                                                            <  2:595166/149‑1 (MQ=255)
         cTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTc                                                                                                                                            >  2:550394/1‑149 (MQ=255)
         cTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTc                                                                                                                                            >  3:369532/1‑149 (MQ=255)
             cAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGCGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGc                                                                                                                                        <  1:603100/149‑1 (MQ=255)
                  ccTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGc                                                                                                                                   >  3:851900/1‑149 (MQ=255)
                    tGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTg                                                                                                                                 >  4:429570/1‑149 (MQ=255)
                    tGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTg                                                                                                                                 >  3:864385/1‑149 (MQ=255)
                      aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGtt                                                                                                                               >  3:76995/1‑149 (MQ=255)
                      aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaa                                                                                                                                              <  3:143361/134‑3 (MQ=255)
                      aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaa                                                                                                                                              >  4:143361/1‑132 (MQ=255)
                        ttACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCa                                                                                                                             >  1:828378/1‑149 (MQ=255)
                           ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                          >  1:173790/1‑149 (MQ=255)
                           ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                          <  4:76995/149‑1 (MQ=255)
                           ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                          <  1:610663/149‑1 (MQ=255)
                               ctgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaa                                                                                                                      >  3:886140/1‑149 (MQ=255)
                                tgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                                                                                                     <  1:187489/149‑1 (MQ=255)
                                tgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                                                                                                     <  3:863317/149‑1 (MQ=255)
                                     ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                                <  2:619193/149‑1 (MQ=255)
                                     ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                                <  4:202241/149‑1 (MQ=255)
                                     ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                                >  3:37117/1‑149 (MQ=255)
                                      cGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAg                                                                                                               <  1:834386/149‑1 (MQ=255)
                                        gATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGa                                                                                                             <  4:369532/149‑1 (MQ=255)
                                         aTGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGAt                                                                                                            >  4:226933/1‑149 (MQ=255)
                                          tGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAGGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATc                                                                                                           <  4:303502/149‑1 (MQ=255)
                                              gtagCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGaa                                                                                                       >  1:304215/1‑149 (MQ=255)
                                                   gcATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTg                                                                                                                        >  2:316374/1‑127 (MQ=255)
                                                   gcATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTg                                                                                                                        <  1:316374/127‑1 (MQ=255)
                                                   gcATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCg                                                                                                  >  4:679839/1‑149 (MQ=255)
                                                    cATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGa                                                                                                 <  2:73963/149‑1 (MQ=255)
                                                     atatTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGAc                                                                                                <  4:886140/149‑1 (MQ=255)
                                                        tttCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgt                                                                                             <  1:342357/149‑1 (MQ=255)
                                                          tCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGTTGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtgt                                                                                           <  1:436844/149‑1 (MQ=255)
                                                          tCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtgt                                                                                           <  3:429570/149‑1 (MQ=255)
                                                          tCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtgt                                                                                           <  2:315354/149‑1 (MQ=255)
                                                          tCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGCGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCgtgt                                                                                           <  1:687343/149‑1 (MQ=255)
                                                            cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                          <  3:261191/116‑1 (MQ=255)
                                                            cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                          >  4:261191/1‑116 (MQ=255)
                                                            cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                         >  4:522614/1‑149 (MQ=255)
                                                            cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                         <  2:25542/149‑1 (MQ=255)
                                                            cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                         <  4:636021/149‑1 (MQ=255)
                                                             gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgcg                                                                                        >  1:678757/1‑149 (MQ=255)
                                                             gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaaa                                                                                                                                             >  3:764240/1‑93 (MQ=255)
                                                             gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaaa                                                                                                                                             <  4:764240/96‑4 (MQ=255)
                                                                     gtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCcgcg                                                                                >  1:11085/1‑149 (MQ=255)
                                                                         tATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTg                                                                            >  1:670186/1‑149 (MQ=255)
                                                                          aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg                                                                           >  1:709862/1‑149 (MQ=255)
                                                                          aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg                                                                           >  3:353341/1‑149 (MQ=255)
                                                                           tGCCAACGCCAATAAAGCAGCCACGCCGCTGGTCAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGa                                                                          <  4:813207/149‑1 (MQ=255)
                                                                              cAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCa                                                                       >  3:388324/1‑149 (MQ=255)
                                                                                aCGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCaaa                                                                     <  2:604362/149‑1 (MQ=255)
                                                                                aCGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCaaa                                                                     <  2:640557/149‑1 (MQ=255)
                                                                                    cAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGa                                                                 >  1:339616/1‑149 (MQ=255)
                                                                                     aaTAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGaa                                                                >  2:789436/1‑149 (MQ=255)
                                                                                     aaTAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGaa                                                                >  1:439819/1‑149 (MQ=255)
                                                                                      aTAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAag                                                               >  4:148684/1‑149 (MQ=255)
                                                                                        aaaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAagag                                                             <  2:452678/149‑1 (MQ=255)
                                                                                         aaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                            <  4:702236/149‑1 (MQ=255)
                                                                                         aaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                            <  2:173790/149‑1 (MQ=255)
                                                                                         aaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                            <  3:139663/149‑1 (MQ=255)
                                                                                          agcagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCg                                                           <  3:679839/149‑1 (MQ=255)
                                                                                          agcagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCg                                                           >  4:10553/1‑149 (MQ=255)
                                                                                            cagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGa                                                         <  2:897103/149‑1 (MQ=255)
                                                                                              gcCACGCCGCTGGTGAGTGCGGAAGCCCGTGCGAACCCAGGTATTTCTCCGCCTGCGGATGTTCGTGCGAAGCTGTTCCCTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGaaa                                                       <  4:37117/149‑1 (MQ=255)
                                                                                               ccACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGaaaa                                                      <  1:789436/149‑1 (MQ=255)
                                                                                                 aCGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAaata                                                    >  1:459947/1‑149 (MQ=255)
                                                                                                     cgcTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATcc                                                >  2:50148/1‑149 (MQ=255)
                                                                                                           tGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAaga                                                              >  4:371798/1‑129 (MQ=255)
                                                                                                           tGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAaga                                                              <  3:371798/129‑1 (MQ=255)
                                                                                                             aGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGt                                        >  1:734118/1‑149 (MQ=255)
                                                                                                                 cGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGAt                                    <  4:864385/149‑1 (MQ=255)
                                                                                                                       tCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGa                              <  3:226933/149‑1 (MQ=255)
                                                                                                                        ccGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAg                             <  1:50148/149‑1 (MQ=255)
                                                                                                                         cGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAgg                            <  2:761776/149‑1 (MQ=255)
                                                                                                                          gTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAggg                           <  1:660768/149‑1 (MQ=255)
                                                                                                                          gTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAggg                           <  1:372364/149‑1 (MQ=255)
                                                                                                                             agaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCAc                                                                                     >  3:473853/1‑88 (MQ=255)
                                                                                                                             agaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCAc                                                                                     <  4:473853/88‑1 (MQ=255)
                                                                                                                               aaCCCAGGTAGTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGGAAGTGCAGGATCCGAAAATCGACCGTGGGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAAGCCGCAGTCGTAGATGCCGGAGGGGCGca                      <  4:353341/149‑1 (MQ=255)
                                                                                                                                  ccAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacca                   >  3:175964/1‑149 (MQ=255)
                                                                                                                                  ccAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacca                   >  1:342130/1‑149 (MQ=255)
                                                                                                                                   cAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcaccac                  >  4:921058/1‑149 (MQ=255)
                                                                                                                                   cAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcaccac                  >  1:572334/1‑149 (MQ=255)
                                                                                                                                    aGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAACTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCTTAGATGCCGGAGGGGCGCACcaca                 <  3:578552/149‑1 (MQ=255)
                                                                                                                                      gTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACAcc               >  3:246935/1‑149 (MQ=255)
                                                                                                                                      gTATTTACCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACAcc               <  3:587318/149‑1 (MQ=255)
                                                                                                                                        gaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCaaa                                                                     <  4:330834/89‑1 (MQ=255)
                                                                                                                                        gaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCaaa                                                                     >  3:330834/5‑93 (MQ=255)
                                                                                                                                          gcaTCCGCAGGCGCATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCGCCAGCGCGTGGACCAAATTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACccgcc           >  4:598967/3‑149 (MQ=255)
                                                                                                                                            aTCCGCCTGCGGGTGTTCGTGCGAAGCTGTTCCCTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg         <  1:502887/149‑1 (MQ=255)
                                                                                                                                            aTCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg         >  4:32578/1‑149 (MQ=255)
                                                                                                                                             tCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGc        >  3:576087/1‑149 (MQ=255)
                                                                                                                                                gCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAt     <  2:590499/149‑1 (MQ=255)
                                                                                                                                                gCCTGCGGATGTTCGTGCGAAGCTGTTCAATCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAt     >  1:743180/1‑149 (MQ=255)
                                                                                                                                                   tGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCg  >  3:875483/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCG  >  NC_000913/894650‑894945

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