Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,182,778 T→C 5.8% M18V (ATG→GTG)  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,7780TC 5.8% 238.3 / 3.6 69M18V (ATG→GTG) potDspermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATG  >  NC_000913/1182632‑1182916
                                                                                                                                                  |                                                                                                                                          
gTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                                           <  3:483130/148‑1 (MQ=255)
        tttCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGc                                                                                                                                   <  2:117637/148‑1 (MQ=255)
           cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                                      <  1:132473/142‑1 (MQ=255)
           cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                                      >  2:132473/1‑142 (MQ=255)
           cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                <  1:191188/148‑1 (MQ=255)
           cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                <  1:783943/148‑1 (MQ=255)
           cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGAGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                <  2:609711/148‑1 (MQ=255)
             ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                <  1:315142/146‑1 (MQ=255)
             ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                >  2:315142/1‑146 (MQ=255)
                        aGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCcgcg                                                                                                                   >  4:111305/1‑148 (MQ=255)
                                 gATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTgg                                                                                                          <  3:111305/148‑1 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgccg                                                                                                                     <  2:129464/136‑1 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgccg                                                                                                                     >  1:129464/1‑136 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                                             <  4:1027750/144‑1 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                                             >  3:1027750/1‑144 (MQ=255)
                                        ttAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGAcc                                                                                                   >  3:450471/1‑148 (MQ=255)
                                         tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGa                                                                                                     <  4:753379/145‑1 (MQ=255)
                                         tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGa                                                                                                     >  3:753379/1‑145 (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)
                                                       ttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                                           >  1:1092622/1‑93 (MQ=255)
                                                       ttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATg                                                                                                                                           <  2:1092622/93‑1 (MQ=255)
                                                              aCTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCt                                                                             >  2:428529/1‑148 (MQ=255)
                                                                      aGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                                     >  1:395517/1‑84 (MQ=255)
                                                                      aGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                                     <  2:395517/84‑1 (MQ=255)
                                                                           tCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                                                <  1:428529/148‑1 (MQ=255)
                                                                           tCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                                                >  3:1031418/1‑148 (MQ=255)
                                                                               ggcggcACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                                      <  2:33794/74‑1 (MQ=255)
                                                                               ggcggcACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAg                                                                                                                                      >  1:33794/1‑74 (MQ=255)
                                                                                        tACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGc                                                   <  4:1031418/148‑1 (MQ=255)
                                                                                         aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGcc                                                                                                                                         <  3:35875/61‑1 (MQ=255)
                                                                                         aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGcc                                                                                                                                         >  4:35875/1‑61 (MQ=255)
                                                                                              gTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAg                                             <  2:1122210/148‑1 (MQ=255)
                                                                                                     tGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                                      <  4:182296/132‑1 (MQ=255)
                                                                                                     tGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                                      >  3:182296/1‑132 (MQ=255)
                                                                                                       tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           <  4:535931/125‑1 (MQ=255)
                                                                                                       tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           >  3:535931/1‑125 (MQ=255)
                                                                                                       tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           >  3:288316/1‑125 (MQ=255)
                                                                                                       tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           <  4:288316/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                                                        >  2:717200/1‑106 (MQ=255)
                                                                                                                             cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCa                                                        <  1:717200/106‑1 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      >  1:958663/1‑121 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      <  2:429253/121‑1 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      >  1:429253/1‑121 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      <  2:958663/121‑1 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg                     >  3:1092994/1‑138 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg                     <  4:1092994/138‑1 (MQ=255)
                                                                                                                                  gcatGAGCGCCGCTCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                             <  3:268175/93‑1 (MQ=255)
                                                                                                                                  gcatGAGCGCCGCTCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                             >  4:268175/4‑96 (MQ=255)
                                                                                                                                      gAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAat     >  3:1029941/1‑148 (MQ=255)
                                                                                                                                         cggcgCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATg  <  4:1029941/148‑1 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           <  4:860829/84‑1 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           >  3:860829/1‑84 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                                      <  2:774607/89‑1 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                                      >  1:774607/1‑89 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAg                                             >  1:605400/1‑98 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCAt                           >  4:987768/1‑116 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCAt                           <  3:987768/116‑1 (MQ=255)
                                                                                                                                                cATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGAAATAAg                                             <  2:605400/98‑1 (MQ=255)
                                                                                                                                                  tGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           >  1:557055/1‑82 (MQ=255)
                                                                                                                                                  tGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           <  2:557055/82‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATG  >  NC_000913/1182632‑1182916

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