Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,796 C→A 6.4% P338H (CCT→CAT) ‡ 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,7960CA 6.4% 369.1 / 9.5 106P338H (CCT→CAT) ‡potFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (49/50);  new base A (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.95e-01

Only 100 of 107 total aligned reads displayed.

AAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAT  >  NC_000913/894648‑894942
                                                                                                                                                    |                                                                                                                                                  
aaaGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGcc                                                                                                                                                    >  3:203083/1‑149 (MQ=255)
  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                                                                                                                              >  3:864385/1‑149 (MQ=255)
                      tGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTg                                                                                                                              >  4:429570/1‑149 (MQ=255)
                        aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGtt                                                                                                                            >  3:76995/1‑149 (MQ=255)
                        aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaa                                                                                                                                           >  4:143361/1‑132 (MQ=255)
                        aaTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaa                                                                                                                                           <  3:143361/134‑3 (MQ=255)
                          ttACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCa                                                                                                                          >  1:828378/1‑149 (MQ=255)
                             ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       <  4:76995/149‑1 (MQ=255)
                             ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       <  1:610663/149‑1 (MQ=255)
                             ccTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       >  1:173790/1‑149 (MQ=255)
                                 ctgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaa                                                                                                                   >  3:886140/1‑149 (MQ=255)
                                  tgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                                                                                                  <  3:863317/149‑1 (MQ=255)
                                  tgcGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGaaa                                                                                                                  <  1:187489/149‑1 (MQ=255)
                                       ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                             <  2:619193/149‑1 (MQ=255)
                                       ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                             >  3:37117/1‑149 (MQ=255)
                                       ccGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCa                                                                                                             <  4:202241/149‑1 (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                                                                                                                       >  4:261191/1‑116 (MQ=255)
                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCActc                                                                                                                       <  3:261191/116‑1 (MQ=255)
                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                      <  2:25542/149‑1 (MQ=255)
                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                      <  4:636021/149‑1 (MQ=255)
                                                              cGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgc                                                                                      >  4:522614/1‑149 (MQ=255)
                                                               gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTgcg                                                                                     >  1:678757/1‑149 (MQ=255)
                                                               gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaaa                                                                                                                                          <  4:764240/96‑4 (MQ=255)
                                                               gACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCAGGCGGAtaaa                                                                                                                                          >  3:764240/1‑93 (MQ=255)
                                                                       gtTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCcgcg                                                                             >  1:11085/1‑149 (MQ=255)
                                                                           tATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTg                                                                         >  1:670186/1‑149 (MQ=255)
                                                                            aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg                                                                        >  3:353341/1‑149 (MQ=255)
                                                                            aTGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTgg                                                                        >  1:709862/1‑149 (MQ=255)
                                                                             tGCCAACGCCAATAAAGCAGCCACGCCGCTGGTCAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGa                                                                       <  4:813207/149‑1 (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                                                         <  3:139663/149‑1 (MQ=255)
                                                                                           aaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                         <  4:702236/149‑1 (MQ=255)
                                                                                           aaGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGc                                                         <  2:173790/149‑1 (MQ=255)
                                                                                            agcagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCg                                                        >  4:10553/1‑149 (MQ=255)
                                                                                            agcagcCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCg                                                        <  3:679839/149‑1 (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                                                           <  3:371798/129‑1 (MQ=255)
                                                                                                             tGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAaga                                                           >  4:371798/1‑129 (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:372364/149‑1 (MQ=255)
                                                                                                                            gTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAggg                        <  1:660768/149‑1 (MQ=255)
                                                                                                                              gagaACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGgc                      <  3:391707/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                >  1:342130/1‑149 (MQ=255)
                                                                                                                                    ccAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGcacca                >  3:175964/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                                                                  >  3:330834/5‑93 (MQ=255)
                                                                                                                                          gaacaTCCGCAGGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCaaa                                                                  <  4:330834/89‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
AAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGAGAACCCAGGTATTTATCCGCCTGCGGATGTTCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAAT  >  NC_000913/894648‑894942

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: