Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,110 G→A 10.2% T26I (ACC→ATC)  lptE ← LPS assembly OM complex LptDE, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,1100GA10.2% 290.4 / 19.3 98T26I (ACC→ATC) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base A (4/6);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAG  >  NC_000913/671962‑672256
                                                                                                                                                    |                                                                                                                                                  
cAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTgg                                                                                                                                                    <  3:590318/149‑1 (MQ=255)
 aaGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGt                                                                                                                                                   >  3:244198/1‑149 (MQ=255)
 aaGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGt                                                                                                                                                   >  2:493744/1‑149 (MQ=255)
  aGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATa                                                                                                                                                  <  4:60210/149‑1 (MQ=255)
    gATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                >  3:10960/1‑149 (MQ=255)
        gAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGc                                                                                                                                            >  3:1587/1‑149 (MQ=255)
         aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                           <  4:81540/149‑1 (MQ=255)
         aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                           <  3:316356/149‑1 (MQ=255)
         aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCa                                                                                                                                           <  1:608037/149‑1 (MQ=255)
           cGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGa                                                                                                                                         <  4:10960/149‑1 (MQ=255)
           cGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAgg                                                                                                                                         <  2:693043/149‑2 (MQ=255)
           cGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAgg                                                                                                                                         <  2:560987/149‑2 (MQ=255)
            gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                        <  2:296939/149‑1 (MQ=255)
            gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGAAGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                        >  2:486675/1‑149 (MQ=255)
               cTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcc                                                                                                                                     >  2:292829/1‑149 (MQ=255)
                ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                    >  2:186291/1‑149 (MQ=255)
                         gTTTCTTTATCAAGCGACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                           >  4:673732/1‑149 (MQ=255)
                              tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                      <  1:486675/149‑1 (MQ=255)
                               ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                     <  1:493744/149‑1 (MQ=255)
                                 aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                   <  1:640335/149‑1 (MQ=255)
                                  tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagc                                                                                                                                  >  2:620835/1‑133 (MQ=255)
                                  tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagc                                                                                                                                  <  1:620835/133‑1 (MQ=255)
                                  tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                  <  4:126761/149‑1 (MQ=255)
                                   cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                 <  4:344042/149‑1 (MQ=255)
                                   cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                 <  4:181505/149‑1 (MQ=255)
                                    aaGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACATTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTa                                                                                                                >  4:618170/1‑149 (MQ=255)
                                        aaCTCGACACCATTCAGACGTAACTGGTTCCGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                            >  3:413203/1‑149 (MQ=255)
                                        aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                            >  3:407829/1‑149 (MQ=255)
                                        aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAacac                                                                                                            >  1:675014/1‑149 (MQ=255)
                                           tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGt                                                                                                                                                   <  3:440905/107‑1 (MQ=255)
                                           tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGt                                                                                                                                                   >  4:440905/1‑107 (MQ=255)
                                           tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAATCACCGc                                                                                                         <  4:573490/149‑1 (MQ=255)
                                           tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                         >  3:108852/1‑149 (MQ=255)
                                             gACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCa                                                                                                       <  1:45264/149‑1 (MQ=255)
                                              acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTGTTAATCACCGccgg                                                                                                      <  1:48637/149‑3 (MQ=255)
                                              acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCag                                                                                                      >  3:268895/1‑149 (MQ=255)
                                                acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                               <  1:11656/138‑1 (MQ=255)
                                                acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                               >  2:11656/1‑138 (MQ=255)
                                                  cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTa                                                                                                                                                  <  3:276463/101‑1 (MQ=255)
                                                  cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTa                                                                                                                                                  >  4:276463/1‑101 (MQ=255)
                                                  cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGAt                                                                                                  <  1:630378/149‑1 (MQ=255)
                                                   aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAggt                                                                                                                                        >  4:621178/1‑108 (MQ=255)
                                                   aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAggt                                                                                                                                        <  3:621178/110‑3 (MQ=255)
                                                    ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                <  3:509058/149‑1 (MQ=255)
                                                    ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                >  1:378146/1‑149 (MQ=255)
                                                     tCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaac                                                                                               <  4:1587/149‑1 (MQ=255)
                                                       aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaa                                                                                             <  1:353131/149‑1 (MQ=255)
                                                           gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                         <  1:404714/149‑1 (MQ=255)
                                                            tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATg                                                                                        >  4:506190/1‑149 (MQ=255)
                                                               cTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTg                                                                                     >  4:551131/1‑149 (MQ=255)
                                                                tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGACCCTGCGTGGTATCACGCAGATGCCAGCCACAACCGGCGGTGATTGACACCGCGAGAGATAACAACAATGTCGc                                                                                    >  4:84977/1‑149 (MQ=255)
                                                                   ttACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAg                                                                                 >  1:96007/1‑149 (MQ=255)
                                                                    tACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGa                                                                                >  2:478270/1‑149 (MQ=255)
                                                                       gcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGatat                                                                             >  3:548036/1‑149 (MQ=255)
                                                                          ccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGc                                                                          >  1:365252/1‑149 (MQ=255)
                                                                                 gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                   >  2:271416/1‑149 (MQ=255)
                                                                                   ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcct                                                                 >  3:313626/1‑149 (MQ=255)
                                                                                                cGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                 >  1:385339/1‑88 (MQ=255)
                                                                                                cGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                 <  2:385339/88‑1 (MQ=255)
                                                                                                 ggATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCa                                                   >  1:366183/1‑149 (MQ=255)
                                                                                                  gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAgaga                                                                                                   <  2:89823/100‑1 (MQ=255)
                                                                                                  gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAgaga                                                                                                   >  1:89823/1‑100 (MQ=255)
                                                                                                        cTGAGTCCAGGATCATTACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCTCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                            <  4:413203/149‑1 (MQ=255)
                                                                                                        cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAc                                                                                                                                              >  4:424349/1‑51 (MQ=255)
                                                                                                        cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAc                                                                                                                                              <  3:424349/51‑1 (MQ=255)
                                                                                                        cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                            <  2:675014/149‑1 (MQ=255)
                                                                                                            gTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAtt                                        <  2:439913/149‑1 (MQ=255)
                                                                                                              ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                      <  4:108852/149‑1 (MQ=255)
                                                                                                               cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                     <  4:313626/149‑1 (MQ=255)
                                                                                                               cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                     >  1:19155/1‑149 (MQ=255)
                                                                                                                aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                   <  1:83202/118‑1 (MQ=255)
                                                                                                                aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                   >  2:83202/1‑118 (MQ=255)
                                                                                                                aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                    >  2:645077/1‑149 (MQ=255)
                                                                                                                 ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                   >  3:69489/1‑149 (MQ=255)
                                                                                                                 ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                   <  3:143388/149‑1 (MQ=255)
                                                                                                                   aTCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGt                                 <  1:478270/149‑1 (MQ=255)
                                                                                                                    tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                <  1:645077/149‑1 (MQ=255)
                                                                                                                    tCATGACCTTCATAGCGGAAGTAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                >  2:432249/1‑149 (MQ=255)
                                                                                                                     cATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGttt                               >  1:549659/1‑149 (MQ=255)
                                                                                                                       tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                             <  2:378146/149‑1 (MQ=255)
                                                                                                                       tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAc                             <  2:96007/149‑1 (MQ=255)
                                                                                                                        gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                            >  3:514278/1‑149 (MQ=255)
                                                                                                                        gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                            <  4:69489/149‑1 (MQ=255)
                                                                                                                        gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                            >  3:387244/1‑149 (MQ=255)
                                                                                                                         aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                           >  2:173345/1‑149 (MQ=255)
                                                                                                                           cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                         <  1:186291/149‑1 (MQ=255)
                                                                                                                           cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                         >  3:27452/1‑149 (MQ=255)
                                                                                                                             tCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTa                       <  1:173345/149‑1 (MQ=255)
                                                                                                                               aTAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACg                     >  1:47046/1‑149 (MQ=255)
                                                                                                                                  gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaa                  <  2:156684/149‑1 (MQ=255)
                                                                                                                                       aggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCAc             <  1:292829/149‑1 (MQ=255)
                                                                                                                                        ggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt            <  4:407829/149‑1 (MQ=255)
                                                                                                                                        ggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt            <  2:19155/149‑1 (MQ=255)
                                                                                                                                           atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa         >  2:476497/3‑149 (MQ=255)
                                                                                                                                            tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt            >  2:597653/2‑145 (MQ=255)
                                                                                                                                            tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACt            <  1:597653/144‑1 (MQ=255)
                                                                                                                                            tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTAc        >  1:100174/2‑149 (MQ=255)
                                                                                                                                                cGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc    <  3:673732/149‑1 (MQ=255)
                                                                                                                                                  tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg  <  4:268895/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAG  >  NC_000913/671962‑672256

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: