Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,182,764 G→A 7.6% H22H (CAC→CAT potD ← spermidine/putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,182,7640GA 7.6% 256.3 / 7.6 79H22H (CAC→CATpotDspermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base A (3/3);  total (38/41)
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.98e-01

TTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAG  >  NC_000913/1182619‑1182897
                                                                                                                                                 |                                                                                                                                     
ttttCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTg                                                                                                                                     <  2:404947/148‑1 (MQ=255)
 tttCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGa                                                                                                                                    >  2:191188/1‑148 (MQ=255)
 tttCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGa                                                                                                                                    >  2:783943/1‑148 (MQ=255)
   tCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGAGAgc                                                                                                                                  >  1:609711/1‑148 (MQ=255)
             gTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                        <  3:483130/148‑1 (MQ=255)
                     tttCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGc                                                                                                                <  2:117637/148‑1 (MQ=255)
                        cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                   >  2:132473/1‑142 (MQ=255)
                        cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                   <  1:132473/142‑1 (MQ=255)
                        cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                             <  1:191188/148‑1 (MQ=255)
                        cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                             <  1:783943/148‑1 (MQ=255)
                        cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGAGAGCGGCGCTCATGCCCAGTGCCag                                                                                                             <  2:609711/148‑1 (MQ=255)
                          ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                             >  2:315142/1‑146 (MQ=255)
                          ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                             <  1:315142/146‑1 (MQ=255)
                              aCTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgc                                                                                                                                  >  2:527979/1‑121 (MQ=255)
                              aCTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgc                                                                                                                                  <  1:527979/121‑1 (MQ=255)
                                 cGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                          <  4:523758/126‑1 (MQ=255)
                                 cGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                          >  3:523758/1‑126 (MQ=255)
                                     aGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCcgcg                                                                                                >  4:111305/1‑148 (MQ=255)
                                        cGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAg                                                                                                                                   >  3:363262/1‑110 (MQ=255)
                                        cGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAg                                                                                                                                   <  4:363262/110‑1 (MQ=255)
                                        cGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCacg                                                                                                                        <  2:713454/121‑3 (MQ=255)
                                        cGAATAGATAACCTTAATACCGGTTCCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCacg                                                                                                                        >  1:713454/1‑119 (MQ=255)
                                          aaTAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCacg                                                                                                                        <  3:509184/119‑3 (MQ=255)
                                          aaTAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCacg                                                                                                                        >  4:509184/1‑117 (MQ=255)
                                              gATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTgg                                                                                       <  3:111305/148‑1 (MQ=255)
                                               aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgccg                                                                                                  >  1:129464/1‑136 (MQ=255)
                                               aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgccg                                                                                                  <  2:129464/136‑1 (MQ=255)
                                               aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                          >  3:1027750/1‑144 (MQ=255)
                                               aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                          <  4:1027750/144‑1 (MQ=255)
                                                     ttAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTg                                                                                                                                     >  1:1043566/1‑95 (MQ=255)
                                                     ttAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTg                                                                                                                                     <  2:1043566/95‑1 (MQ=255)
                                                     ttAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGAcc                                                                                >  3:450471/1‑148 (MQ=255)
                                                      tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGa                                                                                  >  3:753379/1‑145 (MQ=255)
                                                      tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGa                                                                                  <  4:753379/145‑1 (MQ=255)
                                                      tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccg                                                                                                                     <  2:265691/110‑2 (MQ=255)
                                                      tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccg                                                                                                                     >  1:265691/1‑109 (MQ=255)
                                                          aCCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCAtttt                                                                           >  1:1122210/1‑148 (MQ=255)
                                                            cGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCAtttttt                                                                         <  4:450471/148‑1 (MQ=255)
                                                             ggTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgcggcg                                                                                                                              >  4:948152/1‑94 (MQ=255)
                                                             ggTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgcggcg                                                                                                                              <  3:948152/94‑1 (MQ=255)
                                                              gTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCg                                                                                                                                       >  2:103082/1‑84 (MQ=255)
                                                              gTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCg                                                                                                                                       <  1:103082/84‑1 (MQ=255)
                                                                    ttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                        >  1:1092622/1‑93 (MQ=255)
                                                                    ttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                        <  2:1092622/93‑1 (MQ=255)
                                                                         gAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                          <  1:141138/86‑1 (MQ=255)
                                                                         gAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                          >  2:141138/1‑86 (MQ=255)
                                                                           aCTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCt                                                          >  2:428529/1‑148 (MQ=255)
                                                                                   aGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                  <  2:395517/84‑1 (MQ=255)
                                                                                   aGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                  >  1:395517/1‑84 (MQ=255)
                                                                                    gCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCt                                                                                                                            <  3:975641/73‑1 (MQ=255)
                                                                                    gCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCt                                                                                                                            >  4:975641/1‑73 (MQ=255)
                                                                                     cAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCg                                                                                                                                       <  4:983583/61‑1 (MQ=255)
                                                                                     cAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCg                                                                                                                                       >  3:983583/1‑61 (MQ=255)
                                                                                        tCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                             >  3:1031418/1‑148 (MQ=255)
                                                                                        tCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                             <  1:428529/148‑1 (MQ=255)
                                                                                            ggcggcACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                   <  2:33794/74‑1 (MQ=255)
                                                                                            ggcggcACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                   >  1:33794/1‑74 (MQ=255)
                                                                                                     tACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGc                                <  4:1031418/148‑1 (MQ=255)
                                                                                                      aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGcc                                                                                                                      >  4:35875/1‑61 (MQ=255)
                                                                                                      aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGcc                                                                                                                      <  3:35875/61‑1 (MQ=255)
                                                                                                           gTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAg                          <  2:1122210/148‑1 (MQ=255)
                                                                                                                  tGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                   >  3:182296/1‑132 (MQ=255)
                                                                                                                  tGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                   <  4:182296/132‑1 (MQ=255)
                                                                                                                    tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                        >  3:535931/1‑125 (MQ=255)
                                                                                                                    tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                        >  3:288316/1‑125 (MQ=255)
                                                                                                                    tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                        <  4:288316/125‑1 (MQ=255)
                                                                                                                    tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                        <  4:535931/125‑1 (MQ=255)
                                                                                                                     aGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTaa                                                                  >  3:776404/1‑98 (MQ=255)
                                                                                                                     aGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTaa                                                                  <  4:776404/98‑1 (MQ=255)
                                                                                                                      gAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCaa               >  2:402435/1‑148 (MQ=255)
                                                                                                                       aaaTACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTcc                                                             >  2:194225/1‑101 (MQ=255)
                                                                                                                       aaaTACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTcc                                                             <  1:194225/101‑1 (MQ=255)
                                                                                                                         aTACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa            <  1:402435/148‑1 (MQ=255)
                                                                                                                                       tgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa            <  4:344345/134‑1 (MQ=255)
                                                                                                                                       tgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa            >  3:344345/1‑134 (MQ=255)
                                                                                                                                          cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCa                                     <  1:717200/106‑1 (MQ=255)
                                                                                                                                          cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCa                                     >  2:717200/1‑106 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                   >  1:958663/1‑121 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                   <  2:429253/121‑1 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                   <  2:958663/121‑1 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                   >  1:429253/1‑121 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg  <  4:1092994/138‑1 (MQ=255)
                                                                                                                                             cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg  >  3:1092994/1‑138 (MQ=255)
                                                                                                                                                 |                                                                                                                                     
TTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAG  >  NC_000913/1182619‑1182897

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: