Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,182,771 G→C 11.4% A20G (GCC→GGC)  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,7710GC11.4% 232.0 / 16.6 70A20G (GCC→GGC) potDspermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (4/4);  total (33/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 = 1.00e+00

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                                                                                                                                >  2:315142/1‑146 (MQ=255)
             ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCag                                                                                                                                <  1:315142/146‑1 (MQ=255)
                    cGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                                             <  4:523758/126‑1 (MQ=255)
                    cGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCa                                                                                                                                             >  3:523758/1‑126 (MQ=255)
                        aGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCcgcg                                                                                                                   >  4:111305/1‑148 (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                                                                                                                     <  2:129464/136‑1 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgccg                                                                                                                     >  1:129464/1‑136 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                                             >  3:1027750/1‑144 (MQ=255)
                                  aTAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAgg                                                                                                             <  4:1027750/144‑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                                                                                                                                        >  1:265691/1‑109 (MQ=255)
                                         tAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccg                                                                                                                                        <  2:265691/110‑2 (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)
                                                       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                                                                                                                                               >  4:975641/1‑73 (MQ=255)
                                                                       gCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCt                                                                                                                                               <  3:975641/73‑1 (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                                                                                                                                         <  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:288316/1‑125 (MQ=255)
                                                                                                       tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTa                                                           >  3:535931/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                                                                                <  1:194225/101‑1 (MQ=255)
                                                                                                          aaaTACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTcc                                                                                >  2:194225/1‑101 (MQ=255)
                                                                                                            aTACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa                               <  1:402435/148‑1 (MQ=255)
                                                                                                                          tgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa                               >  3:344345/1‑134 (MQ=255)
                                                                                                                          tgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATaa                               <  4:344345/134‑1 (MQ=255)
                                                                                                                             cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCa                                                        <  1:717200/106‑1 (MQ=255)
                                                                                                                             cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCa                                                        >  2:717200/1‑106 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      <  2:429253/121‑1 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      <  2:958663/121‑1 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      >  1:429253/1‑121 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggctg                                      >  1:958663/1‑121 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg                     >  3:1092994/1‑138 (MQ=255)
                                                                                                                                cGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAg                     <  4:1092994/138‑1 (MQ=255)
                                                                                                                                  gcatGAGCGCCGCTCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                             >  4:268175/4‑96 (MQ=255)
                                                                                                                                  gcatGAGCGCCGCTCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                             <  3:268175/93‑1 (MQ=255)
                                                                                                                                      gAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAat     >  3:1029941/1‑148 (MQ=255)
                                                                                                                                         cggcgCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATg  <  4:1029941/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATG  >  NC_000913/1182632‑1182916

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